10 'CHOLESKY 15 CLS 20 CLEAR 30 PRINT "Metodo de Cholesky"; 35 LOCATE 0,1 40 INPUT "Orden del sistema";N 45 DIM A(N,N),B(N),X(N) 50 FOR I=1 TO N 60 FOR J=I TO N 70 PRINT "A(";I;J;")"; 80 INPUT A(I,J):A(J,I)=A(I,J) 82 NEXT J 84 NEXT I 90 CLS: PRINT "Terminos independientes"; 95 LOCATE 0,1 100 FOR I=1 TO N 110 PRINT "B(";I;")"; 120 INPUT B(I) 130 NEXT I 135 INPUT "Ver datos intermedios";V$ 140 'Descomposicion 150 FOR I=1 TO N 160 FOR J=I TO N 170 SUM=A(I,J) 180 FOR K=1 TO I-1 190 SUM=SUM-A(K,I)*A(K,J) 200 NEXT K 210 IF I<>J THEN A(I,J)=SUM*TEMP ELSE TEMP=1/SQRSUM:A(I,J)=TEMP 215 IF V$="S" THEN PRINT "A'(";I;J;")=";:IF I<>J THEN PRINT A(I,J) ELSE PRINT 1/A(I,J) 220 NEXT J 230 NEXT I 240 'Obtencion de B' 250 FOR I=1 TO N 260 SUM=B(I) 270 FOR K=1 TO I-1 280 SUM=SUM-A(K,I)*X(K) 290 NEXT K 300 X(I)=SUM*A(I,I) 310 IF V$="S" THEN PRINT "B'(";I;")=";X(I) 315 NEXT I 320 'Obtencion de X 330 FOR I=N TO 1 STEP -1 340 SUM=X(I) 350 FOR K=I+1 TO N 360 SUM=SUM-A(I,K)*X(K) 370 NEXT K 380 X(I)=SUM*A(I,I) 390 NEXT I 400 BEEP 0:CLS:PRINT "Resultados"; 405 LOCATE 0,1 410 FOR I=1 TO N 420 PRINT "X(";I;")=";X(I) 430 NEXT I 440 END