5  REM NEWTON, HP 36652A, 6/74
10  REM  ***NEWTON DOES INTERPOLATIONS BY NEWTON'S FORMULA***
20  REM  ***VERSION 1  31 JULY 73  RICHARD A MILEWSKI***
30  B=0
40  DEF FNA()=(-T0)/(T1-T0)
50  DEF FNB()=B1-B0
60  DEF FNC()=(B0+B2)-2*B1
70  DEF FNZ()=B0+FNA()*FNB()+((FNA()*(FNA()-1))/2)*FNC()
80  GOTO 100
90  B=1
100  PRINT "INPUT THREE EQUIDISTANT ARGUMENTS"
110  IF B=0 THEN 130
120  PRINT "            ^^^^^^^^^^^"
130  INPUT T0,T1,T2
140  IF ABS((T1-T0)-(T2-T1))>.00001 THEN 90
150  PRINT "INPUT FUNCTION VALUES FOR THE THREE ARGUMENTS"
160  INPUT B0,B1,B2
170  PRINT "INPUT TABLE INCREMENT SIZE"
180  INPUT K
190  PRINT "INPUT TABLE LIMITS (BETWEEN";T0;" & ";T2;")"
200  INPUT K1,K2
210  IF K1<T0 OR K2>T2 THEN 190
220  PRINT LIN(5)
230  FOR T=K1 TO K2 STEP K
240  PRINT T;TAB(15);FNZ(T)
250  NEXT T
260  PRINT LIN(5)
270  END 
