{VERSION 2 3 "SGI IRIS UNIX" "2.3" }
{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 708 "bisect:=proc(f,xlow
_in,xhi_in,epsilon)\n   local xlow, xmid, xhi, flow, fmid, fhi;\n   \n
   xlow := xlow_in;   flow:= f(xlow);\n   xhi  := xhi_in;    fhi := f(
xhi);\n\n   if (sign(flow) = sign(fhi)) then\n       printf( `[%g,%g] \+
and [%g,%g] dont bracket a zero\\n`,\n               evalf(xlow), eval
f(flow), evalf(xhi), evalf(fhi));\n       RETURN();\n   fi;\n\n   \n  \+
 while( evalf(abs(xhi-xlow)) > evalf(epsilon)) do\n      xmid := evalf
((xlow+xhi)/2);\n      fmid := f(xmid);\n   #   print([xlow,flow], [xm
id,fmid], [xhi,fhi]);\n      if (sign(fmid) = sign(fhi)) then\n       \+
      xhi := xmid; \n             fhi := fmid;\n      else\n          \+
   xlow:= xmid;\n             flow:= fmid;\n      fi; \n   od;\n   RET
URN(xmid);\nend:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}
{MARK "1 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 }