Programme de calcul de la fonction elliptique température—durée de développement.
PROG.
M.
SAOUD
-
DIMENSION A (20, 20), B (20), c (5, 6), D (6), X (20), Y % 20 0, YK % 20 0, Q % 18 0
999
READ 200, (Q (1), I = 1, 18)
READ 100, N
DO 1 I
=
1, N
1
READ 101, X (I), Y (1)
,
DO 2 I
=
1, N
A(I,1) = X(I)**2
A(I,2)=2.*X(I)*Y(I)
‘
‘
A(I,3)=Y(I)**2
‘
.
‘
-
A(I,4) =2.*X(I)
A(I,5)=2.*Y(I)
_
A (I, 6) = — 1.
,
2
B (I) = o.
PRINT 2000, (Q (1), I = 1, 18)
PRINT 1000
'
—
PRINT 1001, (X (1), Y (1), I = 1, N)
DO 3 I
=
1, 5
D0 3 ) = 1, 5
s
=
o.
,
_
DD 4 K
=
1, N
4
s
=
s
+A(K,I)*A(K,I)
3
(: (I, ]) = s
'
DD 7 I
=
1, 5
s
=
o.
DO 8 K
=
1, N
'
8
s
=
s
+
A (K, 1)
A (K,—6)
'
7
C (I, 6) = s
.
,
CALL SM2 (C, D, 5, I]K)
_
PRINT 1035, (D (I), I = 1, 5)
_
. , =
-
_
1035 FORMAT (/10 X, 5 (2 X, E 14. 8)/)
'
‘
'
'