{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 }{CSTYLE "" 0 21 "" 0 1 0 0 0 1 0 0 0 0 2 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 "Normal" -1 256 1 
{CSTYLE "" -1 -1 "Times" 1 18 0 0 0 1 1 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 {EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 256 "" 0 "" 
{TEXT -1 50 "A Mathematical Model for Diffusion - Experiment #2" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart;" }}}{EXCHG {PARA 0 "
> " 0 "" {MPLTEXT 1 0 12 "with(plots):" }}}{EXCHG {PARA 0 "" 0 "" 
{MPLTEXT 0 21 31 "\nPoint plot of the data points\n" }}}{EXCHG {PARA 
0 "> " 0 "" {MPLTEXT 1 0 612 "pointplot([[0.0,\0110],[1.0,\011210.2923
54313],[2.0,\011369.121483645],[3.0,\011501.3306009],[4.0,612.2442899]
,[5.0,\011701.8625506],\n[6.0, 777.2860392],[7.0,840.2842879],[8.0,
\011894.4094414],[9.0,\011944.1016809],[10.0,\011988.4671565],[11.0,10
26.619285],\n[12.0,\0111057.671485],[13.0,1086.071203],[14.0,1109.1368
9],[15.0,1132.209844],[16.0, 1151.7292],[17.0,1170.361973],\n[18.0,118
5.448415],[19.0,1199.641007],[20.0,\0111211.181117],[21.0,1220.940795]
,[22.0,1228.927307],[23.0,1236.01997],\n[24.0,1243.119899],[25.0,1245.
786915],[26.0,1247.560081],[27.0,1251.992995],[28.0,1255.546593],[29.0
,1260.86609],[30.0,1264.419688]],symbol=circle);\n" }}}{EXCHG {PARA 0 
"> " 0 "" {MPLTEXT 1 0 529 "A1:=210.292354313;\nA2:=369.121483645;\nA3
:=501.3306009;\nA4:=612.2442899;\nA5:=701.8625506;\nA6:=777.2860392;\n
A7:=840.2842879;\nA8:=894.4094414;\nA9:=944.1016809;\nA10:=988.4671565
;\nA11:=1026.619285;\nA12:=1057.671485;\nA13:=1086.071203;\nA14:=1109.
13689;\nA15:=1132.209844;\nA16:=1151.7292;\nA17:=1170.361973;\nA18:=11
85.448415;\nA19:=1199.641007;\nA20:=1211.181117;\nA21:=1220.940795;\nA
22:=1228.927307;\nA23:=1236.01997;\nA24:=1243.119899;\nA25:=1245.78691
5;\nA26:=1247.560081;\nA27:=1251.992995;\nA28:=1255.546593;\nA29:=1260
.86609;\nA30:=1264.419688;" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}
{PARA 0 "" 0 "" {MPLTEXT 0 21 27 "Calculating the constants \n" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "A:=(n)->c/2-c/2*a^n;" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "eq1:=A(10)-A10;eq2:=A(15)-A1
5;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "sol:=solve(\{eq1,eq2
\},\{a,c\});" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "a:=.8564952
790;c:=2510.222483;" }}}{EXCHG {PARA 0 "" 0 "" {MPLTEXT 0 21 45 "\nA c
omparison between the data and the model\n" }}}{EXCHG {PARA 0 "> " 0 "
" {MPLTEXT 1 0 628 "P1:=pointplot([[0.0,\0110],[1.0,\011210.292354313]
,[2.0,\011369.121483645],[3.0,\011501.3306009],[4.0,612.2442899],[5.0,
\011701.8625506],[6.0, 777.2860392],[7.0,840.2842879],[8.0,\011894.409
4414],[9.0,\011944.1016809],[10.0,\011988.4671565],[11.0,1026.619285],
\n[12.0,\0111057.671485],[13.0,1086.071203],[14.0,1109.13689],[15.0,11
32.209844],[16.0, 1151.7292],[17.0,1170.361973],\n[18.0,1185.448415],[
19.0,1199.641007],[20.0,\0111211.181117],[21.0,1220.940795],[22.0,1228
.927307],[23.0,1236.01997],\n[24.0,1243.119899],[25.0,1245.786915],[26
.0,1247.560081],[27.0,1251.992995],[28.0,1255.546593],[29.0,1260.86609
],[30.0,1264.419688]],symbol=circle,labels=[n,An]):" }}}{EXCHG {PARA 
0 "> " 0 "" {MPLTEXT 1 0 40 "P2:=plot(c/2-c/2*a^n,n=0..30,color=red):
" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "display(P1,P2);" }}}
{EXCHG {PARA 0 "" 0 "" {MPLTEXT 0 21 42 "\nCalculating errors for a pa
rticular base\n" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "Ap:=(n)-
>c/2-c/2*a^n;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 60 "for n from
 1 to 30 do\nERR||n:=abs((A||n-Ap(n))/A||n*100);od;" }}}{EXCHG {PARA 
0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "" 0 "" {MPLTEXT 0 21 
78 "\nYour turn - try another pair of numbers (you need to replace # w
ith numbers)\n" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "a:='a';c:
='c';n:='n';" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "eq1:=A(#)-A
#;eq2:=A(#)-A#;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "sol:=sol
ve(\{eq1,eq2\},\{a,c\});" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 10 
"a:=#;c:=#;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 628 "P1:=pointpl
ot([[0.0,\0110],[1.0,\011210.292354313],[2.0,\011369.121483645],[3.0,
\011501.3306009],[4.0,612.2442899],[5.0,\011701.8625506],[6.0, 777.286
0392],[7.0,840.2842879],[8.0,\011894.4094414],[9.0,\011944.1016809],[1
0.0,\011988.4671565],[11.0,1026.619285],\n[12.0,\0111057.671485],[13.0
,1086.071203],[14.0,1109.13689],[15.0,1132.209844],[16.0, 1151.7292],[
17.0,1170.361973],\n[18.0,1185.448415],[19.0,1199.641007],[20.0,\01112
11.181117],[21.0,1220.940795],[22.0,1228.927307],[23.0,1236.01997],\n[
24.0,1243.119899],[25.0,1245.786915],[26.0,1247.560081],[27.0,1251.992
995],[28.0,1255.546593],[29.0,1260.86609],[30.0,1264.419688]],symbol=c
ircle,labels=[n,An]):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "P2
:=plot(c/2-c/2*a^n,n=0..26,color=red):" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 15 "display(P1,P2);" }}}{EXCHG {PARA 0 "" 0 "" {MPLTEXT 
0 21 42 "\nCalculating errors for a particular base\n" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "Ap:=(n)->c/2-c/2*a^n;" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 60 "for n from 1 to 30 do\nERR||n:=abs(
(A||n-Ap(n))/A||n*100);od;" }}}}{MARK "17 0 0" 17 }{VIEWOPTS 1 1 0 1 
1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }