{VERSION 4 0 "IBM INTEL NT" "4.0" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 1 }{CSTYLE "2D Math" -1 2 "Times" 0 1 0 0 0 0 0 0 2 0 0 0 0 0 0 1 }{CSTYLE "2D Output" 2 20 "" 0 1 0 0 255 1 0 0 0 0 0 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 "Warning" -1 7 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 2 2 2 2 2 1 1 1 1 1 }1 1 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Maple Output" -1 12 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }1 3 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "R3 Font 0" -1 256 1 {CSTYLE "" -1 -1 "Helvetica" 1 10 0 0 0 1 2 1 2 2 2 2 1 1 1 1 }1 1 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "R3 Fo nt 2" -1 257 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 0 1 2 1 2 2 2 2 1 1 1 1 }1 1 0 0 0 0 1 0 1 0 2 2 0 1 }} {SECT 0 {EXCHG {PARA 0 "" 0 "" {TEXT -1 94 "c International Thomson Pu blishing Bonn filename: procnewt.ms" }}{PARA 0 "" 0 "" {TEXT -1 102 "Autor: Komma \+ Datum: 4.5.94" }} {PARA 0 "" 0 "" {TEXT -1 7 "Index: " }}{PARA 0 "" 0 "" {TEXT -1 123 "T hema: Newtons Physik. Hier nur die procedure zur geschlossenen L\366su ng der Bewegungsgleichung (den .m-file mit read laden)." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 47 "restart;with(linalg):with(student): with(plots):" }}{PARA 7 "" 1 "" {TEXT -1 80 "Warning, the protected na mes norm and trace have been redefined and unprotected\n" }}{PARA 7 " " 1 "" {TEXT -1 50 "Warning, the name changecoords has been redefined \n" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 88 "hprocnewt:=TEXT(`FUNK TION: Berechnung der geschlossenen L\366sung der Bewegungsgleichung.`, " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "`AUFRUF: newton(Kraftgesetz);`, " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "`PARAMETER:`," }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 69 "`Kraftgesetz: Vektor der Kraftkomponenten als Fu nktion von x,v, ...`," }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "`ERGEBNIS: `," }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 65 "`Newton stellt die L\366sungs vektoren rf, vf und af zur Verf\374gung`);" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%*hprocnewtG-%%TEXTG6(%`oFUNKTION:~Berechnung~der~gesc hlossenen~L|azsung~der~Bewegungsgleichung.G%=AUFRUF:~newton(Kraftgeset z);G%+PARAMETER:G%]oKraftgesetz:~Vektor~der~Kraftkomponenten~als~Funkt ion~von~x,v,~...G%*ERGEBNIS:G%hnNewton~stellt~die~L|azsungsvektoren~rf ,~vf~und~af~zur~Verf|gzgungG" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 61 "newton:=proc(F) global x,y,z,r,v,a,xx,yy,zz,rf,vf,af,sol,sys;" } }{PARA 0 "> " 0 "" {MPLTEXT 1 0 96 "#with(linalg):with(student):with(p lots):\n# kein Aufruf in Prozeduren von Packages in Maple 6!!!!" }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "#unassign('x(t)','y(t)','z(t)');" } }{PARA 0 "> " 0 "" {MPLTEXT 1 0 122 "#unassign('x(t)','y(t)','z(t)','x 0','vx0','y0','vy0','z0','vz0'); #witzlos, weil schon beim Auruf Auswe rtung von F erfolgt" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "#print(x(t), F());" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 83 "#x(0):='x0': D(x)(0):='vx0 ': y(0):='y0': D(y)(0):='vy0':z(0):='z0': D(z)(0):='vz0';" }}{PARA 0 " > " 0 "" {MPLTEXT 1 0 28 "r:=vector([x(t),y(t),z(t)]);" }}{PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 17 "v:=map(diff,r,t);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "a:=map(diff,v,t);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "sys:=equate(m*a,F);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "#print(s ys);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 49 "sol:=dsolve(sys,\{x(t),y(t),z(t)\},method=laplace);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "#print(sol);" }}{PARA 0 "> " 0 " " {MPLTEXT 1 0 54 "if sol=NULL then sol:=dsolve(sys,\{x(t),y(t),z(t)\} ) fi;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 53 "if sol=NULL then RETURN(`k eine Loesung gefunden`) fi;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "assi gn(sol);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 67 "xx:=makeproc(x(t),t): \+ yy:=makeproc(y(t),t): zz:=makeproc(z(t),t):" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 86 "rf:=makeproc(map(eval,r),t); vf:=makeproc(map(eval,v) ,t); af:=makeproc(map(eval,a),t);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 82 "x(0):='x0': D(x)(0):='vx0': y(0):='y0': D(y)(0):='vy0':z(0):='z0': D(z)(0):='vz0';" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "RETURN(op(rf)); " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 4 "end;" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%'newtonGR6#%\"FG6\"F(F(C7>%\"rG-%'vectorG6#7%-%\"xG6# %\"tG-%\"yGF2-%\"zGF2>%\"vG-%$mapG6%%%diffGF+F3>%\"aG-F;6%F=F9F3>%$sys G-%'equateG6$*&%\"mG\"\"\"F?FI9$>%$solG-%'dsolveG6%FC<%F6F4F0/%'method G%(laplaceG@$/FL%%NULLG>FL-FN6$FCFP@$FU-%'RETURNG6#%7keine~Loesung~gef undenG-%'assignG6#FL>%#xxG-%)makeprocG6$F0F3>%#yyG-F_o6$F4F3>%#zzG-F_o 6$F6F3>%#rfG-F_o6$-F;6$%%evalGF+F3>%#vfG-F_o6$-F;6$F_pF9F3>%#afG-F_o6$ -F;6$F_pF?F3>-F16#\"\"!.%#x0G>--%\"DG6#F1F^q.%$vx0G>-F5F^q.%#y0G>--Feq 6#F5F^q.%$vy0G>-F7F^q.%#z0G>--Feq6#F7F^q.%$vz0G-Ffn6#-%#opG6#FjoF(60F1 F5F7F+F9F?F]oFboFfoFjoFapFgpFLFCF(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "save \"procnewtr6.m\";" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 25 "komma@oe.uni-tuebingen.de" }}}} {MARK "0 0 0" 22 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }