input "Dimensiones de la rejilla? (N,M)"; N, M
input "Cuáles son las temperaturas de las fronteras de la placa (Temperatura Izquierda, derecha, superior, inferior)?"; TempIzq, TempDer, TempSup, TempInf
input "Cuál es la primera estimación de la temperatura en el interior?"; TempInt

'Inicializar variables
dim T(N,M)
let TempAnt=TempInt
let Convergencia$="No"
let itn=0

'Los indices de las matrices en este paquete de programación están invertidos

'Establece las temperaturas izquierda y derecha de la placa
for j=0 to N
    T(j,0)=TempIzq
    T(j,M)=TempDer
next j

'Establece las temperaturas de arriba y abajo.
for j=1 to M
    T(0,j)=TempSup
    T(N,j)=TempInf
next j

'Establece las temperaturas del interior de la placa
for j=1 to (N-1)
    for k=1 to (M-1)
        T(j,k)=TempInt
    next k
next j
let SumEps = 0.0

' (Re)Calcular la temperatura del interior de la placa
for it=1 to 100 'máximo de 100 iteraciones para alcanzar convergencia
    SumEps=0.0
    for j=1 to N-1
        for k=1 to M-1
            TempAnt=T(j,k)
            T(j,k)=(T(j-1,k) + T(j+1,k) + T(j,k-1) + T(j,k+1))/4.0
            SumEps=SumEps + abs(T(j,k)-TempAnt)
        next k
    next j
    if SumEps<=0.08 then
        Convergencia$="Si" 'El error permitido es de 0.10
        it=100 'terminar con el bucle
    end if
    itn=it
next it

'Salida de resultados
print Convergencia$ + " hay convergencia en ";itn;" iteraciones"
for j=0 to N
    for k=0 to M
        print T(j,k),
    next k
    print
next j
print SumEps

end