A Brief Introduction to Fortran 77 - Polaris
A Brief Introduction to Fortran 77 - Polaris
A Brief Introduction to Fortran 77 - Polaris
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
Examplessubroutine mprove(a,alud,n,np,indx,b,x)parameter (nmax=100)dimension a(np,np), alud(np,np), indx(n),b(n), x(n), r(nmax)real*8 sdpdo i=1,nsdp=-b(i)do j=1,nsdp=sdp+dble(a(i,j))*dble(x(j))end dor(i)=sdpend docall lubksb(alud,n,np,indx,r)do i=1,nx(i)=x(i)-r(i)end doreturnend© 1998 David A. Padua24 of 118