home *** CD-ROM | disk | FTP | other *** search
/ C!T ROM 2 / ctrom_ii_b.zip / ctrom_ii_b / PROGRAM / PASCAL / NRPAS13 / RTSEC.DEM < prev    next >
Text File  |  1991-04-29  |  707b  |  39 lines
  1. PROGRAM d9r6(input,output);
  2. (* driver for routine RTSEC *)
  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.    root,xacc : real;
  13.    xb1,xb2 : glnbmax;
  14.  
  15. (*$I MODFILE.PAS *)
  16. (*$I BESSJ0.PAS *)
  17.  
  18. FUNCTION fx(x: real): real;
  19. BEGIN
  20.    fx := bessj0(x)
  21. END;
  22.  
  23. (*$I ZBRAK.PAS *)
  24.  
  25. (*$I RTSEC.PAS *)
  26.  
  27. BEGIN
  28.    nb := nbmax;
  29.    zbrak(x1,x2,n,xb1,xb2,nb);
  30.    writeln;
  31.    writeln('roots of bessj0:');
  32.    writeln('x':20,'f(x)':15);
  33.    FOR i := 1 to nb DO BEGIN
  34.       xacc := (1.0e-6)*(xb1[i]+xb2[i])/2.0;
  35.       root := rtsec(xb1[i],xb2[i],xacc);
  36.       writeln('  root ',i:2,'  ',root:12:6,fx(root):14:6)
  37.    END
  38. END.
  39.