1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
|
!----------------------------------------------------------------------
!----------------------------------------------------------------------
! abc : The A --> B --> C reaction
!----------------------------------------------------------------------
!----------------------------------------------------------------------
SUBROUTINE FUNC(NDIM,U,ICP,PAR,IJAC,F,DFDU,DFDP)
! ---------- ----
IMPLICIT NONE
INTEGER, INTENT(IN) :: NDIM, ICP(*), IJAC
DOUBLE PRECISION, INTENT(IN) :: U(NDIM), PAR(*)
DOUBLE PRECISION, INTENT(OUT) :: F(NDIM)
DOUBLE PRECISION, INTENT(INOUT) :: DFDU(NDIM,NDIM), DFDP(NDIM,*)
DOUBLE PRECISION X1,X2,X3,D,ALPHA,BETA,B,S,E,X1C
X1=U(1)
X2=U(2)
X3=U(3)
D=PAR(1)
ALPHA=PAR(2)
BETA=PAR(3)
B=PAR(4)
S=PAR(5)
E=DEXP(X3)
X1C=1-X1
F(1)=-X1 + D*X1C*E
F(2)=-X2 + D*E*(X1C - S*X2)
F(3)=-X3 - BETA*X3 + D*B*E*(X1C + ALPHA*S*X2)
END SUBROUTINE FUNC
!----------------------------------------------------------------------
SUBROUTINE STPNT(NDIM,U,PAR,T)
! ---------- -----
IMPLICIT NONE
INTEGER, INTENT(IN) :: NDIM
DOUBLE PRECISION, INTENT(INOUT) :: U(NDIM),PAR(*)
DOUBLE PRECISION, INTENT(IN) :: T
PAR(1)=0.0
PAR(2)=1.0
PAR(3)=1.55
PAR(4)=8.
PAR(5)=0.04
U(1)=0.
U(2)=0.
U(3)=0.
END SUBROUTINE STPNT
!----------------------------------------------------------------------
SUBROUTINE BCND
END SUBROUTINE BCND
SUBROUTINE ICND
END SUBROUTINE ICND
SUBROUTINE FOPT
END SUBROUTINE FOPT
SUBROUTINE PVLS
END SUBROUTINE PVLS
|
