#                      MAPLE fuer ANALYSIS 
#                RWTH -Aachen, Wintersemester 1997/98
#                         E. Goerlich
#
# Maple -Input von Kapitel 2:
#
# r-Zykloide:
> H:=(r,u)->[u*t-r*sin(u*t),1-r*cos(u*t),t=0..4*Pi];
# Kleiner r-Punkt:
> J:=(r,u)->[u*4*Pi-r*sin(u*4*Pi)+0.05*cos(t),1-r*cos(u*4*Pi)+0.05*sin(t
> ),t=0..2*Pi];
> with(plots):
> settime:=time();animate(\{[4*Pi*u+cos(t),1+sin(t),t=0..2*Pi],H(1,u),
> H(1.25,u),H(0.75,u),J(1,u),J(1.25,u),J(0.75,u)\},u=0..1,view=[-1..14,-
> 1..3], 
> scaling=constrained,frames=50);cpu_time:=(time()-settime)*seconds;
> settime:=time();animate(\{[4*Pi*u+cos(t),1+sin(t),t=0..2*Pi],H(1,u),H(
> 1.25,u),H(0.75,u),J(1,u),J(1.25,u),J(0.75,u),[4*Pi*u-t*1.25*sin(u*4*Pi
> ),1-t*1.25*cos(u*4*Pi),t=0..1]\},u=0..1,view=[-1..14,-1..3],scaling=co\
> nstrained,color=red,frames=50);cpu_time:=(time()-settime)*seconds;
# 
#           Bewegung auf einer Flaeche im Raum:
> a:= proc(x,y,phi) local result; if 
> (x-cos(phi))**2+(y-sin(phi))**2 >0.01 then     result:=(x*y)/(x*x+y*y)
> else result := FAIL fi; end;
> animate3d('a(x,y,phi)','x'=-1..1,'y'=-1..1,'phi'=0..2*Pi,frames=16,
> numpoints=1600,orientation=[109,49]);
                    
#                        =========================