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 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)

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

Saved successfully!

Ooh no, something went wrong!