23.10.2014 Views

Numerical Analysis Programs Using Fortran 90 - University of ...

Numerical Analysis Programs Using Fortran 90 - University of ...

Numerical Analysis Programs Using Fortran 90 - University of ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Problem: Write a program to determine the roots <strong>of</strong> the equation<br />

3 2<br />

8x − 20x<br />

− 26x<br />

+ 33 = 0 that are close to x = 0. 8 using Birge-<br />

Vieta method.<br />

dimension a(0:10),b(0:10),c(0:10)<br />

data (a(i),i=0,3)/8,-20,-26,33/<br />

m=3<br />

k=0<br />

xo=0.8<br />

tol=1e-6<br />

5 k=k+1<br />

b(0)=a(0); c(0)=b(0)<br />

do j=1,m<br />

b(j)=a(j)+xo*b(j-1)<br />

c(j)=b(j)+xo*c(j-1)<br />

enddo<br />

x=xo-b(m)/c(m-1)<br />

print *,k,x<br />

if(abs(x-xo)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!