{VERSION 3 0 "IBM INTEL NT" "3.0" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 }{CSTYLE "" -1 256 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 } {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 -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 0 0 0 8 4 0 0 0 0 0 0 -1 0 }{PSTYLE "Heading 2" 3 4 1 {CSTYLE "" -1 -1 "" 1 14 0 0 0 0 0 0 0 0 0 0 0 0 0 }0 0 0 -1 8 2 0 0 0 0 0 0 -1 0 }{PSTYLE "Title" 0 18 1 {CSTYLE "" -1 -1 "" 1 18 0 0 0 0 0 1 1 0 0 0 0 0 0 }3 0 0 -1 12 12 0 0 0 0 0 0 19 0 }{PSTYLE "" 0 256 1 {CSTYLE "" -1 -1 "" 0 1 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 256 "" 0 "" {TEXT -1 0 "" }{TEXT 256 96 "Departme nts of Chemical Engineering and Mathematics\nMaple Session 5\nNumerica l Computing in Maple" }}}{EXCHG {PARA 18 "" 0 "" {TEXT -1 28 "Code for the Cobweb Diagrams" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 4 "" 0 "" {TEXT -1 16 "1. Still version" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "g:=unapply(x+f(x)/2,x) ;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "thisx:=1.0;" }{TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 160 "pointsequence := [thisx,thisx]:\nfor i from 1 to 15 do\n nextx:=g(thisx);\n pointseq uence := \n pointsequence,[thisx,nextx], [nextx,nextx];\n thisx:=n extx;\nod:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 47 "plot1:=pointp lot([pointsequence],connect=true):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 46 "plot2:=plot([g(x),x],x=-2.2..2.2,y=-2.2..2.2):" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 41 "display(plot2,plot1,scaling= constrained);" }{TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 4 "" 0 "" {TEXT -1 19 "2. Anim ated version" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "g:=unapply(x+f(x)/2,x);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "thisx:=1.0;" }{TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "pointsequence := [thisx,thisx];" }{TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 46 "plot1:=plot([g(x),x],x =-2.2..2.2,y=-2.2..2.2):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "plotsequence := plot1:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 390 "for i from 1 to 15 do\n nextx:=g(thisx);\n pointsequence := pointse quence,[thisx,nextx];\n plot2:=pointplot([pointsequence],connect=true );\n plot3:=display(plot2,plot1);\n plotsequence:=plotsequence,plot3 ;\n pointsequence := pointsequence,[nextx,nextx];\n plot2:=pointplot ([pointsequence],connect=true);\n plot3:=display(plot2,plot1);\n plo tsequence:=plotsequence,plot3;\n thisx:=nextx;\nod:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 58 "display(plotsequence,insequence=true,scal ing=constrained);" }{TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}}}{MARK "0 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 }