Numerical Analysis Programs Using Fortran 90 - University of ...
Numerical Analysis Programs Using Fortran 90 - University of ...
Numerical Analysis Programs Using Fortran 90 - University of ...
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)