{VERSION 6 0 "IBM INTEL NT" "6.0" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 1 }{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }1 1 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Headi ng 1" -1 3 1 {CSTYLE "" -1 -1 "Times" 1 18 0 0 0 1 2 1 2 2 2 2 1 1 1 1 }1 1 0 0 8 4 1 0 1 0 2 2 0 1 }} {SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "with(VectorCalculus):" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 48 "differentiaal-mechanica met vrijheidsbepe rkingen" }}{PARA 0 "" 0 "" {TEXT -1 10 "26-06-2006" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 27 "param:=[,u=0. .2*Pi];" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "Norm(diff(subs(u =u(t),param[1]),t),2)=1;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 10 "dsolve(%);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "diff(u(t)^3, t)^" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "indets(param[1],symb ol);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 63 "[x,y]=Jacobian(para m[1],convert(indets(param[1],symbol),list));" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "plot([u^3,u,u=0..2*Pi]);" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 87 " VectorCalculus[PathInt](1,[x,y]=Path( param[1], lhs(param[2])=lhs(rhs( param[2]))..t ));" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "value(% );" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 42 "int(,lhs(param[2])=lh s(rhs(param[2]))..t);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "par am[2];" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{SECT 1 {PARA 3 "" 0 "" {TEXT -1 10 "potentiaal" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "param:=[,u=-2..2];" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 11 "V:=p->p[2];" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 59 "addtime:=a->subs(map(i->i=i(t),indets(param[1],symbol)),a);" } }}{EXCHG }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 63 "vgl:=diff(u(t),t$2 )=-addtime(diff(V(param[1]),u))-diff(u(t),t);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 58 "sol:=dsolve(\{vgl,u(0)=lhs(rhs(param[2])),D(u)(0 )=0\},u(t));" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 42 "subs(sol,addtime(convert(param[1],l ist)));" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "plot([%[],t=0..3 ]);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 39 "map(i->i=i(t),indets(param[1],symbol));" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{SECT 1 {PARA 3 "" 0 "" {TEXT -1 2 "3d" }} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "with(plots):with(plottools): " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 53 "param:=[,u=0..2*Pi,v=-1..1];" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 89 "test:=proc(param,\{V:=(p->p[3]),friction:=1\})\n\nind:=convert (indets(param[1],symbol),list);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 49 " plotparam:=param->plot3d(param[1],param[2..3][]);" }}{PARA 0 "> " 0 " " {MPLTEXT 1 0 22 "plotparam(param);m:=%:" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 39 "addtime:=a->subs(map(i->i=i(t),ind),a);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "subs(param[2..3],u);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 66 "beginvwd:=map(i->[i(0)=lhs(subs(param[2..3],i)),D(i)( 0)=1][],ind);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 84 "vgl:=map(i->diff(i (t),t$2)=-addtime(diff(V(param[1]),i))-friction*diff(i(t),t),ind);" }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "stelsel:=\{vgl[],beginvwd[]\};" }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 56 "sol:=dsolve(\{vgl[],beginvwd[]\},ad dtime(ind));print(sol);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 54 "poi:=poi nt(subs(sol,addtime(convert(param[1],list))));" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 94 "display(spacecurve(subs(sol,addtime(convert(param[1], list))),t=0..2,thickness=4,color=red),m);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "end proc:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 96 "test(param,V=(p->-p[3]));test(param,V=(p->p[3]));test(param,V=(p-> 0),friction=0);display(%,%%%);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 66 "test([,phi=0..2*Pi,theta=0..2*P i]);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 95 "test([<(3+cos(theta))*sin(phi),(3+sin(theta))*sin(p hi),2+cos(phi)>,phi=0..2*Pi,theta=0..2*Pi]);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 57 "map(''evalf'',subs(sol,addtime(convert(param[1], list))));" }}}{EXCHG }{EXCHG }{EXCHG }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 114 "animatereal(display,[point(subs(sol,addtime(convert( param[1],list))),color=red,symbolsize=15,symbol=box)],t=0..1);" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 86 "animate(display,[point([-exp (t), -t, -t],color=red,symbolsize=15,symbol=box)],t=0..1);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}} {EXCHG }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "plot([%[],t=0..3]); " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 39 " map(i->i=i(t),indets(param[1],symbol));" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{SECT 0 {PARA 3 "" 0 "" {TEXT -1 5 "multi" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 57 "restart:with(plots):with(plottools):with(VectorCalculus):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 84 "param:=[,phi=0..Pi,theta=0..2*Pi];" }}}{EXCHG }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 41 "points:=[seq([t,param,t/5,t/5],t=1. .10)];" }}}{EXCHG }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 59 "ind:=para m->convert(indets(param,Or(symbol,indexed)),list);" }}{PARA 0 "> " 0 " " {MPLTEXT 1 0 74 "addlabel:=point->subs(map(i->i=i[point[1]],ind(poin t[2][1][1])),point[2]);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 54 "addtime: =param->subs(map(i->i=i(t),ind(param)),param);" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 39 "labels:=point->ind(addlabel(point)[1]);" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "pos:=point->addlabel(point)[ 1];" }}}{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG } {EXCHG }{EXCHG }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 58 "TotalPotenti al:=add(1/Norm(-pos(i),2),i in points):" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 83 "potential:=i->subs(Equate([x,y,z],pos(i)),Tota lPotential-1/Norm(-pos(i),2)):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "potential(points[1]):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "friction:=0.5;" }}}{EXCHG }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 122 "vgl:=map(point->map(i->diff(i(t),t$2)=-addtime(diff( potential(point),i))-friction*diff(i(t),t),labels(point))[],points)[]: " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 83 "beginvwd:=map(point->ma p(i->[i(0)=point[3],D(i)(0)=0][],labels(point))[],points)[];" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "sol:=dsolve([vgl,beginvwd],n umeric);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "m:=plot3d(point s[1][2][]):" }}}{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 173 "snapshot:=proc(t) local a; \na:=subs(sol (t)[2..-1],map(p->point(convert(addtime(pos(p)),list),symbol=box,color =red,symbolsize=15),points));\ndisplay(m,map(evalf,a));\nend proc:" }} }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 172 "pointshot:=proc(t) local a ; \na:=subs(sol(t)[2..-1],map(p->point(convert(addtime(pos(p)),list),s ymbol=box,color=red,symbolsize=15),points));\ndisplay(map(evalf,a));\n end proc:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "snapshot(0);snapshot(1);" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 44 "animate(pointshot,[t],t=0..10,background=m);" }}} {EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG } {EXCHG }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 85 "odeplot(sol,map(poin t->convert(addtime(pos(point)),list),points),t=0..50,frames=100);" }}} {EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG } {EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG } {EXCHG }}{EXCHG }{EXCHG }{EXCHG }{EXCHG }}{MARK "25" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }