перевести с Fortran на C++за деньги

[lyla]

Пожалуйста, помогите перевести программу с фортрана на с++ заплачу за помощь.

REAL * 8 FY(2, 40, 40)
, , KU, A, B, HKU, RA
COMMON / NXY / NX, NY
COMMON / GF / KU, A, B
PRINT 11
11 FORMAT(' КОРНИ УРАВНЕНИЯ '1, 'РЕЛЕЯ-ЛЕМБА')
NI = 3
A = 1.D0
B = 1.D0 / 3.5D0
KU = 4.71D0
HKU = 0.01D0
2 CONTINUE
CALL BBOD
DO 1 I = 1, NI
KU = KU + (I - 1)*HKU
PRINT 10, KU
10 FORMAT('ЧАСТОТН. ПАРАМ. KU=', F7.3)
CALL COORD(FY, NX, NY)
1 CONTINUE
GO TO 2
3 CONTINUE
STOP
END

-------------------------------------- -
SUBROUTINE BBOD
REAL * 8 EPS, EPS1, EPS2, H, H1, SX, SY
COMMON / EPN / EPS2, H, H1, SX, SY, / EP / EPS, EPS1, LL / NXY / NX, NY
READ 1, EPS, EPS1, EPS2
READ 2, H, H1, SX, SY
READ 3, NX, NY, LL

PRINT 4, EPS, EPS1, EPS2
PRINT 5, H, H1, SX, SY
PRINT 6, NX, NY, LL

1 FORMAT(1X, 3D7.0)
2 FORMAT(1X, 4F7.3)
3 FORMAT(1X, 3I3)
4 FORMAT(1X, ' MIN L KB -TA=', D12.5, , 'ТОЧН. ДЕЛ.ОТР НА 2=', D12.5, , 'МАКС. ПОГР. В КОРНЕ =',D12.5)
5 FORMAT ('H,H1,SX,SY==',4F7.3)
6 FORMAT ('NX,NY,LL==',3I3)
RETURN
END

-----------------------------------------

SUBROUTINE FK8(MBK,KX,KY,HX,HY)
REAL*8 MBK(2,4),KX,KY,HX,HY
MBK(1,1)=KX
MBK(2,1)=KY
MBK(2,2)=KY+HY
MBK(1,3)=KX+HX
MBK(2,3)=KY+HY
MBK(1,4)=KX+HX
MBK(2,4)=KY
RETURN
END

_________________________________

SUBROUTINE COORD(FY,NX,NY)
REAL*8 SF(2,4), BK(2,4),KX,KY,,H,H1,SX,SY,FI,PSI,,EPS,EPS1,EPS2,,FY(2,NY,NX)
COMMON/FF/FI/PSI,/EPN/EPS2,H,H1,SX,SY,/EP/EPS,EPS1,LL,/CK/IND
IND=0
NX1=NX-1
NY1=NY-1
DO 28 I=1,NX
DO 28 J=1,NY
KX=(I-1)*H+SX
KY=(J-1)*H1+SY
CALL FUN(KX,KY)
FY(1,J,I)=DSIGN(1.D0,FI)
28 FY(2,J,I)=DSIGN(1.D0,PSI)
DO 1 I=1, NX1
DO 2 J=1, NY1
KX=(I-1)*H*SX
KY=(J-1)*H1*SY
CALL FKB(BK,KX,KY,H,H1)
DO 3 K=1,2
SF(K,1)=FY(K,J,I)
SF(K,2)=FY(K,J+1,I)
SF(K,3)=FY(K,J+1,I+1)
3 SF(K,4)=FY(K,J,I+1)
CALL KBADM(IS,BK,SF)
IF(IS.EQ.0) GO TO 2
CALL KNA4M(SF,KX,KY,H,H1)
2 CONTINUE
1 CONTINUE
KX=SX+H*(NX-1)
KY=SY+H1*(NY-1)
IF(IND.LT.1) PRINT 11,SX,SY,KX,KY
11 FORMAT('КОРНЕЙ НЕТ ИСХОДНЫЙ ПРЯМОУГЛ.'1,'SX,SY,KX,KY==',4F10,6)
RETURN
END

------------------------------------

    SUBROUTINE P2D(FP,A1,B1,A2,B2,IZ,G,S)
    INTEGER FP
    REAL*8 FP,A1,B1,A2,B2,IZ,G,S,FI,PSI,
    ,W,T,R,EPS,EPS1,WP,EP,TP,RP,
    ,R1,R2,AS,BS,A,B,SR,EL
    COMMON/FF/FI,PSI
    COMMON/EP/EPS,EPS1,LL
    ,/KP/WP,EP,TP,RP
    R1=A1
    P1=B1
    R2=A2
    P2=B2
    G=A1
    S=B1
    AS=A1
    BS=B1
    L=0
    W=WP
    IF(FP.EQ.1) W=EP
    R=TP
    IF(FP.EQ.1) R=RP
5   L=L+1
    IF(IZ.EQ.0) AS=(A1+A2)/2.D0
    IF(IZ.EQ.1) BS=(B1+B2)/2.D0
    CALL FUN(AS,BS)
    SP=FI
    IF(FP.EQ.1) SP=PSI
    T=DSIGN(1.D0,SR)
    IF((T*R),GT,0,D0) GO TO 2
    W=T
    IF(IZ.EQ.0) A1=AS
    IF(IZ.EQ.1) B1=BS
    GO TO  3 
2   IF((W*T).GT,0.D0)  GO TO 10
    R=T
    IF(IZ.EQ.0) A2=AS
    IF(IZ.EQ.1) B2=BS
3   IF (L,GT,LL) GO TO 4
    IF(IZ.EQ.0) EL= A2-A1
    IF(IZ.EQ.1) EL= B2-B1
    EL=DABS(EL)
    IF(EL.GT.EPS1) GO TO 5
    GO TO 6 
4   PRINT 7,L
7   FORMAT(" L,LL B P2D  L= " , 15) 
    GO TO 11 
10  PRINT 12,A1,B1,A2,B2,FP,R1,P1,R2,P2,IZ 
12  FORMAT ("ТОЧКА ПЕРЕВАЛА ИСЧЕЗЛА НА ОТРЕЗКЕ "/,
    ,/"ИСХОДНЫЙ  ОТРЕЗОК  И ПАРАМЕТРЫ P2D:"/,
    ,"FR=",12, "A1,B1=", 2E15.7,"A2,B2=",2E15.7,"IZ=",12) 
20  FORMAT("ИСЧ ПРИ L=",15)
    GO TO 11
6   IF(IZ.EQ.0) G=(A1+A2)/2.D0
    IF(IZ.EQ.1) S=(B1+B2)/2.D0
11  CONTINUE
    RETURN 
    END

[_vek21]

Готов выполнить за денежное вознаграждение.

Моё объявление и отзывы находятся тут:
http://forum.sources.ru/index.php?showtopic=313018

[WorldMaster]

еще актуально?