Lập Trình C# all Chap "NUMERICAL RECIPES IN C" part 2 docx

... Scientific Computing 15 1.3 Error, Accuracy, and Stability 28 2 Solution of Linear Algebraic Equations 32 2.0 Introduction 32 2.1 Gauss-Jordan Elimination 36 2. 2 Gaussian Elimination with Backsubstitution ... Sorting 329 8.0 Introduction 329 8.1 Straight Insertion and Shell’s Method 330 8 .2 Quicksort 3 32 8.3 Heapsort 336 8.4 Indexing and Ranking 338 8.5 Selecting the M th Largest...

Ngày tải lên: 01/07/2014, 09:20

... from NUMERICAL RECIPES IN C: THE ART OF SCIENTIFIC COMPUTING (ISBN 0- 521 -43108-5) Copyright (C) 1988-19 92 by Cambridge University Press.Programs Copyright (C) 1988-19 92 by Numerical Recipes Software. ... a single-screen license and password immediately, on-line, from the On-Line Store, with fees ranging from $50 (PC, Macintosh, educational institutions’ UNIX) to $140 (general UNI...

Ngày tải lên: 01/07/2014, 09:20

... spline 3.6 splin2 two-dimensional spline interpolation 4 .2 trapzd trapezoidal rule 4 .2 qtrap integrate using trapezoidal rule 4 .2 qsimp integrate using Simpson’s rule 4.3 qromb integrate using Romberg ... shootf 18.1 fred2 solve linear Fredholm equations of the second kind 18.1 fredin interpolate solutions obtained with fred2 18 .2 voltra linear Volterra equations of the second kind 18...

Ngày tải lên: 01/07/2014, 09:20

... http://www.nr.com or call 1-800-8 72- 7 423 (North America only),or send email to trade@cup.cam.ac.uk (outside North America). t=c/ 123 6.85; t2=t*t; as=359 .22 42+ 29.105356*c; You aren’t really intended to ... error_text[]); int i; float am,as,c,t,t2,xtra; c=n+nph/4.0; This is how we comment an individual line. 1 2 Chapter 1. Preliminaries Sample page from NUMERICAL RECIPES IN C: THE...

Ngày tải lên: 01/07/2014, 09:20

... dates to be searched. #define IYEND 20 00 int main(void) /* Program badluk */ { void flmoon(int n, int nph, long *jd, float *frac); long julday(int mm, int id, int iyyy); int ic,icon,idwk,im,iyyy,n; float ... not very interesting structurally, but it is occasionally useful: #include <math.h> #define IGREG 22 99161 void caldat(long julian, int *mm, int *id, int *iyyy) Inverse of the fu...

Ngày tải lên: 01/07/2014, 09:20

... [0] [2] [0][3] [0][4] [1][0] [1][1] [1] [2] [1][3] [1][4] [0][0] [0][1] [0] [2] [0][3] [0][4] [1][0] [1][1] [1] [2] [1][3] [1][4] [2] [0] [2] [1] [2] [2] [2] [3] [2] [4] [2] [0] [2] [1] [2] [2] [2] [3] [2] [4] *m[0] *m[1] *m [2] **m **m (a) (b) Figure ... division is interpreted as producing an integer result, throwing away any integer remainder. In floating-pointrepresentation, a number...

Ngày tải lên: 01/07/2014, 09:20

... Cliffs, NJ: Prentice Hall), Chapter 2. Johnson, L.W., and Riess, R.D. 19 82, Numerical Analysis , 2nd ed. (Reading, MA: Addison- Wesley), §1.3. Wilkinson, J.H. 1964, Rounding Errors in Algebraic Processes (Englewood ... be evaluated by summing a finite number of leading terms in its in nite series, rather than all in nity terms. In cases like this, there is an adjustable parameter...

Ngày tải lên: 01/07/2014, 09:20

... function at an intermediate point x and obtains a new, smaller bracketing interval, either (a, x) or (x, b). The process continues untilthe bracketing interval is acceptably small. It is optimal ... new bracketing triplet of points is (a, b, x); 10.0 Introduction 395 Sample page from NUMERICAL RECIPES IN C: THE ART OF SCIENTIFIC COMPUTING (ISBN 0- 521 -43108-5) Copyright (C) 1988-19 92 b...

Ngày tải lên: 01/07/2014, 09:20

... function at an intermediate point x and obtains a new, smaller bracketing interval, either (a, x) or (x, b). The process continues until the bracketing interval is acceptably small. It is optimal ... than) the square root of your machine’s floating-point precision, since smaller values will gain you nothing. It remains to decide on a strategy for choosing the new point x,given(a, b, c). Suppos...

Ngày tải lên: 01/07/2014, 09:20

... possible outcome, SHFT3(x0,x1,x2,R*x1+C*x3) its housekeeping, SHFT2(f1,f2,(*f)(x2)) and a new function evaluation. } else { The other outcome, SHFT3(x3,x2,x1,R*x2+C*x0) SHFT2(f2,f1,(*f)(x1)) and its ... your machine’s floating-point precision. #include <math.h> #include "nrutil.h" #define ITMAX 100 #define CGOLD 0.3819660 #define ZEPS 1.0e-10 Here ITMAX is the maximum allowed n...

Ngày tải lên: 01/07/2014, 09:20