Skip to content

Instantly share code, notes, and snippets.

@rbakbashev
Last active June 2, 2016 20:44
Show Gist options
  • Select an option

  • Save rbakbashev/858ab51e87c49b8d0a03 to your computer and use it in GitHub Desktop.

Select an option

Save rbakbashev/858ab51e87c49b8d0a03 to your computer and use it in GitHub Desktop.
SUBROUTINE F(X, Y, J) F 00010
COMMON/A1/M1,N F 00020
COMMON/A10/NF F 00030
DOUBLE PRECISION X(N), Y(M1) F 00040
F 00050
IF (.NOT.(J.EQ.0.OR.J.EQ.-1)) GO TO 1 F 00060
Y(M1) = 3*(X(1)**2) - 2*X(2) F 00070
NF = NF+1 F 00080
1 IF (.NOT.(J.EQ.1.OR.J.EQ.-1)) GO TO 2 F 00090
Y(1) = 2*X(1) + X(2) - 4 F 00100
NF = NF+1 F 00110
2 IF (.NOT.(J.EQ.2.OR.J.EQ.-1)) GO TO 3 F 00120
Y(2) = X(1)**2 + X(2)**2 - 40 F 00130
NF = NF+1 F 00140
3 IF (.NOT.(J.EQ.3.OR.J.EQ.-1)) GO TO 4 F 00150
Y(3) = -X(1) F 00160
NF = NF+1 F 00170
4 RETURN F 00180
END F 00190
PROGRAM TAKBASHEV
REAL *8 X(2)
REAL *8 A(2),B(2)
REAL *8 P(2)
REAL *8 F
REAL *8 Y(3)
REAL *8 PAR(40)
COMMON /A10/NF /A1/M1,N,L
INTEGER N,L,M1,Q,M,NF
EXTERNAL F,CGR,CGS,UNCONS
COMMON/A5/FUNC /A8/Y1 /A9/Y2 /A13/Y11 /C81/ACTIV /C82/XPR
COMMON/C83/GR /C84/HES /C85/LZ /C86/LZZ /C87/NAPR /C88/DVOY
COMMON/C89/DVPR /C80/ZNFPR /C801/LRAB /C802/MRAB
C PAЗMEPHOCTЬ MACCИBOB GR,XPR = N
REAL*8 XPR(2),GR(2)
C PAЗMEPHOCTЬ MACCИBOB DVOY,DVPR = M
REAL*8 DVOY(2),DVPR(2)
C PAЗMEPHOCTЬ MACCИBOB Y1,Y2,Y11,FUNC,ZNFPR,ACTIV = M+1
REAL*8 FUNC(3),Y1(3),Y2(3),Y11(3),ZNFPR(3)
INTEGER ACTIV(3)
C PAЗMEPHOCTЬ MACCИBA HES = ( N,N )
REAL*8 HES(2,2)
C PAЗMEPHOCTЬ MACCИBA LZZ = ( N+M,N+M )
REAL*8 LZZ(4,4)
C PAЗMEPHOCTЬ MACCИBOB LZ,NAPR,LRAB,MRAB = ( N+M )
REAL*8 LZ(4), NAPR(4),LRAB(4),MRAB(4)
NF = 0
Q = 20
M1 = 3
M = M1-1
L=0 VC800710
N=2 VC800720
X(1) = 0
X(2) = 0
PAR(1) = 0.0001
PAR(2) = 10
PAR(3) = 0
PAR(4) = 0.8
PAR(5) = 0.01
PAR(6) = 0.1
PAR(7) = 0.1
PAR(8) = 0.0001
PAR(9) = 2
PAR(10) = 0.0001
PAR(11) = 1
PAR(12) = 1 ! Частота вывода
PAR(13) = 1 ! Детальность вывода
A(1) = -100000
A(2) = -100000
B(1) = 100000
B(2) = 100000
P(1) = 0
P(2) = 0
Y(1) = 0
Y(2) = 0
Y(3) = 0
CALL C8(N,L,M,X,A,B,P,F,CGR,CGS,Y,PAR,Q,0)
return
end
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment