REM Programa de resoluci¢n de la ecuaci¢n y''(x)+12 x*x=0 mediante
REM m‚todos iterativos en diferencias: Gauss-Seidel (omega=1) y
REM superrelajaci¢n (0<omega<2)
REM                 VEASE LA SECCIàN 6.2 DE KOONIN
CLS
n = 50  'n debe ser par
h = 1 / n
omega = 1.9
DIM y(0 TO n)
DIM f(0 TO n)

FOR i = 0 TO n
        x = i * h
        f(i) = -h * h * 12 * x * x
        y(i) = 0
NEXT i
PRINT "omega="; omega, "n puntos="; n
PRINT "y en x=1/2 exacto", 7 / 16

FOR it = 1 TO 701
        FOR i = 1 TO n - 1
                yp = (y(i - 1) + y(i + 1) - f(i)) / 2
                y(i) = (1 - omega) * y(i) + omega * yp
        NEXT i
        IF (it - 1) MOD 50 <> 0 THEN GOTO 100
        PRINT USING "iteraci¢n=######  y(x=1/2)=#.######"; it; y(n / 2)
        

100  NEXT it