home *** CD-ROM | disk | FTP | other *** search
/ CP/M / CPM_CDROM.iso / simtel / cpmug / cpmug003.ark / FIT.BAS < prev    next >
Encoding:
BASIC Source File  |  1984-04-29  |  517 b   |  22 lines
  1. REM    LEAST SQUARE FIT OF DATA X,Y TO Y=A+BX
  2. 10    INPUT "NO. OF POINTS";NPTS
  3.     DIM X(NPTS),Y(NPTS),P(NPTS),S(NPTS)
  4.  
  5. PRINT "ENTER PAIRS OF VALUES (INDEPENDENT FIRST)"
  6.         FOR I=1 TO NPTS
  7.         INPUT X(I),Y(I)
  8.         T1=T1+X(I)
  9.         T2=T2+Y(I)
  10.         T3=T3+X(I)*Y(I)
  11.         T4=T4+X(I)^2
  12.         NEXT I
  13. REM        CALC COEFFICIENTS A,B
  14.     D=NPTS*T4-T1^2
  15.     A=(T2*T4-T1*T3)/D
  16.     B=(NPTS*T3-T1*T2)/D
  17.     PRINT
  18.     PRINT "LEAST SQUARE LINE IS: Y=";A;" + ";B;" X"
  19.     INPUT "MORE DATA (Y/N)";YN$
  20.     IF LEFT$(YN$,1) NE "Y" THEN STOP
  21. T1=0:T2=0:T3=0:T4=0:GO TO 10
  22.