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 solve the following system <strong>of</strong> equations,<br />
by using Newton-Raphson Method<br />
2<br />
2<br />
x − y − 0.2 = 0 ; y − x − 0.3 = 0<br />
Initiate the computation with guesses <strong>of</strong> x = 1. 2 and y = 1. 2 .<br />
f(x,y)=x**2-y-.2<br />
g(x,y)=y**2-x-.3<br />
fx(x,y)=2*x<br />
fy(x,y)=-1<br />
gx(x,y)=-1<br />
gy(x,y)=2*y<br />
xo=1.2<br />
yo=1.2<br />
tol=0.0001<br />
5 xj=fx(xo,yo)*gy(xo,yo)-fy(xo,yo)*gx(xo,yo)<br />
hh=(g(xo,yo)*fy(xo,yo)-f(xo,yo)*gy(xo,yo))/xj<br />
hk=(f(xo,yo)*gx(xo,yo)-g(xo,yo)*fx(xo,yo))/xj<br />
x=xo+hh<br />
y=yo+hk<br />
print*,x,y<br />
if((abs(x-xo)