Cicloids e Trochoids

from math import *

from Vector import *
from Curve  import Curve

class Trochoid(Curve):
    Name="Trochoid"

    #Parameters
    a=1.0
    b=0.0

    NP=400
    t1=0
    t2=8.0*pi
    
    def __init__(self,a=1.0,b=0.0,NP=0):
        self.a=a
        self.b=b
        if (NP): self.NP=NP
        
        return
    
    def R(self,t):
        return Vector([
            self.a*t-(self.a+self.b)*sin(t),
            self.a-(self.a+self.b)*cos(t)
        ])

    def PMin(self):
        return Vector([self.t1,-1.5])
    
    def PMax(self):
        return Vector([self.t2,4.0 ])

from Vector   import *
from Trochoid import Trochoid

a=1.0

b1=-2.0
b2=2.0

N=100
NP=400

db=1.0*(b2-b1)/( 1.0*(N-1) )

cycloid=Trochoid(a,0.0,NP)
cycloid.Rs()

palette=[
    "red","yellow","oragne","blue","green","black",
    "magenta","cyan"
]

b=b1

files=[]

ssvg=[]
for n in range(N):
    curve=Trochoid(a,b,NP)
    curve.Rs()

    curve.Init_Canvas()

    #Transfer Canvas
    cycloid.Canvas=curve.Canvas

    fname="Trochoid/"+("%03d" % n)+".svg"
    options={
        "stroke": 'blue',
        "style": {
            "stroke-width": 2,
        }
    }

    svg=[]
    svg=svg+curve.Curve_SVG_Line(options)
    
    options={
        "stroke": palette[ (n % len(palette)) ],
        "style": {
            "stroke-width": 2,
        }
    }

    ssvg=ssvg+curve.Curve_SVG_Line(options)
    
    options[ "stroke" ]='green'
    svg=svg+cycloid.Curve_SVG_Line(options)

    curve.SVG_2_File(fname,svg)
    files.append(fname)

    b+=db
curve.Curve_Carousel(files)

curve.SVG_2_File("Trochoid/Trochoids.svg",ssvg)