home *** CD-ROM | disk | FTP | other *** search
/ ProfitPress Mega CDROM2 …eeware (MSDOS)(1992)(Eng) / ProfitPress-MegaCDROM2.B6I / PROG / PASCAL / NRPAS13.ZIP / ZBRAK.DEM < prev    next >
Encoding:
Text File  |  1991-04-29  |  684 b   |  36 lines
  1. PROGRAM d9r3(input,output);
  2. (* driver for routine ZBRAK *)
  3. CONST
  4.    n=100;
  5.    nbmax=20;
  6.    x1=1.0;
  7.    x2=50.0;
  8. TYPE
  9.    glnbmax = ARRAY [1..nbmax] OF real;
  10. VAR
  11.    i,nb : integer;
  12.    xb1,xb2 : glnbmax;
  13.  
  14. (*$I MODFILE.PAS *)
  15. (*$I BESSJ0.PAS *)
  16.  
  17. FUNCTION fx(x: real): real;
  18. BEGIN
  19.    fx := bessj0(x)
  20. END;
  21.  
  22. (*$I ZBRAK.PAS *)
  23.  
  24. BEGIN
  25.    nb := nbmax;
  26.    zbrak(x1,x2,n,xb1,xb2,nb);
  27.    writeln;
  28.    writeln('brackets for roots of bessj0:');
  29.    writeln('lower':22,'upper':10,
  30.          'f(lower)':16,'f(upper)':10);
  31.    FOR i := 1 to nb DO BEGIN
  32.       writeln('  root ',i:2,' ':4,xb1[i]:10:4,xb2[i]:10:4,
  33.          ' ':4,fx(xb1[i]):10:4,fx(xb2[i]):10:4)
  34.    END
  35. END.
  36.