Нужно найти ошибку в программе

CLS
 
gamma = 1 'koeffizient zatuhaniya kolebanii
q0 = 100 'raspredelennaya nagruzka
bb = 100
hh = 100
E = 360000
c00 = 0
c10 = 0
nu = 0
L0 = 400 'dlina balki
mas = 100 'massa
g = 0.01 'vert masshtab
g1 = 10000
g2 = 0.000005
q = 500 / LL ' goriz masshtab
LL = 1 'kol-vo shagov
b1 = 4
b = L0 / b1 'shag integrirovaniya
EI = E * bb * hh * hh * hh / 12
c0 = c00 * bb
c1 = c10 * bb
h0 = 0
r11 = 24 * EI / b / b / b
mas1 = mas / (b1 + 1)
omega = SQR(r11 / mas1)
CONST pi = 3.141592653589793
h1 = 2 * pi / omega
h = h1 / 100
DIM Iv(L0) AS DOUBLE, k1(2 * n, 5) AS DOUBLE, z(2 * n) AS DOUBLE, z0(2 * n) AS DOUBLE, x0(2 * n) AS DOUBLE, x(2 * n) AS DOUBLE
DIM Jp(n, n) AS DOUBLE, Ip(n, n) AS DOUBLE, s1(b1 + 2), m1(L0), p(b1 + 2) AS DOUBLE, r(b1 + 2, b1 + 2), mom(b1 + 2) AS DOUBLE, w2(b1 + 2) AS DOUBLE, w(b1 + 2) AS DOUBLE, qu(b1 + 2) AS DOUBLE, r0(b1 + 2) AS DOUBLE, r1(b1 + 2) AS DOUBLE
q = 500 / LL
DIM r2(b1 + 2), r3(b1 + 2), f(b1 + 2), a(b1 + 2, LL), d(b1 + 2), v0(b1 + 2), v(b1 + 2), fs(b1 + 2) AS DOUBLE
SCREEN 12
LINE (10, 0)-(10, 400), 12
LINE (0, 200)-(600, 200), 12
LOCATE 12, 1: PRINT "0"
OPEN "w(x)" FOR INPUT AS #2
OPEN "r" FOR OUTPUT AS #1
OPEN "prog" FOR OUTPUT AS #3
 
FOR s = 1 TO b1 + 1
    w(s) = 0
    r0(s) = 0
    v0(s) = 0
NEXT s
FOR L = 1 TO LL
    v(1) = 0
    v(b1 + 1) = 0
    w(1) = 0
    w(b1 + 1) = 0
    FOR s = 1 TO b1 + 1
        s1(s) = b * s
        IF s = 1 THEN w(s) = 0
        IF s = b1 + 1 THEN w(s) = 0
        w2(s) = (w(s - 1) - 2 * w(s) + w(s + 1)) / b
        IF s = 3 THEN PRINT w(s - 1), w(s), w(s + 1)
        IF s = 4 THEN PRINT w(s - 1), w(s), w(s + 1)
        IF s = 1 THEN w2(s) = 0
        IF s = b1 + 1 THEN w2(s) = 0
        mom(s) = w2(s) * EI
        IF L = LL THEN LINE (10 + s1(s - 1) * 500 / L0, mom(s - 1) * g + 200)-(10 + s1(s) * 500 / L0, mom(s) * g + 200), 5
        IF s = 1 THEN qu(s) = qu(s + 1) ELSE qu(s) = (mom(s) - mom(s - 1)) / b
        p(s) = qu(s + 1) - qu(s)
        IF s = 1 THEN p(s) = qu(s)
        IF s = b1 + 1 THEN p(s) = qu(s)
        FOR m = 1 TO b1 + 1
            m1(m) = b * m
            r(s, m) = EXP(-nu * ABS(s1(s) - m1(m)))
            r1(s) = r0(s) + r(s, m)
            r0(s) = r1(s)
        NEXT m
        r2(s) = w(s) * r1(s) * c1
        r3(s) = r2(s) + w(s) * c0
        'LINE (10 + s1(s - 1) * 500 / L0, r3(s - 1) * g1 * 0.01 + 200)-(10 + s1(s) * 500 / L0, r3(s) * g1 * 0.01 + 200), 25
        mas1 = mas / (b1 + 1)
        fs(s) = -v0(s) * gamma * mas1
        IF s = 1 THEN fs(s) = -v0(s) * gamma * mas1 / 2
        IF s = b1 + 1 THEN fs(s) = -v0(s) * gamma * mas1 / 2
        f(s) = p(s) + r3(s) + fs(s)
        IF s = ((b1 / 2) + 1) THEN f(s) = f(s) - 10000 
        85 a(s, L) = f(s) / mas1
        IF s = 1 THEN a(s, L) = f(s) * 2 / mas1
        IF s = b1 + 1 THEN a(s, L) = f(s) * 2 / mas1
 
        IF L = 1 THEN d(s) = v0(s) * h + (323 * a(s, L)) * h * h / 360
        IF L = 1 THEN v(s) = v0(s) + (55 * a(s, L)) * h / 24
        IF L = 1 THEN GOTO 200
        IF L = 2 THEN d(s) = v0(s) * h + (323 * a(s, L) - 264 * a(s, L - 1)) * h * h / 360
        IF L = 2 THEN v(s) = v0(s) + (55 * a(s, L) - 59 * a(s, L - 1)) * h / 24
        IF L = 2 THEN GOTO 200
        IF L = 3 THEN d(s) = v0(s) * h + (323 * a(s, L) - 264 * a(s, L - 1) + 159 * a(s, L - 2)) * h * h / 360
        IF L = 3 THEN v(s) = v0(s) + (55 * a(s, L) - 59 * a(s, L - 1) + 37 * a(s, L - 2)) * h / 24
        IF L = 3 THEN GOTO 200
        d(s) = v0(s) * h + (323 * a(s, L) - 264 * a(s, L - 1) + 159 * a(s, L - 2) - 38 * a(s, L - 3)) * h * h / 360
        v(s) = v0(s) + (55 * a(s, L) - 59 * a(s, L - 1) + 37 * a(s, L - 2) - 9 * a(s, L - 3)) * h / 24
        200 v0(s) = v(s)
        w(s) = d(s)
        IF L = LL THEN PRINT #3, w(s)
        IF L = LL THEN LINE (10 + s1(s - 1) * 500 / L0, w(s - 1) * g1 + 200)-(10 + s1(s) * 500 / L0, w(s) * g1 + 200), 2
     NEXT s
    LINE (10 + (L - 1) * 500 / LL, d((b1 / 2) + 1) * g1 + 200)-(10 + L * 500 / LL, d((b1 / 2) + 1) * g1 + 200), 60
NEXT L
PRINT "w="; d((b1 / 2) + 1)
CLOSE #1
CLOSE #2
CLOSE #3
END