C DERIP2 SOURCE CHAT 05/01/12 22:42:57 5004 SUBROUTINE DERIP2(P, IVAR, DP) C C Procedure de multiplication de 2 polynomes de 2 variables C Entree : P(9) : polynome a deriver : les coefficients sont C ranges comme suit : C indice : 1 2 3 4 5 6 7 8 9 C terme : 1 T2 T1 T1*T2 T2^2 T1^2 T1*T2^2 T1^2*T2 T1^2*T2^2 C IVAR : numero de la variable par rapport a laquelle il C faut deriver C Sortie : DP(9) : polynome resultat de la derivation C IMPLICIT INTEGER(I-N) IMPLICIT REAL*8(A-H,O-Z) DIMENSION P(9), DP(9), INXD1(9), INXD2(9), C1(9), C2(9) DATA INXD1/ 1, 1, 1, 2, 1, 3, 5, 4, 7/ DATA C1/0., 0., 1., 1., 0., 2., 1., 2., 2./ DATA INXD2/ 1, 1, 1, 3, 2, 1, 4, 6, 8/ DATA C2/0., 1., 0., 1., 2., 0., 2., 1., 2./ C C Initialisation du polynome resultat C DO 70 I = 1, 9 70 DP(I) = 0.D0 C C Derivation de P C IF (IVAR .EQ. 1) THEN DO 80 I = 2, 9 80 DP(INXD1(I)) = DP(INXD1(I)) + C1(I)*P(I) ELSE IF (IVAR .EQ. 2) THEN DO 90 I = 2, 9 90 DP(INXD2(I)) = DP(INXD2(I)) + C2(I)*P(I) ENDIF ENDIF RETURN END