{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 }{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 }}
{SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart;" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 85 "gauss:=proc()\n  RETURN(sqrt(-2*ln(
rand()/1e12)) * cos(2*evalf(Pi)*rand()/1e12)):\nend:" }}}{EXCHG {PARA 
0 "> " 0 "" {MPLTEXT 1 0 485 "wiener:=proc(N)\n  local W, t, i, erg, D
eltam, Delta, Mult, k:\n  \n  for i from 0 to 2^N do W[i]:=0.0 od:\n\n
  W[2^N]:=gauss():\n\n  for i from 1 to N do\n    Deltam := 2^i;\n    \+
Delta := 1/Deltam:\n    Mult := 2^(N-i):\n    for k from 1 by 2 to Del
tam do\n      W[k*Mult]:= evalf((W[(k-1)*Mult] + W[(k+1)*Mult])/2 + sq
rt(Delta/2)*gauss()):\n      #print(k*Mult, W[K*Mult]):\n    od:\n  od
:\n\n\n  erg := [0.0, 0.0]:\n\n  for i from 1 to 2^N do\n    erg := er
g, [i/2^N, W[i]]:\n  od:\n\n  RETURN(erg):\n\nend:\n" }}}{EXCHG {PARA 
0 "> " 0 "" {MPLTEXT 1 0 79 "_seed:=2: \nerg:=wiener(9):\nwith(plots):
\nPLOTNUM:=plot([erg]):\ndisplay(PLOTNUM);" }}}{EXCHG {PARA 0 "> " 0 "
" {MPLTEXT 1 0 0 "" }}}}{MARK "5" 0 }{VIEWOPTS 1 1 0 1 1 1803 }