Есть программа, 50 строк кода на Fortrane, нужно написать аналог, но в MATLABе или в C++.
DIMENSION SG(50,10),EI(5),RAR(100),R(8),
D(5),B0(10,2)
DO 1 K=1,8
R(K)=0.2*(K-1)
DO 2 M=1,100
PAR(M)=1.2+0.048*(M-1)
DO 3 K=1,8
SMAX=0.
SMIN=100.
DO 5 N=2,11
L=N/2
C=58.68/N/N
DO 5 M=1,50
P=1.4142*PAR(M)
T=0.
Z=1.
A=T
IF(MOD(N,2).EQ.O) GO TO 8
T=3.83/N
A=ERF(T/P)
Z=EXP(-(T/P) 2)
X=Z
DO 6 I=1,L
B=7.66*I/N+T
EF=ERF(B/P)
EI(I)1.253*(EF-A)
A=EF
EX=EXP(-(B/P) 2)
D(I)=EX-X
X=EX
V=0.
DO 7 I=1,L
E1=EI(I)
DO 7 J=1,L
V=V+E1*D(J) 2/(EI(J)+R(K)*C/E1)
S1=0.25/V
SG(M,N-1)=S1
SMAX=AMAX1(SMAX,S1)
SMAIN=AMIN1(SMIN,S1)
DO 13 I=1,10
B0(I,1)=SMIN
B0(I,2)=SMAX
CALL TAGR(SG,B0,50,10,3)
CALL TAGR(SG,B0,50,10,0)
CONTINUE
STOP
END