home *** CD-ROM | disk | FTP | other *** search
/ ProfitPress Mega CDROM2 …eeware (MSDOS)(1992)(Eng) / ProfitPress-MegaCDROM2.B6I / PROG / PASCAL / NRPAS13.ZIP / BCUINT.DEM < prev    next >
Encoding:
Text File  |  1991-04-29  |  1.2 KB  |  50 lines
  1. PROGRAM d3r11(input,output);
  2. (* driver for routine BCUINT *)
  3. TYPE
  4.    gl4array = ARRAY [1..4] OF real;
  5.    gl4by4 = ARRAY [1..4,1..4] OF real;
  6. VAR
  7.    ansy,ansy1,ansy2,ey,ey1,ey2 : real;
  8.    x1,x1l,x1u,x1x2,x2,x2l,x2u,xxyy : real;
  9.    i : integer;
  10.    xx,y,y1,y12,y2,yy : gl4array;
  11.    wt : ARRAY [1..16,1..16] OF real;
  12.    glflag : boolean;
  13.  
  14. (*$I MODFILE.PAS *)
  15. (*$I BCUCOF.PAS *)
  16.  
  17. (*$I BCUINT.PAS *)
  18.  
  19. BEGIN
  20.    glflag := true;
  21.    xx[1] := 0.0; xx[2] := 2.0; xx[3] := 2.0; xx[4] := 0.0;
  22.    yy[1] := 0.0; yy[2] := 0.0; yy[3] := 2.0; yy[4] := 2.0;
  23.    x1l := xx[1];
  24.    x1u := xx[2];
  25.    x2l := yy[1];
  26.    x2u := yy[4];
  27.    FOR i := 1 to 4 DO BEGIN
  28.       xxyy := xx[i]*yy[i];
  29.       y[i] := sqr(xxyy);
  30.       y1[i] := 2.0*yy[i]*xxyy;
  31.       y2[i] := 2.0*xx[i]*xxyy;
  32.       y12[i] := 4.0*xxyy
  33.    END;
  34.    writeln;
  35.    writeln('x1':6,'x2':8,'y':7,'expect':11,'y1':6,
  36.       'expect':10,'y2':6,'expect':10);
  37.    writeln;
  38.    FOR i := 1 to 10 DO BEGIN
  39.       x1 := 0.2*i;
  40.       x2 := x1;
  41.       bcuint(y,y1,y2,y12,x1l,x1u,x2l,x2u,x1,x2,ansy,ansy1,ansy2);
  42.       x1x2 := x1*x2;
  43.       ey := sqr(x1x2);
  44.       ey1 := 2.0*x2*x1x2;
  45.       ey2 := 2.0*x1*x1x2;
  46.       writeln(x1:8:4,x2:8:4,ansy:8:4,ey:8:4,
  47.          ansy1:8:4,ey1:8:4,ansy2:8:4,ey2:8:4)
  48.    END
  49. END.
  50.