bisection - test=func(xl)*func(xr); fprintf(' %1d %6.5f...

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
function bisection(func,xl,xu,es,maxit) %Problem test=func(xl)*func(xu); if test>0,error('no sign change'),end iteration=0; xr=xl; disp(' i xl xu xr f(xl) f(xr) f(xr)f(xl) ea'); while(1) xrold=xr; xr=(xl+xu)/2; iteration=iteration+1; if xr~=0, ea=abs((xr-xrold)/xr);end
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: test=func(xl)*func(xr); fprintf(' %1d %6.5f %6.5f %6.5f %6.5f %6.5f %6.5f %4.2f\n',iteration, xl,xu,xr,func(xl),func(xr), test, ea) if test<0 xu=xr; elseif test>0 xl=xr; else ea=0; end if ea<=es||iteration>=maxit,break, end end...
View Full Document

Ask a homework question - tutors are online