Code For Producing An Euler Spiral
The following Sage code produces the second graph above. The first four lines express the Euler spiral component. Fresnel functions could not be found. Instead, the integrals of two expanded Taylor series are adopted. The remaining code expresses respectively the tangent and the circle, including the computation for the center coordinates.
var('L') p = integral(taylor(cos(L^2), L, 0, 12), L) q = integral(taylor(sin(L^2), L, 0, 12), L) r1 = parametric_plot(, (L, 0, 1), color = 'red') r2 = line(, rgbcolor = 'blue') x1 = p.subs(L = 1) y1 = q.subs(L = 1) R = 0.5 x2 = x1 - R*sin(1.0) y2 = y1 + R*cos(1.0) r3 = circle((x2, y2), R, rgbcolor = 'green') show(r1 + r2 + r3, aspect_ratio = 1, axes=false)The following is Mathematica code for the Euler spiral component (it works directly in wolframalpha.com):
ParametricPlot[ {FresnelC] t]/Sqrt], FresnelS] t]/Sqrt]}, {t, -10, 10}]Read more about this topic: Euler Spiral
Famous quotes containing the words code, producing and/or spiral:
“Acknowledge your will and speak to us all, “This alone is what I will to be!” Hang your own penal code up above you: we want to be its enforcers!”
—Friedrich Nietzsche (1844–1900)
“The writer who aims at producing the platitudes which are “not for an age, but for all time” has his reward in being unreadable in all ages.... The man who writes about himself and his own time is the only sort of man who writes about all people and about all time.”
—George Bernard Shaw (1856–1950)
“Year after year beheld the silent toil
That spread his lustrous coil;
Still as the spiral grew,
He left the past year’s dwelling for the new,
Stole with soft step its shining archway through,
Built up its idle door,
Stretched in his last-found home, and knew the old no more.”
—Oliver Wendell Holmes, Sr. (1809–1894)