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

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

View Full Document Right Arrow Icon
function falseposition(func,xl,xu,es,maxit) %Problem 5.8 falseposition test=func(xl)*func(xu); if test>0,error('no sign change'),end iteration=0; xr=xl; disp(' i xl xu f(xl) f(xu) xr f(xr) f(xl)f(xr) ea'); while(1) xrold=xr; xr=xu-(func(xu)*(xl-xu))/(func(xl)-func(xu)); 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 %2.1f %6.5f %6.5f %6.5f %6.5f %6.5f %6.5f %5.3f\n',iteration, xl,xu,func(xl),func(xu),xr,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

This note was uploaded on 09/27/2011 for the course EGM 3344 taught by Professor Raphaelhaftka during the Spring '09 term at University of Florida.

Ask a homework question - tutors are online