перевести из фортрана в матлаб

[Юлия050485]

DIMENSION V1(72), V2(72), V3(72),
%S1(141), X(71), D(141), R(101),VXX(101,
%71), VX( 101,71), VT T( 101,71), V0(72),
%V(101,71)
С ВХОДНЫЕ ДАННЫЕ
H1=0.04 : T=0.0196 : NK1 = 141
NK0=NK1 -1 : NK2=(NK1 -1 )/l 0 : M1 =2 M2=M1 : NI1=72 : NI3=NI1-2
NI4=NI1-1 : NI2=NI3/10: DO 101 I=1, N14
101 X(I)=100*ALOG(1+0.01*(I-1)*H1/2.0)
С ОПРЕДЕЛЕНИЕ ШАГА ПО X
DO 103 I=1, N13 : 103 X(I)=X(I+1)-X(I) X(NI4)=X(NI3)
С РЕШЕНИЕ ИНТ. УРАВНЕНИЯ Т.Е.
С ДАННЫЕ НА ХАРАКТЕРИСТИКАХ
GAMMA=2.0 : S1(1)=GAMMA : LN=1
DO 107 K=1,NK0: IF (K-LN) 111, 109, 113
D(K)=X((K+l)/2)/2 : LN=LN+2
GOTO 113
111 D(K)=X(K/2)/2
113 CONTINUE
107 S1(K+1)=
%S1(K)+D(K)*S 1 (K)/2.00+(D(K)*S1
%(K)/2)*(0.001*H1 *(2+0.01 *(I-1 )*H1))/
%(0.125+0.001*(I-1)*H1)
PRINT 199, ‘(S1(K)’ K=1,NK1, NI2)
199 FORMAT (5X, S(K) / (11F10.5)