Из Basic в Pascal

05 PRINT'Решение системы дифференциальных уравненийЖ
10 PRINT'Методом Рунге-Кутта-Фельберга с Автоматическим'
15 PRINT'Выбором шага':Print!2.0!
40 INPUT'Введите число уравнений N='N
45 DIM Y(N),L(N),W(N),A(N),B(N),C(N),D(N),E(N),F(N)
50 INPUT'ЗАДАЙТЕ ПОГРЕШНОСТЬ ВЫЧИСЛЕНИЙ E='E1
60 INPUT'ЗАДАЙТЕ НАЧАЛЬНЫЙ ШАГ H='H
70 INPUT'ЗАДАЙТЕ НАЧАЛЬНОЕ X0='X
80 FOR J=1 TO N:PRINT'ВВЕДИТЕ НАЧАЛЬНОЕ Y0('J')='
90 INPUT W(J):LETY(J)=W(J):NEXT J
110 LETE3=0:GOSUB 400:LETD=0:FOR J=1 TO N
120 LETA(J)=F(J)*H:LETY(J)=W(J)+2*A(J)/9:NEXT J
130 LETX=X+2*H/9:GOSUB 400:FOR J=1 TO N
140 LETB(J)=H*F(J):LETY(J)=W(J)+A(J)/12+B(J)/4:NEXT J
150 LETX=X+H/9:GOSUB 400:FOR J=1 TO N
160 LETC(J)=H*F(J):LETY(J)=W(J)+(34.5*A(J)-121.5*B(J)+135*C(J))/64
170 NEXT J:LETX=X+H/2.4:GOSUB 400:FOR J=1 TO N
180 LETD(J)=H*F(J):LETQ=W(J)-17*A(J)/12+6.75*B(J)
185 LETY(J)=Q-5.4*C(J)+16*D(J)/15:NEXT J
190 LETX=X+H/4:GOSUB 400:FOR J=1 TO N
195 LETE(J)=H*F(J):LETQ=W(J)+65*A(J)/432-.3125*B(J)
200 LETY(J)=Q+.8125*C(J)+4*D(J)/27+5*E(J)/144:NEXT J
205 LETX=X-H/6:GOSUB 400:FOR J=1 TO N
210 LETL(J)=H*F(J):LETY(J)=W(J)+A(J)/9+.45*C(J)+16*D(J)/45+E(J)/12
220 LETQ=A(J)/150-.03*C(J)+16*D(J)/75+E(J)/20
230 LETE2=ABS(Q-.24*L(J)):IF E2<=E1 THEN 250
240 LETE3=1:GOTO 260
250 IF E2<E1/20 THEN LETD=D+1
260 NEXT J:LETX=X+H/6:GOSUB 400:IF E3=0 THEN 290
270 LETX=X-H:FOR J=1 TO N:LETY(J)=W(J):NEXT J
280 LETH=H/2:GOTO 110
290 IF D=N THEN LETH=H+H
300 PRINT !F1.9!'ДЛЯ X=' X:FOR J=1 TO N
310 PRINT !2.0!'Y('J')='!F1.9!Y(J)
320 LETW(J)=Y(J):NEXT J:GOTO 110
400 LETF(1)=Y(1)+Y(2)-X*X+X-2
410 LETF(2)=-2*Y(1)+4*Y(2)+2*X*X-4*X-7
420 RETURN:END