{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 }{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Heading 1" 0 3 1 {CSTYLE "" -1 -1 "" 1 18 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }1 0 0 0 6 6 0 0 0 0 0 0 -1 0 }} {SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 161 "restart;\nwith(plot s):\nalpha:=2;beta:=0.01;c:=1;delta:=0.01;\nsys:=diff(b(t),t)=-alpha*b (t)+beta*b(t)*l(t),diff(l(t),t)=c*l(t)-delta*b(t)*l(t);\n\nfcns:=\{b(t ),l(t)\};" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 219 "sol1:=dsolve( \{sys,b(0)=15,l(0)=150\},fcns,numeric,method=classical[rk4]);\nsol2:=d solve(\{sys,b(0)=150,l(0)=150\},fcns,numeric,method=classical[rk4]);\n sol3:=dsolve(\{sys,b(0)=150,l(0)=50\},fcns,numeric,method=classical[rk 4]);" }}}{PARA 3 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 273 "p1:=odeplot(sol1,[b(t),l(t)],0..6,labels=[Nagyhal,Ki shal],numpoints=500,color=blue):\np2:=odeplot(sol2,[b(t),l(t)],0..6,la bels=[Nagyhal,Kishal],numpoints=500,color=red):\np3:=odeplot(sol3,[b(t ),l(t)],0..6,labels=[Nagyhal,Kishal],numpoints=500,color=green):\ndisp lay(p1,p2,p3);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 162 "p1:=odep lot(sol2,[t,b(t)],0..20,labels=[t,hal],numpoints=100,color=blue):\np2: =odeplot(sol2,[t,l(t)],0..20,labels=[t,hal],numpoints=100,color=red): \ndisplay(p1,p2);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}} {MARK "9" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }