home *** CD-ROM | disk | FTP | other *** search
- PROCEDURE polin2(x1a: glmarray; x2a: glnarray; ya: glmbyn;
- m,n: integer; x1,x2: real; VAR y,dy: real);
- (* Programs using POLIN2 must define the types
- TYPE
- glmarray = ARRAY [1..m] OF real;
- glnarray = ARRAY [1..n] OF real;
- glmbyn = ARRAY [1..m,1..n] OF real;
- in the main routine. *)
- VAR
- k,j: integer;
- ymtmp: glmarray;
- yntmp: glnarray;
- BEGIN
- FOR j := 1 TO m DO BEGIN
- FOR k := 1 TO n DO BEGIN
- yntmp[k] := ya[j,k]
- END;
- polint(x2a,yntmp,n,x2,ymtmp[j],dy)
- END;
- polint(x1a,ymtmp,m,x1,y,dy)
- END;
-