home *** CD-ROM | disk | FTP | other *** search
/ ProfitPress Mega CDROM2 …eeware (MSDOS)(1992)(Eng) / ProfitPress-MegaCDROM2.B6I / PROG / PASCAL / NRPAS13.ZIP / BESSI0.DEM < prev    next >
Encoding:
Text File  |  1991-04-29  |  552 b   |  25 lines
  1. PROGRAM d6r21(input,output,dfile);
  2. (* driver for routine BESSI0 *)
  3. VAR
  4.    i,nval : integer;
  5.    val,x : real;
  6.    txt : string[27];
  7.    dfile : text;
  8.  
  9. (*$I MODFILE.PAS *)
  10. (*$I BESSI0.PAS *)
  11.  
  12. BEGIN
  13.    glopen(dfile,'fncval.dat');
  14.    REPEAT readln(dfile,txt) UNTIL (txt = 'Modified Bessel Function I0');
  15.    readln(dfile,nval);
  16.    writeln(txt);
  17.    writeln;
  18.    writeln('x':5,'actual':16,'bessi0(x)':17);
  19.    FOR i := 1 to nval DO BEGIN
  20.       readln(dfile,x,val);
  21.       writeln(x:6:2,val:16:7,bessi0(x):16:7)
  22.    END;
  23.    close(dfile)
  24. END.
  25.