C DEINQT SOURCE CHAT 05/01/12 22:41:15 5004 SUBROUTINE DERINVARIPQT (SIG,P_S,Q_S,T_S,N) c la derivada t_s falta dividirla por sin(3t) c no se hace pq luego eso cancela con las derivadas de g(t), MODELO MRSLade IMPLICIT INTEGER(I-N) integer n,i real*8 sig(n),p_s(n),q_s(n),t_s(n),j2,j3,aux0,aux1,aux2 call DerInvariIJJ(sig,p_s,q_s,t_s,n) do i=1,3 p_s(i)=-p_s(i)/3.D0 enddo aux0=0.D0 aux1=0.D0 aux2=0.D0 call Invari_j2(sig,n,j2) if (j2.ge.1.D-10) then call Invari_j3(sig,n,j3) aux0=dsqrt(3.D0/j2)/2.D0 aux1=aux0/j2 aux2=1.5D0*j3/j2 endif do i=1,n t_s(i)=aux1*(aux2*q_s(i)-t_s(i)) q_s(i)=aux0*q_s(i) enddo end