DEFDBL A-Z
DIM s(50), i(50)
a = .00001#
b = 1 / 14
s(0) = 45400
i(0) = 2100
FOR k = 1 TO 50
        sum = 0
        FOR j = 0 TO k - 1
                sum = sum + s(j) * i(k - j - 1)
        NEXT j
        s(k) = -a * sum / k
        i(k) = -s(k) - b * i(k - 1) / k
NEXT k