{VERSION 3 0 "IBM RISC UNIX" "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 "2D Math" -1 2 "Times" 0 1 0 0 0 0 0 0 2 0 0 0 0 0 0 }{CSTYLE "2D Output" 2 20 "" 0 1 0 0 255 1 0 0 0 0 0 0 0 0 0 } {PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 0 2 2 2 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple Output" 0 11 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 3 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 11 12 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }1 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple P lot" 0 13 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 }{PSTYLE "R3 Font 0" -1 256 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 0 2 1 2 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "R3 Font 2" -1 257 1 {CSTYLE "" -1 -1 "Courier" 1 11 0 0 0 0 2 1 2 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }} {SECT 0 {EXCHG {PARA 0 "" 0 "" {TEXT -1 162 " \+ \n\n \+ A SMALL EXAMPLE FOR THE USE OF \"STOCHPLOT\" \+ " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 140 "CAL LING SEQUENCE: stochplot(x,t,p,z)\n\nstochplot plots a numerical solut ion using one of the schemes from the stochastic package for 1d SDEs\n " }}{PARA 0 "" 0 "" {TEXT -1 356 "INPUT:\n\nx: numerical scheme (as pr ovided by the output of the \"Euler\", \"Milstein\" etc. Routines from the stochastic Package\n Att ention: the \"milcomm\" ouput currently is not compatible with \"stoch plot\", use the pre-processing as in the last example below )\nt: time -step\np: number of steps to compute\nz: initial value" }}{PARA 0 "" 0 "" {TEXT -1 2 "\n\n" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "res tart;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "read stochastic5; " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "read stochplot5;" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "with(stochastic);" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 39 "Co nstruction of the Numerical Scheme..." }}{PARA 0 "" 0 "" {TEXT -1 0 " " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 50 "Euler([-(a*sin(x[1])+2*b*sin(2* x[1]))],[[sigma]]);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 48 "Eule r([-(1*sin(x[1])+2*2*sin(2*x[1]))],[[0.5]]);" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#-%&TABLEG6#7#/\"\"\"/&%#Y1G6#,&%\"nGF(F(F(,(&F+6#F.F(*& ,&-%$sinG6#F0!\"\"-F56#,$F0\"\"#!\"%F(&%&DeltaGF1F(F(*&F>F(&%#W1GF1F($ \"+++++]!#5" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 56 "\n...and output of a point-trajectory using \"stochplot\"\n\n" }{MPLTEXT 1 0 71 "stochpl ot(Euler([-(1*sin(x[1])+2*2*sin(2*x[1]))],[[0.5]]), 0.1, 10, 1);" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 34 "\nSmaller time step and more steps \n" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 73 "st ochplot(Euler([-(1*sin(x[1])+2*2*sin(2*x[1]))],[[0.5]]), 0.01, 100, 1) ;" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 38 "Same times but different initial value" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 74 "stochplot(Euler([ -(1*sin(x[1])+2*2*sin(2*x[1]))],[[0.5]]), 0.01, 100, 10);" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 34 "milcomm n eeds some pre-processing:" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 50 "milcomm([-(1*sin(x[1])+2*2*sin(2*x[ 1]))],[[0.5]]);" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#-%&TABLEG6#7#/\"\" \"/&%#Y1G6#,&%\"nGF(F(F(,(&F+6#F.F(*&,&-%$sinG6#F0!\"\"-F56#,$F0\"\"#! \"%F(&%&DeltaGF1F(F(&%(DeltaW1GF1$\"+++++]!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 100 "stochplot(subs(Delta.W.1[n]=W1[n],milcomm([-(1* sin(x[1])+2*2*sin(2*x[1]))],[[0.5]])), 0.01, 100, 1);" }}{PARA 13 "" 1 "" {INLPLOT "6%-%'POINTSG6bq7$$\"#z!\"#$\"+2%zU4\"!#57$$\"#=F)$\"+1N 9MYF,7$$\"#JF)$\"+#yrOI&!#67$$\"#fF)$!+3s,iQF77$$\"#oF)$\"+Ro9]9F,7$$ \"#\")F)$!+oDGV[F77$$\"#kF)$!+o.]nXF77$$\"\"$F)$\"+dhmq\"*F,7$$\"#;F)$ \"*_XUu*F77$$\"#WF)$!+wy)e;$F,7$$\"#bF)$\"+4*yF!HF,7$$\"#mF)$!+\"4rB'R F77$$\"#^F)$!+)\\Gap%F77$$\"##*F)$!+XNZ@7F,7$$\"\"\"F)$\"+:\"zjR)F,7$$ \"#HF)$\"+]&=)[OF,7$$\"#SF)$!+wOh/LF77$$\"#`F)$!+v7(Gy(!#77$$\"#OF)$\" +n%3m8#F,7$$\"#xF)$!+t2-dGF,7$$\"#!*F)$\"+&o003)F77$$\"#9F)$\"+pSr*R&F ,7$$\"#dF)$!+9lVpfF77$$\"#DF)$\"+G$3!)p#F,7$$\"#QF)$!*qu+f(F77$$\"#@F) $\"+,X%*p@F,7$$\"#\\F)$\"+'HTb+#F,7$$\"#vF)$\"+b7pN>F,7$$\"#))F)$!+Vz# 4\"GF77$$\"#UF)$\"+JV%GT\"F,7$$\"#5F)$\"+w4@mHF,7$$\"#BF)$!+,Jv(\\\"F, 7$$\"\"'F)$\"+YS4xxF,7$$\"#MF)$\"+qi`(4\"F,7$$\"#iF)$\"+solYLF,7$$\"#t F)$\"+_*oKd%F77$$\"#FF)$\"+O2rpBF,7$$\"#**F)$!+I^#QU\"F,7$$\"\")F)$\"+ nrjuuF,7$$\"#&*F)$!+8Z3JPF77$$\"#>F)$\"+^cJBAF,7$$\"#ZF)$\"+\"H!oj:F77 $$\"#gF)$!+AF,7$$\"#XF)$\"*`87 w&F77$$\"#')F)$!+)z#\\K;F77$$\"#pF)$\"+KA'H#=F,7$$\"##)F)$\"+9,rS7F,7$ $\"#$*F)$!+6QvQDF77$$\"#YC&F77$$\"\"#F)$\"+Cy;y&*F,7$$\"#:F)$\"+cPWUSF,7$$\"#e F)$\"+6@m.*p=HF77$$\"#aF)$\"+x?>![(F77$$\"#lF)$\"+,oj E=F,7$$\"#]F)$!+j&GoS$F,7$$\"#\"*F)$!+)H?K-$F,7$\"\"!F[p7$$\"#VF)$\"+@ xig=F,7$$\"#EF)$\"+wK$)QgF77$$\"#RF)$\"+>$yOT\"F,7$$\"#_F)$!+DY'o)QF77 $$\"#NF)$\"+$4)HY9F,7$$\"#6F)$\"+\"GGjK%F,7$$\"#CF)$\"+!=7V0\"F,7$$\"#PF)$!+i >zjDF,7$$\"#?F)$\"++By]RF,7$$\"#jF)$!+wJEjIF,7$$\"#uF)$\"+0b#>=#F]q7$$ \"#8F)$\"+WVZ)f#F,7$$\"$+\"F)$!+93SO\")F77$$\"\"*F)$\"+2k%Gf'F,7$$\"#A F)$\"+Lw'pn%F,7$$\"\"&F)$\"+4oX]mF,7$$\"#[F)$\"+9^Z ZF,$\"+'GoS8\"F]\\m" 2 601 601 601 2 0 1 0 2 9 0 3 2 1.000000 45.000000 45.000000 10030 10061 10056 10074 0 0 0 20520 0 12020 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 }}}}{MARK "0 2 0" 139 }{VIEWOPTS 1 1 0 1 1 1803 }