home *** CD-ROM | disk | FTP | other *** search
/ ProfitPress Mega CDROM2 …eeware (MSDOS)(1992)(Eng) / ProfitPress-MegaCDROM2.B6I / PROG / PASCAL / NRPAS13.ZIP / PZEXTR.DEM < prev    next >
Encoding:
Text File  |  1991-04-29  |  1.0 KB  |  45 lines
  1. PROGRAM d15r8(input,output);
  2. (* driver for routine PZEXTR *)
  3. CONST
  4.    glimax=11;
  5.    glnmax=10;
  6.    glncol=7;
  7.    nv=4;
  8.    nuse=5;
  9. TYPE
  10.    glyarray = ARRAY [1..nv] OF real;
  11. VAR
  12.    glx : ARRAY [1..glimax] OF real;
  13.    glqcol : ARRAY [1..glnmax,1..glncol] OF real;
  14.    dum,xest : real;
  15.    i,iest,j : integer;
  16.    dy,yest,yz : glyarray;
  17.  
  18. (*$I MODFILE.PAS *)
  19. (*$I PZEXTR.PAS *)
  20.  
  21. BEGIN
  22. (* feed values from a rational function *)
  23. (* fn(x) := (1-x+x**3)/(x+1)**n *)
  24.    FOR i := 1 to 10 DO BEGIN
  25.       iest := i;
  26.       xest := 1.0/i;
  27.       dum := 1.0-xest+xest*xest*xest;
  28.       FOR j := 1 to nv DO BEGIN
  29.          dum := dum/(xest+1.0);
  30.          yest[j] := dum
  31.       END;
  32.       pzextr(iest,xest,yest,yz,dy,nv,nuse);
  33.       writeln;
  34.       writeln('i  :=  ',i:2);
  35.       write('Extrap. function:');
  36.       FOR j := 1 to nv DO write(yz[j]:12:6);
  37.       writeln;
  38.       write('Estimated error: ');
  39.       FOR j := 1 to nv DO write(dy[j]:12:6);
  40.       writeln
  41.    END;
  42.    writeln;
  43.    writeln('actual values:   ',1.0:12:6,1.0:12:6,1.0:12:6,1.0:12:6)
  44. END.
  45.