home *** CD-ROM | disk | FTP | other *** search
- PROGRAM d4r7(input,output);
- (* driver for routine QGAUS *)
- CONST
- x1=0.0;
- x2=5.0;
- nval=10;
- VAR
- dx,ss,x : real;
- i : integer;
-
- FUNCTION func(x: real): real;
- BEGIN
- func := x*exp(-x)
- END;
-
- (*$I MODFILE.PAS *)
- (*$I QGAUS.PAS *)
-
- BEGIN
- dx := (x2-x1)/nval;
- writeln;
- writeln('0.0 to','qgaus':10,'expected':13);
- writeln;
- FOR i := 1 to nval DO BEGIN
- x := x1+i*dx;
- qgaus(x1,x,ss);
- writeln(x:5:2,ss:12:6,(-(1.0+x)*exp(-x)
- +(1.0+x1)*exp(-x1)):12:6)
- END
- END.
-