{VERSION 5 0 "IBM INTEL LINUX" "5.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 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Text Output" -1 2 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 0 0 0 0 0 1 3 0 3 0 }1 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 2 6 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 2 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "W arning" 2 7 1 {CSTYLE "" -1 -1 "" 0 1 0 0 255 1 0 0 0 0 0 0 1 0 0 0 } 0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple Output" 0 11 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 3 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 11 12 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 }1 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple Plo t" 0 13 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }} {SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart;" }}{PARA 0 " > " 0 "" {MPLTEXT 1 0 16 "with( DEtools ):" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "with( plots ):" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 15 " with( linalg ):" }}{PARA 7 "" 1 "" {TEXT -1 50 "Warning, the name chan gecoords has been redefined\n" }}{PARA 7 "" 1 "" {TEXT -1 45 "Warning, the name adjoint has been redefined\n" }}{PARA 7 "" 1 "" {TEXT -1 80 "Warning, the protected names norm and trace have been redefined and u nprotected\n" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 81 "Setup the ODE. D on't worry about constants, we'll put those in in the next step:" }} {PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "ode:=diff(v(t),t)=-g+(k/m)*v(t)^2;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$odeG/-%%diffG6$-%\"vG6#%\"tGF,,&%\"gG!\"\"*(%\"kG\"\"\"%\"mGF /F)\"\"#F2" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 48 "substitute the cons tants with the \"subs\" routine" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 35 "ode:=subs(\{m=100,k=0.2,g=9. 8\},ode);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$odeG/-%%diffG6$-%\"vG6 #%\"tGF,,&$\"#)*!\"\"F0*&$\"+++++?!#7\"\"\")F)\"\"#F5F5" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 91 "we can also add the initial condition so \+ that we don't have to solve for the constant later" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "ic:=v(0)=0; " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#icG/-%\"vG6#\"\"!F)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "ode;\nic;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%%diffG6$-%\"vG6#%\"tGF*,&$\"#)*!\"\"F.*&$\"+++++?!#7 \"\"\")F'\"\"#F3F3" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"vG6#\"\"!F' " }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 63 "Turn the verbosity up to 3 so we can see what \"dsolve\" is doing" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "infolevel[dsolve]:=3;" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>&%*infolevelG6#%'dsolveG\"\"$" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 79 "Check the odeadvisor. We know it \+ is separable, but it's a nice command to know" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "odeadvisor( ode);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7#%,_quadratureG" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "odeadvisor(ode,[separable]);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#7#%+_separableG" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 48 "Solve....and tell it to use separable methods..." }} }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 27 "sol:=dsolve(\{ode,ic\},v(t) );" }}{PARA 6 "" 1 "" {TEXT -1 29 "Methods for first order ODEs:" }} {PARA 6 "" 1 "" {TEXT -1 37 "--- Trying classification methods ---" }} {PARA 6 "" 1 "" {TEXT -1 19 "trying a quadrature" }}{PARA 6 "" 1 "" {TEXT -1 23 "trying 1st order linear" }}{PARA 6 "" 1 "" {TEXT -1 16 "t rying Bernoulli" }}{PARA 6 "" 1 "" {TEXT -1 16 "trying separable" }} {PARA 6 "" 1 "" {TEXT -1 23 "<- separable successful" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$solG/-%\"vG6#%\"tG*&,&*&\"#q\"\"\"-%$expG6#,$*(\" \"(F.\"#D!\"\"F)F.F.F.F.F-F6F.,&F.F6F/F6F6" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 71 "Now we can check on the limit of this equation with the \+ \"limit\" routine" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 38 "limit_vel:=limit(rhs(sol),t=infinity);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%*limit_velG!#q" }}}{EXCHG {PARA 0 " " 0 "" {TEXT -1 85 "If you don't substitute at the beginning, substitu te here and you get the same answer" }}{PARA 0 "" 0 "" {TEXT -1 0 "" } }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 45 "value(subs(\{m=100,k=0.2,g =9.8\},sqrt(m*g/k)));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+++++q!\") " }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 54 "To evaluate the solution at a specific time use \"eval\"" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "eval(sol,t=0);" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#/-%\"vG6#\"\"!F'" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "mynumber:=eval(sol,t=10);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%)mynumberG/-%\"vG6#\"#5*&,&*&\"#q\"\"\"-%$expG6##\"#9 \"\"&F.F.F-!\"\"F.,&F.F5F/F5F5" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 53 "One may have to force maple to give a numerical value" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "v10:=e valf(mynumber);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$v10G/-%\"vG6#\"# 5$!+Q:Y(>'!\")" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 87 "We want to inte grate v(t) to get the distance (see the Fundamental Theorem of Calculu s)" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "sola:=int(rhs(sol),t=0..10);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%solaG,(*&\"$+&\"\"\"-%#lnG6#,&F(F(-%$expG6##\"#9\"\" &F(F(!\"\"\"$+(F(*&F'F(-F*6#\"\"#F(F(" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 39 "Add to 1200 since the value is negative" }}{PARA 0 "" 0 " " {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "dist1:=ev alf(sola);\n1200+dist1;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&dist1G$! +F#G%HQ!\"(" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+t " 0 "" {MPLTEXT 1 0 15 "eva l(sol,t=16);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"vG6#\"#;*&,&*&\"# q\"\"\"-%$expG6##\"$7\"\"#DF,F,F+!\"\"F,,&F,F3F-F3F3" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "evalf(%);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"vG6#\"#;$!+7.6Vo!\")" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 95 "L ook what happens to the velocity as t increases. At t=16 we are almos t at terminal velocity!!" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "plot(rhs(sol),t=0..16);" }}{PARA 13 "" 1 "" {GLPLOT2D 318 234 234 {PLOTDATA 2 "6%-%'CURVESG6$7S7$$\"\"! F)F(7$$\"3^mmmmYa([$!#=$!3y@.M?/3:M!#<7$$\"3cLLLL'Q?_'F-$!3G\"eFi6%*QP 'F07$$\"3ommmmWkM**F-$!3H)RA!ySlt'*F07$$\"3jmmmEr)pL\"F0$!3KsVcFq:&H\" !#;7$$\"3cch\"F@7$$\"3cmmmwG&e*>F0$!3U^]9wMd1>F@7 $$\"3v******HH1CBF0$!3G7P%GIt/?#F@7$$\"3SmmmO#)\\jEF0$!3?Fln@th&\\#F@7 $$\"3h******p\\%=+$F0$!3m%4:,2O+y#F@7$$\"3aLLLth()\\LF0$!3OA.V@(Q;1$F@ 7$$\"3rmmmYAUcOF0$!3kn[$)H]'**H$F@7$$\"3p******>0_,SF0$!3G@)[3^Zpb$F@7 $$\"3*)*******zN![VF0$!3rL\"4h&Gg-QF@7$$\"3&*******zu'>o%F0$!31F4$)R7T FSF@7$$\"3hmmmO%4_)\\F0$!3l(\\X<`i8A%F@7$$\"36LLL`MzX`F0$!3b+#[O\\+%RW F@7$$\"3\"GLLLTb7l&F0$!3(=')z04ANh%F@7$$\"3g*******G!e1gF0$!3KViKEp2/[ F@7$$\"3eKLL8I5@jF0$!3nj\\N$o[A'\\F@7$$\"3M+++!H%=mmF0$!3ZuL#HA!zC^F@7 $$\"37+++qKy%*pF0$!3$RX_G+Z#p_F@7$$\"3`LLLL=kPtF0$!3`e$zwRx'4aF@7$$\"3 ELLLBI\\_wF0$!3:TSX![#yHbF@7$$\"3_mmmmD5#*zF0$!3g>l)F0$!34Jpn\"QB%feF@7$$\"3$emm mK\"f$)*)F0$!3]z67go+`fF@7$$\"3W******f0AE$*F0$!3/w_0=9BUgF@7$$\"3M)** ***>kTh'*F0$!3P/H))GilAhF@7$$\"3v)*****\\ct&)**F0$!3cF4x!eMW>'F@7$$\"3 ')*****fo$eM5F@$!3VS:vOhlniF@7$$\"3ALLL\"QSp1\"F@$!3!yC5!pT/GjF@7$$\"3 (*******f!)[,6F@$!3/1Gi1%)G(Q'F@7$$\"3fmmm\"R$zK6F@$!3MxkBv)RmV'F@7$$ \"3(******zQ=q;\"F@$!3?dGfa8B'['F@7$$\"3;LLLU9A*>\"F@$!3)HjVUuR!HlF@7$ $\"3!******H\"H)GB\"F@$!3KNx9u-5qlF@7$$\"3FLLL`Jzl7F@$!3ENsvw=)og'F@7$ $\"3$******\\7Z-I\"F@$!3!)**o+a.6UmF@7$$\"3rmmm%RIML\"F@$!3r'zx@aBJn'F @7$$\"3immm!3ltO\"F@$!34RGiq?7-nF@7$$\"3JLLLq(=5S\"F@$!3h\"pU7h*QGnF@7 $$\"3'******f,V>V\"F@$!3)))\\R@)o^]nF@7$$\"3KLLL\"p&Qn9F@$!3@!H%y " 0 "" {MPLTEXT 1 0 28 "solb:=int(rhs(sol),t=0..16);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%solbG,(*&\"$+&\"\"\"-%#lnG6#,&F(F(-%$expG6## \"$7\"\"#DF(F(!\"\"\"%?6F(*&F'F(-F*6#\"\"#F(F(" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 14 "1200+evalf(%);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+VvQ4U!\"(" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 49 "Now set up another ODE for the post-ripcord event" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 98 "ode2:=diff(v2(t),t)=-g +(k/m)*v2(t)^2;\node2:=subs(\{m=100,k=7*0.2,g=9.8\},ode2);\nic2:=v2(0) =rhs(v10);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%ode2G/-%%diffG6$-%#v2 G6#%\"tGF,,&%\"gG!\"\"*(%\"kG\"\"\"%\"mGF/F)\"\"#F2" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#>%%ode2G/-%%diffG6$-%#v2G6#%\"tGF,,&$\"#)*!\"\"F0*&$ \"+++++9!#6\"\"\")F)\"\"#F5F5" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$ic 2G/-%#v2G6#\"\"!$!+P:Y(>'!\")" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 11 " Again solve" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 43 "sol2:=dsolve(\{ode2,ic2\},v2(t),[separable]);" }} {PARA 6 "" 1 "" {TEXT -1 33 "Classification methods on request" }} {PARA 6 "" 1 "" {TEXT -1 35 "Methods to be used are: [separable]" }} {PARA 6 "" 1 "" {TEXT -1 28 "----------------------------" }}{PARA 6 " " 1 "" {TEXT -1 38 "* Tackling ODE using method: separable" }}{PARA 6 "" 1 "" {TEXT -1 37 "--- Trying classification methods ---" }}{PARA 6 "" 1 "" {TEXT -1 16 "trying separable" }}{PARA 6 "" 1 "" {TEXT -1 23 " <- separable successful" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%sol2G/-% #v2G6#%\"tG,$*(\"#5\"\"\"\"\"(#F-\"\"#-%%tanhG6#,&**F.F-\"#]!\"\"F.F/F )F-F-*&F/F--%#lnG6#*&,&*&\"+++++5F-F.F/F-\"+P:Y(>'F-F-,&*&F?F-F.F/F-\" +P:Y(>'F7F7F-F-F-F7" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 53 "And find t he limit velocity because it is interesting" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "plot(rhs(sol2),t= 0..8);" }}{PARA 13 "" 1 "" {GLPLOT2D 236 218 218 {PLOTDATA 2 "6%-%'CUR VESG6$7Y7$$\"\"!F)$!+P:Y(>'!\")7$$\"+N3VfV!#6$!+j2x7gF,7$$\"+q;')=()F0 $!+;V:TeF,7$$\"+]#HyI\"!#5$!+Z=N\"o&F,7$$\"+MBxVbH&F,7$$\"+<$>5E$F;$!+'e&)R3&F,7$$\"+v2<9TF;$!+m$H=([F,7$$ \"+MAKn\\F;$!+&p)*Ho%F,7$$\"+N*Gh#eF;$!+&=3J^%F,7$$\"+Nc$\\o'F;$!+i/qg VF,7$$\"+=bQ%R)F;$!+#=E15%F,7$$\"+&Qk#z**F;$!+$R<7!RF,7$$\"+l9.i6!\"*$ !++f5GPF,7$$\"+>\"\\Zg#\\#Fho$!+ ![fP+$F,7$$\"+Fn*Gn#Fho$!+'Grj&HF,7$$\"+2xiDGFho$!+=0T@HF,7$$\"+Y,H.IF ho$!+mL(e)GF,7$$\"+2:bgJFho$!+pdSeGF,7$$\"+Y@4LLFho$!+\"p\"*>$GF,7$$\" +O;R(\\$Fho$!+ti)*4GF,7$$\"+<4#)oOFho$!+5c')*y#F,7$$\"+7lCEQFho$!+8$=O x#F,7$$\"+%G^g*RFho$!+UH=eFF,7$$\"+>2VsTFho$!+1a:WFF,7$$\"+O&pfK%Fho$! +f/SLFF,7$$\"+kcz\"\\%Fho$!+!3?Js#F,7$$\"+\"G5Jm%Fho$!+f+y8FF,7$$\"+6# 32$[Fho$!+mpu0FF,7$$\"+Ey'G*\\Fho$!+$f()))p#F,7$$\"+J%=H<&Fho$!+3U>#p# F,7$$\"+3>qM`Fho$!+4y!po#F,7$$\"+,.W2bFho$!+5/$>o#F,7$$\"+fp'Rm&Fho$!+ -_%zn#F,7$$\"+T>4NeFho$!+9:4uEF,7$$\"+8s5'*fFho$!+P(*)3n#F,7$$\"+mXTkh Fho$!+d0$zm#F,7$$\"+od'*GjFho$!+NdPlEF,7$$\"+EcB,lFho$!+Go,jEF,7$$\"+v >:nmFho$!+cT,hEF,7$$\"+0a#o$oFho$!+rn?fEF,7$$\"+`Q40qFho$!+riidEF,7$$ \"+\"3:(frFho$!+%\\Qjl#F,7$$\"+e%GpL(Fho$!+yO.bEF,7$$\"+:-V&\\(Fho$!+) p.Sl#F,7$$\"+ZhUkwFho$!+o7.`EF,7$$\"+F9l#F,-%'COLOURG6&%$RGBG$\"#5!\"\"F(F(-%+AXESLABELSG6$Q\"t6\"Q!Fc]l -%%VIEWG6$;F(Fd\\l%(DEFAULTG" 1 2 0 1 10 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "limit_vel2:=limit(rhs(sol2),t=infinity);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%+limit_vel2G,$*&\"#5\"\"\"\"\"(#F(\"\"#!\"\" " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "evalf(%);" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#$!+68vXE!\")" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "eval(sol2,t=0):evalf(%);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%#v2G6#\"\"!^$$!+P:Y(>'!\")$\"+ma%[V#!#<" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 59 "Now find the time it takes to go from a h eight of AGL+dist1" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 15 "First find y(t)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "sol2a:=int(rhs(sol2),t=0..s);" }}{PARA 12 "" 1 "" {XPPMATH 20 "6 #>%&sol2aG,0*&#\"$]#\"\"(\"\"\"-%#lnG6#,&-%%tanhG6#,&**F)F*\"#]!\"\"F) #F*\"\"#%\"sGF*F**&F6F*-F,6#*&,&*&\"+++++5F*F)F6F*\"+P:Y(>'F*F*,&*&F?F *F)F6F*\"+P:Y(>'F5F5F*F*F*F*F5F*F**&F'F*-F,6#,&F/F*F*F*F*F**&#\"%+XF)F *-F,6#F7F*F**&FIF*-F,6#\"\"&F*F**&F'F*-F,6#F)F*F**&#F(F)F*-F,6#,&*&F?F *F)F6F5F@F*F*F5*&#F(F)F*-F,6#F=F*F5" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 58 "Yikes! Plot it to see if it makes sense for the distance:" }} {PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "plot(sol2a,s=0..30);" }}{PARA 13 "" 1 "" {GLPLOT2D 366 227 227 {PLOTDATA 2 "6%-%'CURVESG6$7bo7$$\"+-K[V?!#6$!*[mtD\"!\")7$$\"+0k'p3%F *$!*)\\+(\\#F-7$$\"+2'\\/8'F*$!*J=&>PF-7$$\"+5G$R<)F*$!*j$[D\\F-7$$\"+ A**3E7!#5$!*qV**G(F-7$$\"+ilyM;F@$!*_'[%f*F-7$$\"+V)z@X#F@$!+I*[QS\"F- 7$$\"+DJdpKF@$!+)R&QG=F-7$$\"+(ofV!\\F@$!+m]XEEF-7$$\"+]i9RlF@$!+O_!zO $F-7$$\"+XV)RQ*F@$!+')fV`XF-7$$\"+WA)GA\"!\"*$!+!R\\%QcF-7$$\"+Qeui=F^ o$!+dqmNyF-7$$\"+i3&o]#F^o$!+)>c\"Q)*F-7$$\"+pX*y9$F^o$!+%Hg7<\"!\"(7$ $\"+WTAUPF^o$!+*)[vQ8F`p7$$\"+%*zhdVF^o$!+hoF3:F`p7$$\"+%>fS*\\F^o$!+U K*4o\"F`p7$$\"+>$f%GcF^o$!+C>_^=F`p7$$\"+Dy,\"G'F^o$!+j6)e-#F`p7$$\"+7 &)\\\"F-$! +Q/_JVF`p7$$\"+>:mk:F-$!+Pe_1XF`p7$$\"+w&QAi\"F-$!+o?%)eYF`p7$$\"+uLU% o\"F-$!+PDPB[F`p7$$\"+bjm[F-$!+JjC*\\&F`p7$$\"+:K^+?F-$!+Td*)fc F`p7$$\"+7,Hl?F-$!+j==+DF-$!+,yx%*pF`p7$$ \"+E&4Qc#F-$!+bk%\\6(F`p7$$\"+%>5pi#F-$!+!=%Q(H(F`p7$$\"+bJ*[o#F-$!+[J >UuF`p7$$\"+L+?$p#F-F]\\l7$$\"+5p],FF-F]\\l7$$\"+\\.m0FF-F]\\l7$$\"+(y 8)4FF-F]\\l7$$\"+20*=r#F-F]\\l7$$\"+Es'Rr#F-$!+8du*o(F`p7$$\"+XR/;FF-F a]l7$$\"+k17=FF-Fa]l7$$\"+=WtMFF-Fa]l7$$\"+r\"[8v#F-Fa]l7$$\"+^s1\"y#F -Fa]l7$$\"+Ijy5GF-Fa]l7$$\"+*pHm#GF-Fa]l7$$\"+nIZUGF-Fa]l7$$\"+_ZR]GF- Fa]l7$$\"+OkJeGF-Fa]l7$%*undefinedGFb_l-%'COLOURG6&%$RGBG$\"#5!\"\"$\" \"!F[`lFj_l-%+AXESLABELSG6$Q\"s6\"Q!F``l-%%VIEWG6$;Fj_l$\"#IF[`l%(DEFA ULTG" 1 2 0 1 10 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1 " }}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 101 "Maple has trouble evaluatin g with such large numbers. Let's see if it can get us the correct ans wer:" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 29 "S olve for the time of impact:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 55 "t_impact:=solve(int(rhs(sol2),t=0..s)=-(1200+dist1),s);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%)t_impactG^$$\"+E4a\\H!\")$!+TH]\"[) !\"*" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 20 "Check the real part:" }} {PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 44 "sol2b:=int(rhs(sol2),t=0..29.4954);evalf(%);" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#>%&sol2bG^$$!+A$p0<)!\"($\"+H'))G\"**!#@" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#^$$!+A$p0<)!\"($\"+H'))G\"**!#@" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 49 "Besides the (small) imaginary term, it lo oks good" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "t_impact:=Re(so l2b);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%)t_impactG$!+A$p0<)!\"(" }} }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 51 "eval(sol2,t=t_impact):\nv_i mpact:=Re(rhs(evalf(%)));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%)v_impa ctG$\"+68vXE!\")" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}} {MARK "34 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }