конвертация кода

PROGRAM oscillators ! вычисляется смещение N связанных осциляторов
DIM u(0 to 21),vel(20)
CALL initial(N,u,vel,kc,k,dt,tmax)
CALL screen(N,nplot,tmax,dy)
CALL move(N,u,vel,kc,k,dt,nplot,tmax,dy)
END
 
SUB initial(N,u(),vel(),kc,k,dt,tmax)
    INPUT prompt "число частиц = ":N
    INPUT prompt "шаг по времени = ":dt
    INPUT prompt "продолжительность = ":tmax
    INPUT prompt "силовая постоянная kc = ":kc ! внутреннее пружинки
    LET k = 1                  ! пристеночные пружинки
    DATA 0.5,0,0                   ! данные для N = 2
    FOR i = 1 to N                 ! начальные условия
        LET vel(i) = 0
        READ u(i)
    NEXT i
END SUB
 
SUB screen(N,nplot,tmax,dy)
    LET dy = 2      ! расстояние между графиками на экране
    SET WINDOW -1,tmax + 1,-dy,3*dy
    LET ntick = 100
    LET dx = tmax/ntick ! расстояние между делениями
    LET Ly = 0.1*dx     ! "высота" делений
    LET nplot = min(4,N)    ! число осцилляторов для рисования графиков
    LET row = -1
    FOR iplot = 1 to nplot
        PLOT LINES: 0,row; tmax,row ! отмечаем равновесное расстояние
        FOR itick = 1 to ntick      ! рисуем деление оси
            LET col = itick*dx
            PLOT LINES: col,row; col,row + Ly   ! рисуем деления
        NEXT itick
        LET row = row + dy
    NEXT iplot
END SUB
 
SUB move(N,u(),vel(),kc,k,dt,nplot,tmax,dy)
    DIM a(20)
    DO
    LET t = t +dt       ! время
    FOR i = 2 to (N - 1)    ! ускорение точек, не связанных со стенками
        LET a(i) = kc*(u(i+1) + u(i-1) - 2*u(i))
    NEXT i
    LET a(1) = kc*(u(2) - u(1)) - k*u(1)    ! ускорение концевых масс
    LET a(N) = kc*(u(N-1) - u(N)) - k*u(N)
    FOR i = 1 to N      ! алгоритм Эйлера-Кромера
        LET vel(i) = vel(i) + a(i)*dt
        LET u(i) = u(i) + vel(i)*dt
    NEXT i
    LET row = -1
    FOR iplot = 1 to nplot
        PLOT POINTS: t,u(iplot) + row
        LET row = row + dy
    NEXT iplot
    LOOP until t > tmax
END SUB