5 ' HERMITE 10 CLEAR :CLS 20 PRINT "***INTERPOLACION CUBICA***" 25 LOCATE 0,1 30 PRINT "**HERMITE**" 40 INPUT "NUMERO PUNTOS";N 50 INPUT "DISTANCIA ENTRE PUNTOS";H 60 DIM A(4,N-1),S(4,4),F(2*N) 70 S(1,1)=1:S(2,2)=1:S(3,1)=1:S(4,2)=1:S(3,2)=H:S(3,3)=H*H:S(3,4)=H^3:S(4,3)=2*H:S(4,4)=3*H*H 80 CLS:PRINT "ENTRADA DATOS" 100 FOR I=0 TO N-1 105 CLS:PRINT "PUNTO";I; 110 LOCATE 0,1 115 K=2*I+1 120 INPUT "F";F(K) 125 K=K+1 130 INPUT "F'";F(K) 140 NEXT I 150 FOR I=1 TO N-1 160 FOR J=1 TO 4 170 A(J,I)=F(2*(I-1)+J) 180 NEXT J 190 NEXT I 200 'RESOLUCION SISTEMA 210 DIM P(4) 220 FOR I=1 TO 4 230 P(I)=I 240 NEXT I 250 FOR K=1 TO 3 260 M=0:KO=4 270 FOR I=K TO 4 280 W=0 290 FOR J=K TO 4 300 W=W+ABSS(P(I),J) 310 NEXT J 320 H=ABS(S(P(I),K)/W) 330 IF H>=M THEN M=H:KO=I 340 NEXT I 350 IF KO=K THEN 365 360 IF P(K)<>KO THEN M=P(KO):P(KO)=P(K):P(K)=M 365 FOR I=K+1 TO 4 370 S(P(I),K)=S(P(I),K)/S(P(K),K) 375 FOR J=K+1 TO 4 380 S(P(I),J)=S(P(I),J)-S(P(K),J)*S(P(I),K) 390 NEXT J 400 NEXT I 410 NEXT K 420 FOR R=1 TO N-1 430 DIM W(4),C(4) 440 FOR I=1 TO 4 450 C(I)=A(I,R) 460 NEXT I 470 FOR K=1 TO 4 480 S=0 490 FOR J=1 TO K-1 500 S=S+S(P(K),J)*W(J) 510 NEXT J 520 W(K)=C(P(K))-S 530 NEXT K 540 FOR K=4 TO 1 STEP-1 550 S=0 560 FOR J=4 TO K+1 STEP-1 570 S=S+S(P(K),J)*C(J) 580 NEXT J 590 C(K)=(W(K)-S)/S(P(K),K) 600 NEXT K 610 BEEP0:BEEP1 620 CLS:PRINT "COEFICIENTES POLINOMIO";R-1; 630 LOCATE 0,1 640 FOR I=1 TO 4 650 PRINT "C";I;R-1;"=";C(I) 660 NEXT I:STOP 670 NEXT R 675 CLS 680 END