double dAdAlnpsi(nu1,nu2,yx1,yx2,yx3,yy1,yy2,yy3,yz1,yz2,yz3,x,y,z,t)
double yx1,yx2,yx3,yy1,yy2,yy3,yz1,yz2,yz3,nu1,nu2,x,y,z,t;
{
 double Pi=3.14159265358979324;
 double sinh(),cosh(),cos(),sqrt();
 double psi,psih,Lpsi,LLpsi,dpsi2,Ldpsi2p2dpsidLpsi;
 double L1psi,L2psi,L3psi,dpsiddpsi2;
 double d100psi,d010psi,d001psi;
 double d100psi2,d010psi2,d001psi2;
 double d110psi,d200psi,d020psi,d011psi,d101psi,d002psi;
 double d110psi2,d200psi2,d020psi2,d011psi2,d101psi2,d002psi2;
 double d120psi,d210psi,d300psi,d030psi,d102psi,d201psi,d111psi;
 double d012psi,d021psi,d003psi;
 double d130psi,d310psi,d220psi,d400psi,d040psi,d103psi,d301psi,d202psi;
 double d013psi,d031psi,d022psi,d004psi,d112psi,d211psi,d121psi;
 double dApsi,dipsi2,d00dAlnpsi,diidAlnpsi;
 double c0,s02,c1,c2,c3,s1,s2,s3,x1,x2,x3;
 double r1,r2,r3,r1r2h,r1r3h,r2r3h,nu3;
 double r12,r22,r32,r1r2h2,r2r3h2,r1r3h2,m1,m2,m3;
 double u1,u11,u12,u13,u14,u2,u21,u22,u23,u24;
 double v1,v11,v12,v13,v14,v2,v21,v22,v23,v24;
 double u3,u31,u32,u33,u34,v3,v31,v32,v33,v34;
 double c11,c12,c13,c14,c21,c22,c23,c24,c31,c32,c33,c34;
     
 nu3=1-nu1-nu2;
 r1=sqrt((x-yx1)*(x-yx1)+(y-yy1)*(y-yy1)+(z-yz1)*(z-yz1));
 r2=sqrt((x-yx2)*(x-yx2)+(y-yy2)*(y-yy2)+(z-yz2)*(z-yz2));
 r3=sqrt((x-yx3)*(x-yx3)+(y-yy3)*(y-yy3)+(z-yz3)*(z-yz3));
 r12=r1*r1;
 r22=r2*r2;
 r32=r3*r3;
 r1r2h=((x-yx1)*(x-yx2)+(y-yy1)*(y-yy2)+(z-yz1)*(z-yz2))/r1/r2;
 r1r3h=((x-yx1)*(x-yx3)+(y-yy1)*(y-yy3)+(z-yz1)*(z-yz3))/r1/r3;
 r2r3h=((x-yx2)*(x-yx3)+(y-yy2)*(y-yy3)+(z-yz2)*(z-yz3))/r2/r3;
 r1r2h2=r1r2h*r1r2h;
 r1r3h2=r1r3h*r1r3h;
 r2r3h2=r2r3h*r2r3h;
 x1=sqrt((yx1-yx2)*(yx1-yx2)+(yy1-yy2)*(yy1-yy2)+(yz1-yz2)*(yz1-yz2));
 x2=sqrt((yx3-yx2)*(yx3-yx2)+(yy3-yy2)*(yy3-yy2)+(yz3-yz2)*(yz3-yz2));
 x3=sqrt((yx1-yx3)*(yx1-yx3)+(yy1-yy3)*(yy1-yy3)+(yz1-yz3)*(yz1-yz3));
 m1=2*Pi*nu1;
 m2=2*Pi*nu2;
 m3=2*Pi*nu3;
 c0=cos(2*Pi*t);
 s02=1.0-c0*c0;
 c1=cosh(m1*r1);
 c2=cosh(m2*r2);
 c3=cosh(m3*r3);
 s1=sinh(m1*r1);
 s2=sinh(m2*r2);
 s3=sinh(m3*r3);
     
 u1=s1/r1;
 u11=(m1*c1-u1)/r1;
 u12=(m1*m1*s1-2*u11)/r1;
 u13=(m1*m1*m1*c1-3*u12)/r1;
 u14=(m1*m1*m1*m1*s1-4*u13)/r1;
 u2=s2/r2;
 u21=(m2*c2-u2)/r2;
 u22=(m2*m2*s2-2*u21)/r2;
 u23=(m2*m2*m2*c2-3*u22)/r2;
 u24=(m2*m2*m2*m2*s2-4*u23)/r2;
 u3=s3/r3;
 u31=(m3*c3-u3)/r3;
 u32=(m3*m3*s3-2*u31)/r3;
 u33=(m3*m3*m3*c3-3*u32)/r3;
 u34=(m3*m3*m3*m3*s3-4*u33)/r3;
 v1=s1*r1;
 v11=m1*c1*r1+s1;
 v12=m1*(m1*s1*r1+2*c1);
 v13=m1*m1*(m1*c1*r1+3*s1);
 v14=m1*m1*m1*(m1*s1*r1+4*c1);
 v2=s2*r2;
 v21=m2*c2*r2+s2;
 v22=m2*(m2*s2*r2+2*c2);
 v23=m2*m2*(m2*c2*r2+3*s2);
 v24=m2*m2*m2*(m2*s2*r2+4*c2);
 v3=s3*r3;
 v31=m3*c3*r3+s3;
 v32=m3*(m3*s3*r3+2*c3);
 v33=m3*m3*(m3*c3*r3+3*s3);
 v34=m3*m3*m3*(m3*s3*r3+4*c3);
 c11=m1*s1;
 c12=m1*m1*c1;
 c13=m1*m1*m1*s1;
 c14=m1*m1*m1*m1*c1;
 c21=m2*s2;
 c22=m2*m2*c2;
 c23=m2*m2*m2*s2;
 c24=m2*m2*m2*m2*c2;
 c31=m3*s3;
 c32=m3*m3*c3;
 c33=m3*m3*m3*s3;
 c34=m3*m3*m3*m3*c3;
 
 psih=c1*c2*c3+c3*u2*v1/2+c2*u3*v1/2+c3*u1*v2/2+c1*u3*v2/2+
  c2*u1*v3/2+c1*u2*v3/2+c2*c3*u1*x1/2+c1*c3*u2*x1/2+
  c1*c2*u3*x1/2+u1*u2*v3*x1/2+c2*c3*u1*x2/2+c1*c3*u2*x2/2+
  c1*c2*u3*x2/2+u2*u3*v1*x2/2+c3*u1*u2*x1*x2/2+
  c1*u2*u3*x1*x2/2+c2*c3*u1*x3/2+c1*c3*u2*x3/2+c1*c2*u3*x3/2+
  u1*u3*v2*x3/2+c3*u1*u2*x1*x3/2+c2*u1*u3*x1*x3/2+
  c2*u1*u3*x2*x3/2+c1*u2*u3*x2*x3/2+u1*u2*u3*x1*x2*x3/2;
 d100psi=c11*c2*c3+c3*u2*v11/2+c2*u3*v11/2+c3*u11*v2/2+c11*u3*v2/2+
  c2*u11*v3/2+c11*u2*v3/2+c2*c3*u11*x1/2+c11*c3*u2*x1/2+
  c11*c2*u3*x1/2+u11*u2*v3*x1/2+c2*c3*u11*x2/2+c11*c3*u2*x2/2+
  c11*c2*u3*x2/2+u2*u3*v11*x2/2+c3*u11*u2*x1*x2/2+
  c11*u2*u3*x1*x2/2+c2*c3*u11*x3/2+c11*c3*u2*x3/2+
  c11*c2*u3*x3/2+u11*u3*v2*x3/2+c3*u11*u2*x1*x3/2+
  c2*u11*u3*x1*x3/2+c2*u11*u3*x2*x3/2+c11*u2*u3*x2*x3/2+
  u11*u2*u3*x1*x2*x3/2;
 d010psi=c1*c21*c3+c3*u21*v1/2+c21*u3*v1/2+c3*u1*v21/2+c1*u3*v21/2+
  c21*u1*v3/2+c1*u21*v3/2+c21*c3*u1*x1/2+c1*c3*u21*x1/2+
  c1*c21*u3*x1/2+u1*u21*v3*x1/2+c21*c3*u1*x2/2+c1*c3*u21*x2/2+
  c1*c21*u3*x2/2+u21*u3*v1*x2/2+c3*u1*u21*x1*x2/2+
  c1*u21*u3*x1*x2/2+c21*c3*u1*x3/2+c1*c3*u21*x3/2+
  c1*c21*u3*x3/2+u1*u3*v21*x3/2+c3*u1*u21*x1*x3/2+
  c21*u1*u3*x1*x3/2+c21*u1*u3*x2*x3/2+c1*u21*u3*x2*x3/2+
  u1*u21*u3*x1*x2*x3/2;
 d001psi=c1*c2*c31+c31*u2*v1/2+c2*u31*v1/2+c31*u1*v2/2+c1*u31*v2/2+
  c2*u1*v31/2+c1*u2*v31/2+c2*c31*u1*x1/2+c1*c31*u2*x1/2+
  c1*c2*u31*x1/2+u1*u2*v31*x1/2+c2*c31*u1*x2/2+c1*c31*u2*x2/2+
  c1*c2*u31*x2/2+u2*u31*v1*x2/2+c31*u1*u2*x1*x2/2+
  c1*u2*u31*x1*x2/2+c2*c31*u1*x3/2+c1*c31*u2*x3/2+
  c1*c2*u31*x3/2+u1*u31*v2*x3/2+c31*u1*u2*x1*x3/2+
  c2*u1*u31*x1*x3/2+c2*u1*u31*x2*x3/2+c1*u2*u31*x2*x3/2+
  u1*u2*u31*x1*x2*x3/2;
 d110psi=c11*c21*c3+c3*u21*v11/2+c21*u3*v11/2+c3*u11*v21/2+
  c11*u3*v21/2+c21*u11*v3/2+c11*u21*v3/2+c21*c3*u11*x1/2+
  c11*c3*u21*x1/2+c11*c21*u3*x1/2+u11*u21*v3*x1/2+
  c21*c3*u11*x2/2+c11*c3*u21*x2/2+c11*c21*u3*x2/2+
  u21*u3*v11*x2/2+c3*u11*u21*x1*x2/2+c11*u21*u3*x1*x2/2+
  c21*c3*u11*x3/2+c11*c3*u21*x3/2+c11*c21*u3*x3/2+
  u11*u3*v21*x3/2+c3*u11*u21*x1*x3/2+c21*u11*u3*x1*x3/2+
  c21*u11*u3*x2*x3/2+c11*u21*u3*x2*x3/2+u11*u21*u3*x1*x2*x3/2;
 d200psi=c12*c2*c3+c3*u2*v12/2+c2*u3*v12/2+c3*u12*v2/2+c12*u3*v2/2+
  c2*u12*v3/2+c12*u2*v3/2+c2*c3*u12*x1/2+c12*c3*u2*x1/2+
  c12*c2*u3*x1/2+u12*u2*v3*x1/2+c2*c3*u12*x2/2+c12*c3*u2*x2/2+
  c12*c2*u3*x2/2+u2*u3*v12*x2/2+c3*u12*u2*x1*x2/2+
  c12*u2*u3*x1*x2/2+c2*c3*u12*x3/2+c12*c3*u2*x3/2+
  c12*c2*u3*x3/2+u12*u3*v2*x3/2+c3*u12*u2*x1*x3/2+
  c2*u12*u3*x1*x3/2+c2*u12*u3*x2*x3/2+c12*u2*u3*x2*x3/2+
  u12*u2*u3*x1*x2*x3/2;
 d020psi=c1*c22*c3+c3*u22*v1/2+c22*u3*v1/2+c3*u1*v22/2+c1*u3*v22/2+
  c22*u1*v3/2+c1*u22*v3/2+c22*c3*u1*x1/2+c1*c3*u22*x1/2+
  c1*c22*u3*x1/2+u1*u22*v3*x1/2+c22*c3*u1*x2/2+c1*c3*u22*x2/2+
  c1*c22*u3*x2/2+u22*u3*v1*x2/2+c3*u1*u22*x1*x2/2+
  c1*u22*u3*x1*x2/2+c22*c3*u1*x3/2+c1*c3*u22*x3/2+
  c1*c22*u3*x3/2+u1*u3*v22*x3/2+c3*u1*u22*x1*x3/2+
  c22*u1*u3*x1*x3/2+c22*u1*u3*x2*x3/2+c1*u22*u3*x2*x3/2+
  u1*u22*u3*x1*x2*x3/2;
 d101psi=c11*c2*c31+c31*u2*v11/2+c2*u31*v11/2+c31*u11*v2/2+
  c11*u31*v2/2+c2*u11*v31/2+c11*u2*v31/2+c2*c31*u11*x1/2+
  c11*c31*u2*x1/2+c11*c2*u31*x1/2+u11*u2*v31*x1/2+
  c2*c31*u11*x2/2+c11*c31*u2*x2/2+c11*c2*u31*x2/2+
  u2*u31*v11*x2/2+c31*u11*u2*x1*x2/2+c11*u2*u31*x1*x2/2+
  c2*c31*u11*x3/2+c11*c31*u2*x3/2+c11*c2*u31*x3/2+
  u11*u31*v2*x3/2+c31*u11*u2*x1*x3/2+c2*u11*u31*x1*x3/2+
  c2*u11*u31*x2*x3/2+c11*u2*u31*x2*x3/2+u11*u2*u31*x1*x2*x3/2;
 d011psi=c1*c21*c31+c31*u21*v1/2+c21*u31*v1/2+c31*u1*v21/2+
  c1*u31*v21/2+c21*u1*v31/2+c1*u21*v31/2+c21*c31*u1*x1/2+
  c1*c31*u21*x1/2+c1*c21*u31*x1/2+u1*u21*v31*x1/2+
  c21*c31*u1*x2/2+c1*c31*u21*x2/2+c1*c21*u31*x2/2+
  u21*u31*v1*x2/2+c31*u1*u21*x1*x2/2+c1*u21*u31*x1*x2/2+
  c21*c31*u1*x3/2+c1*c31*u21*x3/2+c1*c21*u31*x3/2+
  u1*u31*v21*x3/2+c31*u1*u21*x1*x3/2+c21*u1*u31*x1*x3/2+
  c21*u1*u31*x2*x3/2+c1*u21*u31*x2*x3/2+u1*u21*u31*x1*x2*x3/2;
 d002psi=c1*c2*c32+c32*u2*v1/2+c2*u32*v1/2+c32*u1*v2/2+c1*u32*v2/2+
  c2*u1*v32/2+c1*u2*v32/2+c2*c32*u1*x1/2+c1*c32*u2*x1/2+
  c1*c2*u32*x1/2+u1*u2*v32*x1/2+c2*c32*u1*x2/2+c1*c32*u2*x2/2+
  c1*c2*u32*x2/2+u2*u32*v1*x2/2+c32*u1*u2*x1*x2/2+
  c1*u2*u32*x1*x2/2+c2*c32*u1*x3/2+c1*c32*u2*x3/2+
  c1*c2*u32*x3/2+u1*u32*v2*x3/2+c32*u1*u2*x1*x3/2+
  c2*u1*u32*x1*x3/2+c2*u1*u32*x2*x3/2+c1*u2*u32*x2*x3/2+
  u1*u2*u32*x1*x2*x3/2;
 d300psi=c13*c2*c3+c3*u2*v13/2+c2*u3*v13/2+c3*u13*v2/2+c13*u3*v2/2+
  c2*u13*v3/2+c13*u2*v3/2+c2*c3*u13*x1/2+c13*c3*u2*x1/2+
  c13*c2*u3*x1/2+u13*u2*v3*x1/2+c2*c3*u13*x2/2+c13*c3*u2*x2/2+
  c13*c2*u3*x2/2+u2*u3*v13*x2/2+c3*u13*u2*x1*x2/2+
  c13*u2*u3*x1*x2/2+c2*c3*u13*x3/2+c13*c3*u2*x3/2+
  c13*c2*u3*x3/2+u13*u3*v2*x3/2+c3*u13*u2*x1*x3/2+
  c2*u13*u3*x1*x3/2+c2*u13*u3*x2*x3/2+c13*u2*u3*x2*x3/2+
  u13*u2*u3*x1*x2*x3/2;
 d030psi=c1*c23*c3+c3*u23*v1/2+c23*u3*v1/2+c3*u1*v23/2+c1*u3*v23/2+
  c23*u1*v3/2+c1*u23*v3/2+c23*c3*u1*x1/2+c1*c3*u23*x1/2+
  c1*c23*u3*x1/2+u1*u23*v3*x1/2+c23*c3*u1*x2/2+c1*c3*u23*x2/2+
  c1*c23*u3*x2/2+u23*u3*v1*x2/2+c3*u1*u23*x1*x2/2+
  c1*u23*u3*x1*x2/2+c23*c3*u1*x3/2+c1*c3*u23*x3/2+
  c1*c23*u3*x3/2+u1*u3*v23*x3/2+c3*u1*u23*x1*x3/2+
  c23*u1*u3*x1*x3/2+c23*u1*u3*x2*x3/2+c1*u23*u3*x2*x3/2+
  u1*u23*u3*x1*x2*x3/2;
 d003psi=c1*c2*c33+c33*u2*v1/2+c2*u33*v1/2+c33*u1*v2/2+c1*u33*v2/2+
  c2*u1*v33/2+c1*u2*v33/2+c2*c33*u1*x1/2+c1*c33*u2*x1/2+
  c1*c2*u33*x1/2+u1*u2*v33*x1/2+c2*c33*u1*x2/2+c1*c33*u2*x2/2+
  c1*c2*u33*x2/2+u2*u33*v1*x2/2+c33*u1*u2*x1*x2/2+
  c1*u2*u33*x1*x2/2+c2*c33*u1*x3/2+c1*c33*u2*x3/2+
  c1*c2*u33*x3/2+u1*u33*v2*x3/2+c33*u1*u2*x1*x3/2+
  c2*u1*u33*x1*x3/2+c2*u1*u33*x2*x3/2+c1*u2*u33*x2*x3/2+
  u1*u2*u33*x1*x2*x3/2;
 d120psi=c11*c22*c3+c3*u22*v11/2+c22*u3*v11/2+c3*u11*v22/2+
  c11*u3*v22/2+c22*u11*v3/2+c11*u22*v3/2+c22*c3*u11*x1/2+
  c11*c3*u22*x1/2+c11*c22*u3*x1/2+u11*u22*v3*x1/2+
  c22*c3*u11*x2/2+c11*c3*u22*x2/2+c11*c22*u3*x2/2+
  u22*u3*v11*x2/2+c3*u11*u22*x1*x2/2+c11*u22*u3*x1*x2/2+
  c22*c3*u11*x3/2+c11*c3*u22*x3/2+c11*c22*u3*x3/2+
  u11*u3*v22*x3/2+c3*u11*u22*x1*x3/2+c22*u11*u3*x1*x3/2+
  c22*u11*u3*x2*x3/2+c11*u22*u3*x2*x3/2+u11*u22*u3*x1*x2*x3/2;
 d210psi=c12*c21*c3+c3*u21*v12/2+c21*u3*v12/2+c3*u12*v21/2+
  c12*u3*v21/2+c21*u12*v3/2+c12*u21*v3/2+c21*c3*u12*x1/2+
  c12*c3*u21*x1/2+c12*c21*u3*x1/2+u12*u21*v3*x1/2+
  c21*c3*u12*x2/2+c12*c3*u21*x2/2+c12*c21*u3*x2/2+
  u21*u3*v12*x2/2+c3*u12*u21*x1*x2/2+c12*u21*u3*x1*x2/2+
  c21*c3*u12*x3/2+c12*c3*u21*x3/2+c12*c21*u3*x3/2+
  u12*u3*v21*x3/2+c3*u12*u21*x1*x3/2+c21*u12*u3*x1*x3/2+
  c21*u12*u3*x2*x3/2+c12*u21*u3*x2*x3/2+u12*u21*u3*x1*x2*x3/2;
 d102psi=c11*c2*c32+c32*u2*v11/2+c2*u32*v11/2+c32*u11*v2/2+
  c11*u32*v2/2+c2*u11*v32/2+c11*u2*v32/2+c2*c32*u11*x1/2+
  c11*c32*u2*x1/2+c11*c2*u32*x1/2+u11*u2*v32*x1/2+
  c2*c32*u11*x2/2+c11*c32*u2*x2/2+c11*c2*u32*x2/2+
  u2*u32*v11*x2/2+c32*u11*u2*x1*x2/2+c11*u2*u32*x1*x2/2+
  c2*c32*u11*x3/2+c11*c32*u2*x3/2+c11*c2*u32*x3/2+
  u11*u32*v2*x3/2+c32*u11*u2*x1*x3/2+c2*u11*u32*x1*x3/2+
  c2*u11*u32*x2*x3/2+c11*u2*u32*x2*x3/2+u11*u2*u32*x1*x2*x3/2;
 d201psi=c12*c2*c31+c31*u2*v12/2+c2*u31*v12/2+c31*u12*v2/2+
  c12*u31*v2/2+c2*u12*v31/2+c12*u2*v31/2+c2*c31*u12*x1/2+
  c12*c31*u2*x1/2+c12*c2*u31*x1/2+u12*u2*v31*x1/2+
  c2*c31*u12*x2/2+c12*c31*u2*x2/2+c12*c2*u31*x2/2+
  u2*u31*v12*x2/2+c31*u12*u2*x1*x2/2+c12*u2*u31*x1*x2/2+
  c2*c31*u12*x3/2+c12*c31*u2*x3/2+c12*c2*u31*x3/2+
  u12*u31*v2*x3/2+c31*u12*u2*x1*x3/2+c2*u12*u31*x1*x3/2+
  c2*u12*u31*x2*x3/2+c12*u2*u31*x2*x3/2+u12*u2*u31*x1*x2*x3/2;
 d012psi=c1*c21*c32+c32*u21*v1/2+c21*u32*v1/2+c32*u1*v21/2+
  c1*u32*v21/2+c21*u1*v32/2+c1*u21*v32/2+c21*c32*u1*x1/2+
  c1*c32*u21*x1/2+c1*c21*u32*x1/2+u1*u21*v32*x1/2+
  c21*c32*u1*x2/2+c1*c32*u21*x2/2+c1*c21*u32*x2/2+
  u21*u32*v1*x2/2+c32*u1*u21*x1*x2/2+c1*u21*u32*x1*x2/2+
  c21*c32*u1*x3/2+c1*c32*u21*x3/2+c1*c21*u32*x3/2+
  u1*u32*v21*x3/2+c32*u1*u21*x1*x3/2+c21*u1*u32*x1*x3/2+
  c21*u1*u32*x2*x3/2+c1*u21*u32*x2*x3/2+u1*u21*u32*x1*x2*x3/2;
 d021psi=c1*c22*c31+c31*u22*v1/2+c22*u31*v1/2+c31*u1*v22/2+
  c1*u31*v22/2+c22*u1*v31/2+c1*u22*v31/2+c22*c31*u1*x1/2+
  c1*c31*u22*x1/2+c1*c22*u31*x1/2+u1*u22*v31*x1/2+
  c22*c31*u1*x2/2+c1*c31*u22*x2/2+c1*c22*u31*x2/2+
  u22*u31*v1*x2/2+c31*u1*u22*x1*x2/2+c1*u22*u31*x1*x2/2+
  c22*c31*u1*x3/2+c1*c31*u22*x3/2+c1*c22*u31*x3/2+
  u1*u31*v22*x3/2+c31*u1*u22*x1*x3/2+c22*u1*u31*x1*x3/2+
  c22*u1*u31*x2*x3/2+c1*u22*u31*x2*x3/2+u1*u22*u31*x1*x2*x3/2;
 d111psi=c11*c21*c31+c31*u21*v11/2+c21*u31*v11/2+c31*u11*v21/2+
  c11*u31*v21/2+c21*u11*v31/2+c11*u21*v31/2+c21*c31*u11*x1/2+
  c11*c31*u21*x1/2+c11*c21*u31*x1/2+u11*u21*v31*x1/2+
  c21*c31*u11*x2/2+c11*c31*u21*x2/2+c11*c21*u31*x2/2+
  u21*u31*v11*x2/2+c31*u11*u21*x1*x2/2+c11*u21*u31*x1*x2/2+
  c21*c31*u11*x3/2+c11*c31*u21*x3/2+c11*c21*u31*x3/2+
  u11*u31*v21*x3/2+c31*u11*u21*x1*x3/2+c21*u11*u31*x1*x3/2+
  c21*u11*u31*x2*x3/2+c11*u21*u31*x2*x3/2+u11*u21*u31*x1*x2*x3/2;
 d400psi=c14*c2*c3+c3*u2*v14/2+c2*u3*v14/2+c3*u14*v2/2+c14*u3*v2/2+
  c2*u14*v3/2+c14*u2*v3/2+c2*c3*u14*x1/2+c14*c3*u2*x1/2+
  c14*c2*u3*x1/2+u14*u2*v3*x1/2+c2*c3*u14*x2/2+c14*c3*u2*x2/2+
  c14*c2*u3*x2/2+u2*u3*v14*x2/2+c3*u14*u2*x1*x2/2+
  c14*u2*u3*x1*x2/2+c2*c3*u14*x3/2+c14*c3*u2*x3/2+
  c14*c2*u3*x3/2+u14*u3*v2*x3/2+c3*u14*u2*x1*x3/2+
  c2*u14*u3*x1*x3/2+c2*u14*u3*x2*x3/2+c14*u2*u3*x2*x3/2+
  u14*u2*u3*x1*x2*x3/2;
 d040psi=c1*c24*c3+c3*u24*v1/2+c24*u3*v1/2+c3*u1*v24/2+c1*u3*v24/2+
  c24*u1*v3/2+c1*u24*v3/2+c24*c3*u1*x1/2+c1*c3*u24*x1/2+
  c1*c24*u3*x1/2+u1*u24*v3*x1/2+c24*c3*u1*x2/2+c1*c3*u24*x2/2+
  c1*c24*u3*x2/2+u24*u3*v1*x2/2+c3*u1*u24*x1*x2/2+
  c1*u24*u3*x1*x2/2+c24*c3*u1*x3/2+c1*c3*u24*x3/2+
  c1*c24*u3*x3/2+u1*u3*v24*x3/2+c3*u1*u24*x1*x3/2+
  c24*u1*u3*x1*x3/2+c24*u1*u3*x2*x3/2+c1*u24*u3*x2*x3/2+
  u1*u24*u3*x1*x2*x3/2;
 d004psi=c1*c2*c34+c34*u2*v1/2+c2*u34*v1/2+c34*u1*v2/2+c1*u34*v2/2+
  c2*u1*v34/2+c1*u2*v34/2+c2*c34*u1*x1/2+c1*c34*u2*x1/2+
  c1*c2*u34*x1/2+u1*u2*v34*x1/2+c2*c34*u1*x2/2+c1*c34*u2*x2/2+
  c1*c2*u34*x2/2+u2*u34*v1*x2/2+c34*u1*u2*x1*x2/2+
  c1*u2*u34*x1*x2/2+c2*c34*u1*x3/2+c1*c34*u2*x3/2+
  c1*c2*u34*x3/2+u1*u34*v2*x3/2+c34*u1*u2*x1*x3/2+
  c2*u1*u34*x1*x3/2+c2*u1*u34*x2*x3/2+c1*u2*u34*x2*x3/2+
  u1*u2*u34*x1*x2*x3/2;
 d130psi=c11*c23*c3+c3*u23*v11/2+c23*u3*v11/2+c3*u11*v23/2+
  c11*u3*v23/2+c23*u11*v3/2+c11*u23*v3/2+c23*c3*u11*x1/2+
  c11*c3*u23*x1/2+c11*c23*u3*x1/2+u11*u23*v3*x1/2+
  c23*c3*u11*x2/2+c11*c3*u23*x2/2+c11*c23*u3*x2/2+
  u23*u3*v11*x2/2+c3*u11*u23*x1*x2/2+c11*u23*u3*x1*x2/2+
  c23*c3*u11*x3/2+c11*c3*u23*x3/2+c11*c23*u3*x3/2+
  u11*u3*v23*x3/2+c3*u11*u23*x1*x3/2+c23*u11*u3*x1*x3/2+
  c23*u11*u3*x2*x3/2+c11*u23*u3*x2*x3/2+u11*u23*u3*x1*x2*x3/2;
 d310psi=c13*c21*c3+c3*u21*v13/2+c21*u3*v13/2+c3*u13*v21/2+
  c13*u3*v21/2+c21*u13*v3/2+c13*u21*v3/2+c21*c3*u13*x1/2+
  c13*c3*u21*x1/2+c13*c21*u3*x1/2+u13*u21*v3*x1/2+
  c21*c3*u13*x2/2+c13*c3*u21*x2/2+c13*c21*u3*x2/2+
  u21*u3*v13*x2/2+c3*u13*u21*x1*x2/2+c13*u21*u3*x1*x2/2+
  c21*c3*u13*x3/2+c13*c3*u21*x3/2+c13*c21*u3*x3/2+
  u13*u3*v21*x3/2+c3*u13*u21*x1*x3/2+c21*u13*u3*x1*x3/2+
  c21*u13*u3*x2*x3/2+c13*u21*u3*x2*x3/2+u13*u21*u3*x1*x2*x3/2;
 d220psi=c12*c22*c3+c3*u22*v12/2+c22*u3*v12/2+c3*u12*v22/2+
  c12*u3*v22/2+c22*u12*v3/2+c12*u22*v3/2+c22*c3*u12*x1/2+
  c12*c3*u22*x1/2+c12*c22*u3*x1/2+u12*u22*v3*x1/2+
  c22*c3*u12*x2/2+c12*c3*u22*x2/2+c12*c22*u3*x2/2+
  u22*u3*v12*x2/2+c3*u12*u22*x1*x2/2+c12*u22*u3*x1*x2/2+
  c22*c3*u12*x3/2+c12*c3*u22*x3/2+c12*c22*u3*x3/2+
  u12*u3*v22*x3/2+c3*u12*u22*x1*x3/2+c22*u12*u3*x1*x3/2+
  c22*u12*u3*x2*x3/2+c12*u22*u3*x2*x3/2+u12*u22*u3*x1*x2*x3/2;
 d103psi=c11*c2*c33+c33*u2*v11/2+c2*u33*v11/2+c33*u11*v2/2+
  c11*u33*v2/2+c2*u11*v33/2+c11*u2*v33/2+c2*c33*u11*x1/2+
  c11*c33*u2*x1/2+c11*c2*u33*x1/2+u11*u2*v33*x1/2+
  c2*c33*u11*x2/2+c11*c33*u2*x2/2+c11*c2*u33*x2/2+
  u2*u33*v11*x2/2+c33*u11*u2*x1*x2/2+c11*u2*u33*x1*x2/2+
  c2*c33*u11*x3/2+c11*c33*u2*x3/2+c11*c2*u33*x3/2+
  u11*u33*v2*x3/2+c33*u11*u2*x1*x3/2+c2*u11*u33*x1*x3/2+
  c2*u11*u33*x2*x3/2+c11*u2*u33*x2*x3/2+u11*u2*u33*x1*x2*x3/2;
 d301psi=c13*c2*c31+c31*u2*v13/2+c2*u31*v13/2+c31*u13*v2/2+
  c13*u31*v2/2+c2*u13*v31/2+c13*u2*v31/2+c2*c31*u13*x1/2+
  c13*c31*u2*x1/2+c13*c2*u31*x1/2+u13*u2*v31*x1/2+
  c2*c31*u13*x2/2+c13*c31*u2*x2/2+c13*c2*u31*x2/2+
  u2*u31*v13*x2/2+c31*u13*u2*x1*x2/2+c13*u2*u31*x1*x2/2+
  c2*c31*u13*x3/2+c13*c31*u2*x3/2+c13*c2*u31*x3/2+
  u13*u31*v2*x3/2+c31*u13*u2*x1*x3/2+c2*u13*u31*x1*x3/2+
  c2*u13*u31*x2*x3/2+c13*u2*u31*x2*x3/2+u13*u2*u31*x1*x2*x3/2;;
 d202psi=c12*c2*c32+c32*u2*v12/2+c2*u32*v12/2+c32*u12*v2/2+
  c12*u32*v2/2+c2*u12*v32/2+c12*u2*v32/2+c2*c32*u12*x1/2+
  c12*c32*u2*x1/2+c12*c2*u32*x1/2+u12*u2*v32*x1/2+
  c2*c32*u12*x2/2+c12*c32*u2*x2/2+c12*c2*u32*x2/2+
  u2*u32*v12*x2/2+c32*u12*u2*x1*x2/2+c12*u2*u32*x1*x2/2+
  c2*c32*u12*x3/2+c12*c32*u2*x3/2+c12*c2*u32*x3/2+
  u12*u32*v2*x3/2+c32*u12*u2*x1*x3/2+c2*u12*u32*x1*x3/2+
  c2*u12*u32*x2*x3/2+c12*u2*u32*x2*x3/2+u12*u2*u32*x1*x2*x3/2;
 d013psi=c1*c21*c33+c33*u21*v1/2+c21*u33*v1/2+c33*u1*v21/2+
  c1*u33*v21/2+c21*u1*v33/2+c1*u21*v33/2+c21*c33*u1*x1/2+
  c1*c33*u21*x1/2+c1*c21*u33*x1/2+u1*u21*v33*x1/2+
  c21*c33*u1*x2/2+c1*c33*u21*x2/2+c1*c21*u33*x2/2+
  u21*u33*v1*x2/2+c33*u1*u21*x1*x2/2+c1*u21*u33*x1*x2/2+
  c21*c33*u1*x3/2+c1*c33*u21*x3/2+c1*c21*u33*x3/2+
  u1*u33*v21*x3/2+c33*u1*u21*x1*x3/2+c21*u1*u33*x1*x3/2+
  c21*u1*u33*x2*x3/2+c1*u21*u33*x2*x3/2+u1*u21*u33*x1*x2*x3/2;
 d031psi=c1*c23*c31+c31*u23*v1/2+c23*u31*v1/2+c31*u1*v23/2+
  c1*u31*v23/2+c23*u1*v31/2+c1*u23*v31/2+c23*c31*u1*x1/2+
  c1*c31*u23*x1/2+c1*c23*u31*x1/2+u1*u23*v31*x1/2+
  c23*c31*u1*x2/2+c1*c31*u23*x2/2+c1*c23*u31*x2/2+
  u23*u31*v1*x2/2+c31*u1*u23*x1*x2/2+c1*u23*u31*x1*x2/2+
  c23*c31*u1*x3/2+c1*c31*u23*x3/2+c1*c23*u31*x3/2+
  u1*u31*v23*x3/2+c31*u1*u23*x1*x3/2+c23*u1*u31*x1*x3/2+
  c23*u1*u31*x2*x3/2+c1*u23*u31*x2*x3/2+u1*u23*u31*x1*x2*x3/2;
 d022psi=c1*c22*c32+c32*u22*v1/2+c22*u32*v1/2+c32*u1*v22/2+
  c1*u32*v22/2+c22*u1*v32/2+c1*u22*v32/2+c22*c32*u1*x1/2+
  c1*c32*u22*x1/2+c1*c22*u32*x1/2+u1*u22*v32*x1/2+
  c22*c32*u1*x2/2+c1*c32*u22*x2/2+c1*c22*u32*x2/2+
  u22*u32*v1*x2/2+c32*u1*u22*x1*x2/2+c1*u22*u32*x1*x2/2+
  c22*c32*u1*x3/2+c1*c32*u22*x3/2+c1*c22*u32*x3/2+
  u1*u32*v22*x3/2+c32*u1*u22*x1*x3/2+c22*u1*u32*x1*x3/2+
  c22*u1*u32*x2*x3/2+c1*u22*u32*x2*x3/2+u1*u22*u32*x1*x2*x3/2;
 d211psi=c12*c21*c31+c31*u21*v12/2+c21*u31*v12/2+c31*u12*v21/2+
  c12*u31*v21/2+c21*u12*v31/2+c12*u21*v31/2+c21*c31*u12*x1/2+
  c12*c31*u21*x1/2+c12*c21*u31*x1/2+u12*u21*v31*x1/2+
  c21*c31*u12*x2/2+c12*c31*u21*x2/2+c12*c21*u31*x2/2+
  u21*u31*v12*x2/2+c31*u12*u21*x1*x2/2+c12*u21*u31*x1*x2/2+
  c21*c31*u12*x3/2+c12*c31*u21*x3/2+c12*c21*u31*x3/2+
  u12*u31*v21*x3/2+c31*u12*u21*x1*x3/2+c21*u12*u31*x1*x3/2+
  c21*u12*u31*x2*x3/2+c12*u21*u31*x2*x3/2+u12*u21*u31*x1*x2*x3/2;
 d121psi=c11*c22*c31+c31*u22*v11/2+c22*u31*v11/2+c31*u11*v22/2+
  c11*u31*v22/2+c22*u11*v31/2+c11*u22*v31/2+c22*c31*u11*x1/2+
  c11*c31*u22*x1/2+c11*c22*u31*x1/2+u11*u22*v31*x1/2+
  c22*c31*u11*x2/2+c11*c31*u22*x2/2+c11*c22*u31*x2/2+
  u22*u31*v11*x2/2+c31*u11*u22*x1*x2/2+c11*u22*u31*x1*x2/2+
  c22*c31*u11*x3/2+c11*c31*u22*x3/2+c11*c22*u31*x3/2+
  u11*u31*v22*x3/2+c31*u11*u22*x1*x3/2+c22*u11*u31*x1*x3/2+
  c22*u11*u31*x2*x3/2+c11*u22*u31*x2*x3/2+u11*u22*u31*x1*x2*x3/2;
 d112psi=c11*c21*c32+c32*u21*v11/2+c21*u32*v11/2+c32*u11*v21/2+
  c11*u32*v21/2+c21*u11*v32/2+c11*u21*v32/2+c21*c32*u11*x1/2+
  c11*c32*u21*x1/2+c11*c21*u32*x1/2+u11*u21*v32*x1/2+
  c21*c32*u11*x2/2+c11*c32*u21*x2/2+c11*c21*u32*x2/2+
  u21*u32*v11*x2/2+c32*u11*u21*x1*x2/2+c11*u21*u32*x1*x2/2+
  c21*c32*u11*x3/2+c11*c32*u21*x3/2+c11*c21*u32*x3/2+
  u11*u32*v21*x3/2+c32*u11*u21*x1*x3/2+c21*u11*u32*x1*x3/2+
  c21*u11*u32*x2*x3/2+c11*u21*u32*x2*x3/2+u11*u21*u32*x1*x2*x3/2;

 psi=psih-c0;
 d100psi2=d100psi*d100psi;
 d010psi2=d010psi*d010psi;
 d001psi2=d001psi*d001psi;
 d200psi2=d200psi*d200psi;
 d020psi2=d020psi*d020psi;
 d002psi2=d002psi*d002psi;
 d110psi2=d110psi*d110psi;
 d011psi2=d011psi*d011psi;
 d101psi2=d101psi*d101psi;

 L1psi=d200psi+2*d100psi/r1;
 L2psi=d020psi+2*d010psi/r2;
 L3psi=d002psi+2*d001psi/r3;
 Lpsi=L1psi+L2psi+L3psi+2*r1r2h*d110psi+2*r1r3h*d101psi+2*r2r3h*d011psi;
     
 dApsi=Lpsi+4*Pi*Pi*c0;
 dipsi2=d100psi*d100psi+d010psi*d010psi+2*r1r2h*d100psi*d010psi+
  d001psi*d001psi+2*r1r3h*d100psi*d001psi+2*r2r3h*d010psi*d001psi;
 dpsi2=dipsi2+4*Pi*Pi*s02;
     
/*     dAlnpsi=dApsi/psi-dpsi2/psi/psi;  */

 d00dAlnpsi  =(Lpsi-4*Pi*Pi*psih)*(2*s02/psi-c0);
 d00dAlnpsi -=2*(dipsi2+4*Pi*Pi*(1-psih*psih))*(3*s02/psi-c0)/psi;
 d00dAlnpsi *=4*Pi*Pi/psi/psi;

 LLpsi=d004psi+2*d022psi+d040psi+2*d202psi+2*d220psi+d400psi+
  (8*d102psi)/r1+(8*d120psi)/r1+(4*d300psi)/r1+4*d112psi*r1r2h+
  4*d130psi*r1r2h+4*d310psi*r1r2h-(8*d110psi*r1r2h)/r12+
  (8*d210psi*r1r2h)/r1+4*d220psi*r1r2h2-(4*d120psi*r1r2h2)/r1+
  4*d103psi*r1r3h+4*d121psi*r1r3h+4*d301psi*r1r3h-
  (8*d101psi*r1r3h)/r12+(8*d201psi*r1r3h)/r1+8*d211psi*r1r2h*r1r3h-
  (8*d111psi*r1r2h*r1r3h)/r1+4*d202psi*r1r3h2-(4*d102psi*r1r3h2)/r1-
  (8*d110psi*r1r2h)/r22+(8*d012psi)/r2+(4*d030psi)/r2+(8*d210psi)/r2+
  (12*d110psi)/(r1*r2)+(8*d120psi*r1r2h)/r2-(4*d210psi*r1r2h2)/r2+
  (4*d110psi*r1r2h2)/(r1*r2)+(16*d111psi*r1r3h)/r2+4*d013psi*r2r3h+
  4*d031psi*r2r3h+4*d211psi*r2r3h+(16*d111psi*r2r3h)/r1+
  8*d121psi*r1r2h*r2r3h+8*d112psi*r1r3h*r2r3h-(8*d011psi*r2r3h)/r22+
  (8*d021psi*r2r3h)/r2-(8*d111psi*r1r2h*r2r3h)/r2+4*d022psi*r2r3h2-
  (4*d012psi*r2r3h2)/r2-(8*d101psi*r1r3h)/r32-(8*d011psi*r2r3h)/r32+
  (4*d003psi)/r3+(8*d021psi)/r3+(8*d201psi)/r3+(12*d101psi)/(r1*r3)+
  (16*d111psi*r1r2h)/r3+(8*d102psi*r1r3h)/r3-(4*d201psi*r1r3h2)/r3+
  (4*d101psi*r1r3h2)/(r1*r3)+(12*d011psi)/(r2*r3)+(8*d012psi*r2r3h)/r3-
  (8*d111psi*r1r3h*r2r3h)/r3-(4*d021psi*r2r3h2)/r3+
  (4*d011psi*r2r3h2)/(r2*r3);
 Ldpsi2p2dpsidLpsi=2*d002psi2+4*d001psi*d003psi+4*d011psi2+4*d010psi*d012psi+
  2*d020psi2+4*d001psi*d021psi+4*d010psi*d030psi+4*d101psi2+
  4*d100psi*d102psi+4*d110psi2+4*d100psi*d120psi+2*d200psi2+
  4*d001psi*d201psi+4*d010psi*d210psi+4*d100psi*d300psi-
  (4*d100psi2)/r12+(4*d002psi*d100psi)/r1+(4*d020psi*d100psi)/r1+
  (16*d001psi*d101psi)/r1+(16*d010psi*d110psi)/r1+
  (8*d100psi*d200psi)/r1+4*d012psi*d100psi*r1r2h+
  4*d030psi*d100psi*r1r2h+8*d011psi*d101psi*r1r2h+
  4*d010psi*d102psi*r1r2h+8*d020psi*d110psi*r1r2h+
  8*d001psi*d111psi*r1r2h+12*d010psi*d120psi*r1r2h+
  8*d110psi*d200psi*r1r2h+12*d100psi*d210psi*r1r2h+
  4*d010psi*d300psi*r1r2h-(8*d010psi*d100psi*r1r2h)/r12+
  (8*d100psi*d110psi*r1r2h)/r1+(8*d010psi*d200psi*r1r2h)/r1+
  4*d110psi2*r1r2h2+8*d100psi*d120psi*r1r2h2+
  4*d020psi*d200psi*r1r2h2+8*d010psi*d210psi*r1r2h2-
  (4*d020psi*d100psi*r1r2h2)/r1-(8*d010psi*d110psi*r1r2h2)/r1+
  4*d003psi*d100psi*r1r3h+4*d021psi*d100psi*r1r3h+
  8*d002psi*d101psi*r1r3h+12*d001psi*d102psi*r1r3h+
  8*d011psi*d110psi*r1r3h+8*d010psi*d111psi*r1r3h+
  4*d001psi*d120psi*r1r3h+8*d101psi*d200psi*r1r3h+
  12*d100psi*d201psi*r1r3h+4*d001psi*d300psi*r1r3h-
  (8*d001psi*d100psi*r1r3h)/r12+(8*d100psi*d101psi*r1r3h)/r1+
  (8*d001psi*d200psi*r1r3h)/r1+8*d101psi*d110psi*r1r2h*r1r3h+
  16*d100psi*d111psi*r1r2h*r1r3h+8*d011psi*d200psi*r1r2h*r1r3h+
  8*d010psi*d201psi*r1r2h*r1r3h+8*d001psi*d210psi*r1r2h*r1r3h-
  (8*d011psi*d100psi*r1r2h*r1r3h)/r1-(8*d010psi*d101psi*r1r2h*r1r3h)/r1-
  (8*d001psi*d110psi*r1r2h*r1r3h)/r1+4*d101psi2*r1r3h2+
  8*d100psi*d102psi*r1r3h2+4*d002psi*d200psi*r1r3h2+
  8*d001psi*d201psi*r1r3h2-(4*d002psi*d100psi*r1r3h2)/r1-
  (8*d001psi*d101psi*r1r3h2)/r1-(4*d010psi2)/r22-
  (8*d010psi*d100psi*r1r2h)/r22+(4*d002psi*d010psi)/r2+
  (16*d001psi*d011psi)/r2+(8*d010psi*d020psi)/r2+
  (16*d100psi*d110psi)/r2+(4*d010psi*d200psi)/r2+
  (4*d010psi*d100psi)/(r1*r2)+(8*d020psi*d100psi*r1r2h)/r2+
  (8*d010psi*d110psi*r1r2h)/r2-(8*d100psi*d110psi*r1r2h2)/r2-
  (4*d010psi*d200psi*r1r2h2)/r2+(4*d010psi*d100psi*r1r2h2)/(r1*r2)+
  (16*d011psi*d100psi*r1r3h)/r2+(8*d010psi*d101psi*r1r3h)/r2+
  (16*d001psi*d110psi*r1r3h)/r2+4*d003psi*d010psi*r2r3h+
  8*d002psi*d011psi*r2r3h+12*d001psi*d012psi*r2r3h+
  8*d011psi*d020psi*r2r3h+12*d010psi*d021psi*r2r3h+
  4*d001psi*d030psi*r2r3h+8*d101psi*d110psi*r2r3h+
  8*d100psi*d111psi*r2r3h+4*d010psi*d201psi*r2r3h+
  4*d001psi*d210psi*r2r3h+(8*d011psi*d100psi*r2r3h)/r1+
  (16*d010psi*d101psi*r2r3h)/r1+(16*d001psi*d110psi*r2r3h)/r1+
  8*d021psi*d100psi*r1r2h*r2r3h+8*d020psi*d101psi*r1r2h*r2r3h+
  8*d011psi*d110psi*r1r2h*r2r3h+16*d010psi*d111psi*r1r2h*r2r3h+
  8*d001psi*d120psi*r1r2h*r2r3h+8*d012psi*d100psi*r1r3h*r2r3h+
  8*d011psi*d101psi*r1r3h*r2r3h+8*d010psi*d102psi*r1r3h*r2r3h+
  8*d002psi*d110psi*r1r3h*r2r3h+16*d001psi*d111psi*r1r3h*r2r3h-
  (8*d001psi*d010psi*r2r3h)/r22+(8*d010psi*d011psi*r2r3h)/r2+
  (8*d001psi*d020psi*r2r3h)/r2-(8*d011psi*d100psi*r1r2h*r2r3h)/r2-
  (8*d010psi*d101psi*r1r2h*r2r3h)/r2-(8*d001psi*d110psi*r1r2h*r2r3h)/r2+
  4*d011psi2*r2r3h2+8*d010psi*d012psi*r2r3h2+
  4*d002psi*d020psi*r2r3h2+8*d001psi*d021psi*r2r3h2-
  (4*d002psi*d010psi*r2r3h2)/r2-(8*d001psi*d011psi*r2r3h2)/r2-
  (4*d001psi2)/r32-(8*d001psi*d100psi*r1r3h)/r32-
  (8*d001psi*d010psi*r2r3h)/r32+(8*d001psi*d002psi)/r3+
  (16*d010psi*d011psi)/r3+(4*d001psi*d020psi)/r3+
  (16*d100psi*d101psi)/r3+(4*d001psi*d200psi)/r3+
  (4*d001psi*d100psi)/(r1*r3)+(16*d011psi*d100psi*r1r2h)/r3+
  (16*d010psi*d101psi*r1r2h)/r3+(8*d001psi*d110psi*r1r2h)/r3+
  (8*d002psi*d100psi*r1r3h)/r3+(8*d001psi*d101psi*r1r3h)/r3-
  (8*d100psi*d101psi*r1r3h2)/r3-(4*d001psi*d200psi*r1r3h2)/r3+
  (4*d001psi*d100psi*r1r3h2)/(r1*r3)+(4*d001psi*d010psi)/(r2*r3)+
  (8*d002psi*d010psi*r2r3h)/r3+(8*d001psi*d011psi*r2r3h)/r3-
  (8*d011psi*d100psi*r1r3h*r2r3h)/r3-(8*d010psi*d101psi*r1r3h*r2r3h)/r3-
  (8*d001psi*d110psi*r1r3h*r2r3h)/r3-(8*d010psi*d011psi*r2r3h2)/r3-
  (4*d001psi*d020psi*r2r3h2)/r3+(4*d001psi*d010psi*r2r3h2)/(r2*r3);
 dpsiddpsi2=2*d001psi2*d002psi+4*d001psi*d010psi*d011psi+2*d010psi2*d020psi+
  4*d001psi*d100psi*d101psi+4*d010psi*d100psi*d110psi+
  2*d100psi2*d200psi+(2*d001psi2*d100psi)/r1+(2*d010psi2*d100psi)/r1+
  4*d001psi*d011psi*d100psi*r1r2h+4*d010psi*d020psi*d100psi*r1r2h+
  4*d001psi*d010psi*d101psi*r1r2h+4*d010psi2*d110psi*r1r2h+
  4*d100psi2*d110psi*r1r2h+4*d010psi*d100psi*d200psi*r1r2h+
  2*d020psi*d100psi2*r1r2h2+4*d010psi*d100psi*d110psi*r1r2h2+
  2*d010psi2*d200psi*r1r2h2-(2*d010psi2*d100psi*r1r2h2)/r1+
  4*d001psi*d002psi*d100psi*r1r3h+4*d010psi*d011psi*d100psi*r1r3h+
  4*d001psi2*d101psi*r1r3h+4*d100psi2*d101psi*r1r3h+
  4*d001psi*d010psi*d110psi*r1r3h+4*d001psi*d100psi*d200psi*r1r3h+
  4*d011psi*d100psi2*r1r2h*r1r3h+4*d010psi*d100psi*d101psi*r1r2h*r1r3h+
  4*d001psi*d100psi*d110psi*r1r2h*r1r3h+
  4*d001psi*d010psi*d200psi*r1r2h*r1r3h-
  (4*d001psi*d010psi*d100psi*r1r2h*r1r3h)/r1+2*d002psi*d100psi2*r1r3h2+
  4*d001psi*d100psi*d101psi*r1r3h2+2*d001psi2*d200psi*r1r3h2-
  (2*d001psi2*d100psi*r1r3h2)/r1+(2*d001psi2*d010psi)/r2+
  (2*d010psi*d100psi2)/r2-(2*d010psi*d100psi2*r1r2h2)/r2+
  (4*d001psi*d010psi*d100psi*r1r3h)/r2+4*d001psi*d002psi*d010psi*r2r3h+
  4*d001psi2*d011psi*r2r3h+4*d010psi2*d011psi*r2r3h+
  4*d001psi*d010psi*d020psi*r2r3h+4*d010psi*d100psi*d101psi*r2r3h+
  4*d001psi*d100psi*d110psi*r2r3h+(4*d001psi*d010psi*d100psi*r2r3h)/r1+
  4*d010psi*d011psi*d100psi*r1r2h*r2r3h+
  4*d001psi*d020psi*d100psi*r1r2h*r2r3h+4*d010psi2*d101psi*r1r2h*r2r3h+
  4*d001psi*d010psi*d110psi*r1r2h*r2r3h+
  4*d002psi*d010psi*d100psi*r1r3h*r2r3h+
  4*d001psi*d011psi*d100psi*r1r3h*r2r3h+
  4*d001psi*d010psi*d101psi*r1r3h*r2r3h+4*d001psi2*d110psi*r1r3h*r2r3h-
  (4*d001psi*d010psi*d100psi*r1r2h*r2r3h)/r2+2*d002psi*d010psi2*r2r3h2+
  4*d001psi*d010psi*d011psi*r2r3h2+2*d001psi2*d020psi*r2r3h2-
  (2*d001psi2*d010psi*r2r3h2)/r2+(2*d001psi*d010psi2)/r3+
  (2*d001psi*d100psi2)/r3+(4*d001psi*d010psi*d100psi*r1r2h)/r3-
  (2*d001psi*d100psi2*r1r3h2)/r3-
  (4*d001psi*d010psi*d100psi*r1r3h*r2r3h)/r3-
  (2*d001psi*d010psi2*r2r3h2)/r3;

 diidAlnpsi  =2*(dApsi-3*dpsi2/psi)*dipsi2/psi+(2*dpsi2/psi-dApsi)*Lpsi;
 diidAlnpsi /=psi*psi;
 diidAlnpsi +=(LLpsi-(Ldpsi2p2dpsidLpsi-4*dpsiddpsi2/psi)/psi)/psi;

 return(d00dAlnpsi+diidAlnpsi); 
}