{VERSION 5 0 "IBM INTEL NT" "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 "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 "Heading 1" 0 3 1 
{CSTYLE "" -1 -1 "" 1 18 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }1 0 0 0 8 4 0 0 
0 0 0 0 -1 0 }{PSTYLE "Maple Output" -1 11 1 {CSTYLE "" -1 -1 "Times" 
1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }3 3 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 "Maple Plot
" -1 13 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }
3 1 0 0 0 0 1 0 1 0 2 2 0 1 }}
{SECT 0 {SECT 1 {PARA 3 "" 0 "" {TEXT -1 19 "Initialize packages" }}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 27 "with(plots):\nwith(DEtools):
" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#%5Math~696.200~and~700G" }}{PARA 
12 "" 1 "" {TEXT -1 0 "" }}}}{SECT 1 {PARA 3 "" 0 "" {TEXT -1 17 "Defi
ne parameters" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 75 "M := t -> 8
0-15*cos(Pi*t/12);\nH := t -> 0;\nU := t -> 0;\nT0 := 75;\nK := 0.1;" 
}}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"MGf*6#%\"tG6\"6$%)operatorG%&arr
owGF(,&\"#!)\"\"\"*&\"#:F.-%$cosG6#,$*&#F.\"#7F.*&%#PiGF.9$F.F.F.F.!\"
\"F(F(F(" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"HG\"\"!" }}{PARA 11 "
" 1 "" {XPPMATH 20 "6#>%\"UG\"\"!" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>
%#T0G\"#v" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"KG$\"\"\"!\"\"" }}}
{PARA 0 "" 0 "" {TEXT -1 28 "Plot outside air temperature" }}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "plot(M(t),t=0..24);" }}{PARA 13 "" 
1 "" {GLPLOT2D 484 484 484 {PLOTDATA 2 "6%-%'CURVESG6$7en7$$\"\"!F)$\"
#lF)7$$\"3*)*******\\ech#!#=$\"3x2uiIb^.l!#;7$$\"3x********pJJ_F/$\"3J
t@!RkXS^'F27$$\"30+++vu=2vF/$\"3GGE(*pr()GlF27$$\"3M+++]z0$y*F/$\"3#e#
Q;&oH*[lF27$$\"3++++qm>!\\\"!#<$\"3_%yT\"G?r7mF27$$\"3=+++!p![0?FE$\"3
GF*G@OR?q'F27$$\"3))*****\\l:$=DFE$\"3yU'*)\\>iV\"oF27$$\"3#)*****\\Jz
P*HFE$\"3uFdfuQhPpF27$$\"3')*****\\R%4'[$FE$\"3u$o5zjJD3(F27$$\"31+++b
tC&*RFE$\"3y*fwjE%Q[sF27$$\"3')*****\\XnF]%FE$\"3b.g3`!zpU(F27$$\"3J++
+gU\"[-&FE$\"3xyuJe?>@wF27$$\"31+++qLj%[&FE$\"3Ce!HjZH#)z(F27$$\"3V+++
!y!G-gFE$\"3cL$e,m&*3+)F27$$\"3%)*******p`?_'FE$\"32>\\b`BP/#)F27$$\"3
Q+++?7&H-(FE$\"3M-j?At#pR)F27$$\"3#******\\:9yZ(FE$\"3?mXIRj'fc)F27$$
\"36+++!=!p=!)FE$\"3'\\I)o!HZjv)F27$$\"35+++?J)oZ)FE$\"3Vwp5@M#f!*)F27
$$\"3%)*****\\Vq)4!*FE$\"3)f!))o^uRj!*F27$$\"3K******>Xl\"[*FE$\"3%*p]
`I1j&=*F27$$\"3_+++NkF****FE$\"3I(4:2,'*))H*F27$$\"31++]!\\<#\\5F2$\"3
WnNwS/k%Q*F27$$\"32+++vik+6F2$\"3mRaEENa\\%*F27$$\"3*****\\U6gU7\"F2$
\"3M)G(G@$31Z*F27$$\"33++]`R(y9\"F2$\"3mtC$G]ag[*F27$$\"3%****\\#pYMt6
F2$\"3#o!G,!=\\j\\*F27$$\"3)******\\Q:))>\"F2$\"3Y#e!R'y#***\\*F27$$\"
3)****\\-Is_A\"F2$\"3I<*=w1=n\\*F27$$\"3)*****\\:#H<D\"F2$\"3O8^*Gqli[
*F27$$\"3)******z.gZF\"F2$\"3s$3q%p9Or%*F27$$\"3'*****\\g3z(H\"F2$\"3%
\\ac-#)46X*F27$$\"31+++*pQvM\"F2$\"3)RF:AL*[*Q*F27$$\"3#******R3L*)R\"
F2$\"3'=mr,#z7,$*F27$$\"3%******HY7#\\9F2$\"3!)RJ9m,\">>*F27$$\"3+++]Z
.'y\\\"F2$\"3!et.6*[em!*F27$$\"3)*******Gb(=b\"F2$\"3*)[8:lzG2*)F27$$
\"3$)*****>d5/g\"F2$\"3'fD\\ME.'[()F27$$\"3/+++!4KAl\"F2$\"3/zE/Bs\"fc
)F27$$\"3!*****\\(3!>*p\"F2$\"3mx\\.Q)*H\"R)F27$$\"39+++#eF0v\"F2$\"3q
hi<<]t$>)F27$$\"3%)****\\j@$))z\"F2$\"3]:k1wge/!)F27$$\"3/++]pVK\\=F2$
\"3,Gb*G^To!yF27$$\"3!*******H(*o)*=F2$\"3]lc3*zVnh(F27$$\"3!*****\\(o
q.&>F2$\"3x!*)*3V-jCuF27$$\"3:+++#fX,+#F2$\"3#GNtW\"\\]\\sF27$$\"3-+++
@w/^?F2$\"3-o$44W(f$3(F27$$\"39++]b\"G:5#F2$\"3Nn()yU\\5NpF27$$\"37+++
CX\"z9#F2$\"3\\[phTK(\\\"oF27$$\"32+++P&y5?#F2$\"33G+8t$\\))p'F27$$\"3
++++k!H'[AF2$\"3&4oiu&)\\ih'F27$$\"31++]VyK*H#F2$\"3k?P#*pmz^lF27$$\"3
))*****RW!fBBF2$\"3ccr+!)>\"*HlF27$$\"31++]WI&yM#F2$\"392!p(pm&R^'F27$
$\"3A++DAl#RP#F2$\"3*R&e+NK\\.lF27$$\"#CF)F*-%'COLOURG6&%$RGBG$\"#5!\"
\"F(F(-%+AXESLABELSG6$Q\"t6\"Q!Fj]l-%%VIEWG6$;F(F]]l%(DEFAULTG" 1 2 0 
1 10 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" }}}}}
{SECT 1 {PARA 3 "" 0 "" {TEXT -1 28 "Set up differential equation" }}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 46 "ode := diff(T(t),t) = K*(M(t
)-T(t))+H(t)+U(t);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$odeG/-%%diffG
6$-%\"TG6#%\"tGF,,($\"#!)!\"\"\"\"\"*&$\"#:F0F1-%$cosG6#,$*(\"#7F0%#Pi
GF1F,F1F1F1F0*&$F1F0F1F)F1F0" }}}{PARA 0 "" 0 "" {TEXT -1 29 "setup in
itial conditions also" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "ic \+
:= T(0)=T0;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#icG/-%\"TG6#\"\"!\"#
v" }}}}{SECT 1 {PARA 3 "" 0 "" {TEXT -1 27 "Solve initial value proble
m" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "dsolve( \{ode, ic\}, T(
t) );" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"TG6#%\"tG,&**\"\"&\"\"\"
-%$expG6#,$*&\"#5!\"\"F'F+F2F+,&\"#sF2*&\"#DF+)%#PiG\"\"#F+F+F+,&\"#OF
+*&F6F+F7F+F+F2F2*(F1F+,*\"$)GF2*&\"$+#F+F7F+F2*&\"#aF+-%$cosG6#,$*(\"
#7F2F8F+F'F+F+F+F+*(\"#XF+F8F+-%$sinGFFF+F+F+F:F2F2" }}}{PARA 0 "" 0 "
" {TEXT -1 24 "find the right hand side" }}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 7 "rhs(%);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,&**\"\"&\"
\"\"-%$expG6#,$*&\"#5!\"\"%\"tGF&F-F&,&\"#sF-*&\"#DF&)%#PiG\"\"#F&F&F&
,&\"#OF&*&F2F&F3F&F&F-F-*(F,F&,*\"$)GF-*&\"$+#F&F3F&F-*&\"#aF&-%$cosG6
#,$*(\"#7F-F4F&F.F&F&F&F&*(\"#XF&F4F&-%$sinGFBF&F&F&F6F-F-" }}}{PARA 
0 "" 0 "" {TEXT -1 35 "create a function for the solution " }}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "sol := unapply(%,t);" }}{PARA 11 "
" 1 "" {XPPMATH 20 "6#>%$solGf*6#%\"tG6\"6$%)operatorG%&arrowGF(,&**\"
\"&\"\"\"-%$expG6#,$*&#F/\"#5F/9$F/!\"\"F/,&\"#sF8*&\"#DF/)%#PiG\"\"#F
/F/F/,&\"#OF/*&F<F/F=F/F/F8F8*(F6F/,*\"$)GF8*&\"$+#F/F=F/F8*&\"#aF/-%$
cosG6#,$*&#F/\"#7F/*&F>F/F7F/F/F/F/F/*(\"#XF/F>F/-%$sinGFLF/F/F/F@F8F8
F(F(F(" }}}}{EXCHG }{SECT 1 {PARA 3 "" 0 "" {TEXT -1 44 "Plot solution
 versus outside air temperature" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 
1 0 29 "plot(\{M(t),sol(t)\}, t=0..24);" }}{PARA 13 "" 1 "" {GLPLOT2D 
315 303 303 {PLOTDATA 2 "6&-%'CURVESG6$7en7$$\"\"!F)$\"#lF)7$$\"3*)***
****\\ech#!#=$\"3x2uiIb^.l!#;7$$\"3x********pJJ_F/$\"3Jt@!RkXS^'F27$$
\"30+++vu=2vF/$\"3GGE(*pr()GlF27$$\"3M+++]z0$y*F/$\"3#e#Q;&oH*[lF27$$
\"3++++qm>!\\\"!#<$\"3_%yT\"G?r7mF27$$\"3=+++!p![0?FE$\"3GF*G@OR?q'F27
$$\"3))*****\\l:$=DFE$\"3yU'*)\\>iV\"oF27$$\"3#)*****\\JzP*HFE$\"3uFdf
uQhPpF27$$\"3')*****\\R%4'[$FE$\"3u$o5zjJD3(F27$$\"31+++btC&*RFE$\"3y*
fwjE%Q[sF27$$\"3')*****\\XnF]%FE$\"3b.g3`!zpU(F27$$\"3J+++gU\"[-&FE$\"
3xyuJe?>@wF27$$\"31+++qLj%[&FE$\"3Ce!HjZH#)z(F27$$\"3V+++!y!G-gFE$\"3c
L$e,m&*3+)F27$$\"3%)*******p`?_'FE$\"32>\\b`BP/#)F27$$\"3Q+++?7&H-(FE$
\"3M-j?At#pR)F27$$\"3#******\\:9yZ(FE$\"3?mXIRj'fc)F27$$\"36+++!=!p=!)
FE$\"3'\\I)o!HZjv)F27$$\"35+++?J)oZ)FE$\"3Vwp5@M#f!*)F27$$\"3%)*****\\
Vq)4!*FE$\"3)f!))o^uRj!*F27$$\"3K******>Xl\"[*FE$\"3%*p]`I1j&=*F27$$\"
3_+++NkF****FE$\"3I(4:2,'*))H*F27$$\"31++]!\\<#\\5F2$\"3WnNwS/k%Q*F27$
$\"32+++vik+6F2$\"3mRaEENa\\%*F27$$\"3*****\\U6gU7\"F2$\"3M)G(G@$31Z*F
27$$\"33++]`R(y9\"F2$\"3mtC$G]ag[*F27$$\"3%****\\#pYMt6F2$\"3#o!G,!=\\
j\\*F27$$\"3)******\\Q:))>\"F2$\"3Y#e!R'y#***\\*F27$$\"3)****\\-Is_A\"
F2$\"3I<*=w1=n\\*F27$$\"3)*****\\:#H<D\"F2$\"3O8^*Gqli[*F27$$\"3)*****
*z.gZF\"F2$\"3s$3q%p9Or%*F27$$\"3'*****\\g3z(H\"F2$\"3%\\ac-#)46X*F27$
$\"31+++*pQvM\"F2$\"3)RF:AL*[*Q*F27$$\"3#******R3L*)R\"F2$\"3'=mr,#z7,
$*F27$$\"3%******HY7#\\9F2$\"3!)RJ9m,\">>*F27$$\"3+++]Z.'y\\\"F2$\"3!e
t.6*[em!*F27$$\"3)*******Gb(=b\"F2$\"3*)[8:lzG2*)F27$$\"3$)*****>d5/g
\"F2$\"3'fD\\ME.'[()F27$$\"3/+++!4KAl\"F2$\"3/zE/Bs\"fc)F27$$\"3!*****
\\(3!>*p\"F2$\"3mx\\.Q)*H\"R)F27$$\"39+++#eF0v\"F2$\"3qhi<<]t$>)F27$$
\"3%)****\\j@$))z\"F2$\"3]:k1wge/!)F27$$\"3/++]pVK\\=F2$\"3,Gb*G^To!yF
27$$\"3!*******H(*o)*=F2$\"3]lc3*zVnh(F27$$\"3!*****\\(oq.&>F2$\"3x!*)
*3V-jCuF27$$\"3:+++#fX,+#F2$\"3#GNtW\"\\]\\sF27$$\"3-+++@w/^?F2$\"3-o$
44W(f$3(F27$$\"39++]b\"G:5#F2$\"3Nn()yU\\5NpF27$$\"37+++CX\"z9#F2$\"3
\\[phTK(\\\"oF27$$\"32+++P&y5?#F2$\"33G+8t$\\))p'F27$$\"3++++k!H'[AF2$
\"3&4oiu&)\\ih'F27$$\"31++]VyK*H#F2$\"3k?P#*pmz^lF27$$\"3))*****RW!fBB
F2$\"3ccr+!)>\"*HlF27$$\"31++]WI&yM#F2$\"392!p(pm&R^'F27$$\"3A++DAl#RP
#F2$\"3*R&e+NK\\.lF27$$\"#CF)F*-%'COLOURG6&%$RGBG$\"#5!\"\"F(F(-F$6$7Y
7$F($\"#vF)7$F4$\"3'3;Ck^t#\\uF27$F>$\"3uB!*QEKO3uF27$FC$\"3/(G^BWvpO(
F27$FI$\"3pH`d4QDJtF27$FN$\"3b+51!3PDI(F27$FS$\"3A`k)4+'o#G(F27$FX$\"3
%Qp-\\[_&psF27$$\"3=+++v3nSPFE$\"3QqkC?l'eE(F27$Fgn$\"3C@gyS'eVE(F27$$
\"3S+++0u+\\UFE$\"3))pg^5p0lsF27$F\\o$\"3)=%H'R;wzE(F27$$\"33++]d3zjZF
E$\"3+]%*pR&4LF(F27$Fao$\"3/`K%R^35G(F27$Ffo$\"3!Q/)>O%*G+tF27$F[p$\"3
hsQMvxbItF27$F`p$\"3oC`K'z:(ptF27$Fep$\"33*G#4AyD:uF27$Fjp$\"3+P:Q,Asi
uF27$F_q$\"3]]K**=V#f_(F27$Fdq$\"3qLFi\\!\\We(F27$Fiq$\"3!e'yxsRCdwF27
$F^r$\"3_rB=?![\\s(F27$Fcr$\"3Ft`^U-h,yF27$Fhr$\"3[Q>n$\\!zvyF27$F]s$
\"3#fF!)[RvJ&zF27$Fgs$\"36'fh.qPJ-)F27$Fat$\"3m[J@4OD'4)F27$F[u$\"3]Pk
BCjOo\")F27$Feu$\"3A%p'RMp(pA)F27$Fju$\"39\"pl'))[(\\G)F27$F_v$\"36!39
?Hw\"Q$)F27$Fdv$\"3[k'eS,%y#Q)F27$Fiv$\"3HV-_.sF=%)F27$F^w$\"3UATENvA[
%)F27$Fcw$\"3QOmG?&RiY)F27$$\"3w*****4L@ji\"F2$\"3EWUtV\\Js%)F27$Fhw$
\"3!G#p')Qv*eZ)F27$$\"3!)***\\()36dn\"F2$\"3sr'p\\(Q)pZ)F27$F]x$\"3!)=
r\")Q#=gZ)F27$$\"3-++vM)e[s\"F2$\"3%*4&o*\\fis%)F27$Fbx$\"3WS!G.E:oY)F
27$Fgx$\"3I\"4[$)G[%\\%)F27$F\\y$\"3%*f'p!zxiA%)F27$Fay$\"3AO#[)p@L)Q)
F27$Ffy$\"3!**=S\"e=bW$)F27$F[z$\"3w&>.\"36b&H)F27$F`z$\"3')>vn#fC%R#)
F27$Fez$\"3Kk4T;#p(y\")F27$Fjz$\"3E'yszF\"f>\")F27$F_[l$\"3Wb(=/22*[!)
F27$Fd[l$\"31?-F'y*>%)zF27$Fi[l$\"3LF[h$fl[\"zF27$Fc\\l$\"3k+bxCXI\\yF
27$F]]l$\"31PyLh*y4y(F2-F`]l6&Fb]lF(Fc]lF(-%+AXESLABELSG6$Q\"t6\"Q!F`i
l-%%VIEWG6$;F(F]]l%(DEFAULTG" 1 2 0 1 10 0 2 9 1 4 2 1.000000 
45.000000 45.000000 0 0 "Curve 1" "Curve 2" }}}}}{EXCHG {PARA 0 "> " 
0 "" {MPLTEXT 1 0 0 "" }}}}{MARK "5" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 
1 1 }{PAGENUMBERS 0 1 2 33 1 1 }