1 'INTERPOL 5 CLEAR :CLS 10 PRINT "Interpolacion por NEWTON"; 20 LOCATE 0,1 25 INPUT "Numero puntos";N 30 GR=N-1 35 PRINT "Interpolacion por un polinomio de grado";GR 40 DIM X(N),D(N) 50 CLS:PRINT"ENTRADA DATOS"; 60 LOCATE 0,1 70 FOR I=0 TO GR 80 INPUT "Punto";X(I) 90 PRINT "F(";X(I);")";:INPUT D(I) 100 NEXT I 110 FOR K=1 TO GR 120 FOR I=0 TO GR-K 130 D(I)=(D(I+1)-D(I))/(X(I+K)-X(I)) 140 NEXT I 145 NEXT K 150 B=D(0) 160 INPUT "Punto para interpolar";Z 170 FOR I=1 TO GR 180 B=D(I)+(Z-X(I))*B 190 NEXT I 200 PRINT "F(";Z;")=";B 210 INPUT "Otro punto S/N";A$ 220 IF A$="S" THEN 150 ELSE IF A$<>"N" THEN 210 230 END 240 CLEAR:N=20 250 DIM D(N),C(N),A(N-1),X(N),B(N) 260 FOR P=2 TO N 270 PRINT "C(";P;")=";:INPUT C(P) 277 NEXT P 278 FOR P=1 TO N 280 PRINT "D(";P;")=";:INPUT D(P) 283 PRINT "B(";P;")=";:INPUT B(P) 285 NEXT P 290 FOR P=1 TO N-1 295 PRINT "A(";P;")=";:INPUT D(P) 297 NEXT P 300 FOR K=2 TO N 303 IF D(K-1)<>0 THEN GOTO 305 ELSE 322 305 M=C(K)/D(K-1) 310 D(K)=D(K)-M*A(K-1) 315 B(K)=B(K)-M*B(K-1) 320 NEXT K 322 IF D(N)<>0 THEN 325 ELSE END 325 X(N)=B(N)/D(N) 330 FOR J=N-1 TO 1 STEP -1 340 X(J)=(B(J)-A(J)*X(J+1))/D(J) 345 NEXT J 350 FOR G=1 TO N 355 PRINT "X(";G;")=";X(G) 360 NEXT G