Numerical Recipes in Fortran 77  The Art of Scientific Computing 

ISBN : 052143064X 

Cover Design  Numerical Recipes in Fortran 77  The Art of Scientific Computing 
For your free electronic copy of this book please verify the numbers below. (We need to do this to make sure you're a person and not a malicious script)  
Sample Chapter From Numerical Recipes in Fortran 77  The Art of Scientific Computing Copyright © William H. Press, Brian P. Flannery, Saul A. Teukolsky, William T. Vetterling 

Chapter 1. Preliminaries1.0 IntroductionThis book, like its predecessor edition, is supposed to teach you methods ofnumerical computing that are practical, efficient, and (insofar as possible) elegant. We presume throughout this book that you, the reader, have particular tasks that you want to get done. We view our job as educating you on how to proceed. Occasionally we may try to reroute you briefly onto a particularly beautiful side road; but by and large, we will guide you along main highways that lead to practical destinations. Throughout this book, you will find us fearlessly editorializing, telling you what you should and shouldn’t do. This prescriptive tone results from a conscious decision on our part, and we hope that you will not find it irritating. We do not claim that our advice is infallible! Rather, we are reacting against a tendency, in the textbook literature of computation, to discuss every possible method that has ever been invented, without ever offering a practical judgment on relative merit. We do, therefore, offer you our practical judgments whenever we can. As you gain experience, you will form your own opinion of how reliable our advice is. We presume that you are able to read computer programs in FORTRAN, that being the language of this version of Numerical Recipes (Second Edition). The book Numerical Recipes in C (Second Edition) is separately available, if you prefer to program in that language. Earlier editions of Numerical Recipes in Pascal and Numerical Recipes Routines and Examples in BASIC are also available; while not containing the additional material of the Second Edition versions in C and FORTRAN, these versions are perfectly serviceable if Pascal or BASIC is your language of choice. When we include programs in the text, they look like this: SUBROUTINE
flmoon(n,nph,jd,frac)
INTEGER jd,n,nph REAL frac,RAD PARAMETER (RAD=3.14159265/180.) Our programs begin with an introductory comment summarizing their purpose and explaining their calling sequence. This routine calculates the phases of the moon. Given an integer n and a code nph for the phase desired (nph = 0 for new moon, 1 for first quarter, 2 for full, 3 for last quarter), the routine returns the Julian Day Number jd, and the fractional part of a day frac to be added to it, of the nth such phase since January, 1900. Greenwich Mean Time is assumed. INTEGER i
REAL am,as,c,t,t2,xtra c=n+nph/4. t=c/1236.85 t2=t**2 as=359.2242+29.105356*c am=306.0253+385.816918*c+0.010730*t2 ! jd=2415020+28*n+7*nph xtra=0.75933+1.53058868*c+(1.178e41.55e7*t)*t2 if(nph.eq.0.or.nph.eq.2)then xtra=xtra+(0.17343.93e4*t)*sin(RAD*as)0.4068*sin(RAD*am) else if(nph.eq.1.or.nph.eq.3)then xtra=xtra+(0.17214.e4*t)*sin(RAD*as)0.6280*sin(RAD*am) else pause ’nph is unknown in flmoon’ endif if(xtra.ge.0.)then i=int(xtra) else i=int(xtra1.) endif jd=jd+i frac=xtrai return END
