{VERSION 5 0 "APPLE_PPC_MAC" "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 }{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 }} {SECT 0 {EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 50 " higher-order undetermined coefficients (#4.3.12)" }}{PARA 0 " " 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 32 "clear workspace & load ODE stuff" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart:" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 14 "with(DEtools):" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 14 "define forcing" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "G := 12*sin(t);" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 20 "define solution form" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "Y := A*cos(t) + B*sin(t);" } }}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 19 " substitute into ODE" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 95 "ode \+ := collect(diff(Y,t$4) + 2*diff(Y,t$3) + diff(Y,t$2) + 8*diff(Y,t) - 1 2*Y,\{sin(t),cos(t)\});" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }} {PARA 0 "" 0 "" {TEXT -1 35 "solve for undetermined coefficients" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 "co := solve(identity(ode=G,t ),\{A,B\});" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "part_sol := \+ subs(co,Y);" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 " " {TEXT -1 5 "voila" }}}}{MARK "13 0 0" 0 }{VIEWOPTS 1 1 0 3 2 1804 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }