home *** CD-ROM | disk | FTP | other *** search
/ ProfitPress Mega CDROM2 …eeware (MSDOS)(1992)(Eng) / ProfitPress-MegaCDROM2.B6I / PROG / PASCAL / NRPAS13.ZIP / FIT.DEM < prev    next >
Encoding:
Text File  |  1991-04-29  |  1.2 KB  |  54 lines
  1. PROGRAM d14r1(input,output);
  2. (* driver for routine FIT *)
  3. CONST
  4.    npt=100;
  5.    spread=0.5;
  6. TYPE
  7.    glndata = ARRAY [1..npt] OF real;
  8. VAR
  9.    gliset : integer;
  10.    glgset : real;
  11.    glinext,glinextp : integer;
  12.    glma : ARRAY [1..55] OF real;
  13.    a,b,chi2,q,siga,sigb : real;
  14.    i,idum,mwt : integer;
  15.    x,y,sig : glndata;
  16.  
  17. (*$I MODFILE.PAS *)
  18. (*$I GAMMLN.PAS *)
  19.  
  20. (*$I GSER.PAS *)
  21.  
  22. (*$I GCF.PAS *)
  23.  
  24. (*$I GAMMQ.PAS *)
  25.  
  26. (*$I RAN3.PAS *)
  27.  
  28. (*$I GASDEV.PAS *)
  29.  
  30. (*$I FIT.PAS *)
  31.  
  32. BEGIN
  33.    gliset := 0;
  34.    idum := -117;
  35.    FOR i := 1 to npt DO BEGIN
  36.       x[i] := 0.1*i;
  37.       y[i] := -2.0*x[i]+1.0+spread*gasdev(idum);
  38.       sig[i] := spread
  39.    END;
  40.    FOR mwt := 0 to 1 DO BEGIN
  41.       fit(x,y,npt,sig,mwt,a,b,siga,sigb,chi2,q);
  42.       writeln;
  43.       IF (mwt = 0) THEN BEGIN
  44.          writeln('ignoring standard deviations')
  45.       END ELSE BEGIN
  46.          writeln('including standard deviation')
  47.       END;
  48.       writeln(' ':5,'a  =  ',a:9:6,' ':6,'uncertainty:',siga:9:6);
  49.       writeln(' ':5,'b  =  ',b:9:6,' ':6,'uncertainty:',sigb:9:6);
  50.       writeln(' ':5,'chi-squared: ',chi2:14:6);
  51.       writeln(' ':5,'goodness-of-fit: ',q:10:6)
  52.    END
  53. END.
  54.