psins.r

r软件 另一款可以计算核估计的软件包 需安装r软件

-3*a11*(4*a13*a15 + a11*a35)
if ((r[1]==4) & (r[2]==0) & (r[3]==1) & (r[4]==0) & (r[5]==0) & (r[6]==1))
psi <-
-3*a11*(4*a13*a16 + a11*a36)
if ((r[1]==4) & (r[2]==0) & (r[3]==0) & (r[4]==2) & (r[5]==0) & (r[6]==0))
psi <-
-3*a11*(4*a14^2 + a11*a44)
if ((r[1]==4) & (r[2]==0) & (r[3]==0) & (r[4]==1) & (r[5]==1) & (r[6]==0))
psi <-
-3*a11*(4*a14*a15 + a11*a45)
if ((r[1]==4) & (r[2]==0) & (r[3]==0) & (r[4]==1) & (r[5]==0) & (r[6]==1))
psi <-
-3*a11*(4*a14*a16 + a11*a46)
if ((r[1]==4) & (r[2]==0) & (r[3]==0) & (r[4]==0) & (r[5]==2) & (r[6]==0))
psi <-
-3*a11*(4*a15^2 + a11*a55)
  if ((r[1]==4) & (r[2]==0) & (r[3]==0) & (r[4]==0) & (r[5]==1) & (r[6]==1))
psi <-
-3*a11*(4*a15*a16 + a11*a56)
if ((r[1]==4) & (r[2]==0) & (r[3]==0) & (r[4]==0) & (r[5]==0) & (r[6]==2))
psi <-
-3*a11*(4*a16^2 + a11*a66)
if ((r[1]==3) & (r[2]==3) & (r[3]==0) & (r[4]==0) & (r[5]==0) & (r[6]==0))
psi <-
-6*a12^3 - 9*a11*a12*a22
if ((r[1]==3) & (r[2]==2) & (r[3]==1) & (r[4]==0) & (r[5]==0) & (r[6]==0))
psi <-
-3*(2*a12^2*a13 + a11*a13*a22 + 2*a11*a12*a23)
if ((r[1]==3) & (r[2]==2) & (r[3]==0) & (r[4]==1) & (r[5]==0) & (r[6]==0))
psi <-
-3*(2*a12^2*a14 + a11*a14*a22 + 2*a11*a12*a24)
if ((r[1]==3) & (r[2]==2) & (r[3]==0) & (r[4]==0) & (r[5]==1) & (r[6]==0))
psi <-
-3*(2*a12^2*a15 + a11*a15*a22 + 2*a11*a12*a25)
if ((r[1]==3) & (r[2]==2) & (r[3]==0) & (r[4]==0) & (r[5]==0) & (r[6]==1))
psi <-
-3*(2*a12^2*a16 + a11*a16*a22 + 2*a11*a12*a26)
if ((r[1]==3) & (r[2]==1) & (r[3]==2) & (r[4]==0) & (r[5]==0) & (r[6]==0))
psi <-
-3*(2*a12*a13^2 + 2*a11*a13*a23 + a11*a12*a33)
if ((r[1]==3) & (r[2]==1) & (r[3]==1) & (r[4]==1) & (r[5]==0) & (r[6]==0))
psi <-
-3*(2*a12*a13*a14 + a11*a14*a23 + a11*a13*a24 + a11*a12*a34)
if ((r[1]==3) & (r[2]==1) & (r[3]==1) & (r[4]==0) & (r[5]==1) & (r[6]==0))
psi <-
-3*(2*a12*a13*a15 + a11*a15*a23 + a11*a13*a25 + a11*a12*a35)
if ((r[1]==3) & (r[2]==1) & (r[3]==1) & (r[4]==0) & (r[5]==0) & (r[6]==1))
psi <-
-3*(2*a12*a13*a16 + a11*a16*a23 + a11*a13*a26 + a11*a12*a36)
if ((r[1]==3) & (r[2]==1) & (r[3]==0) & (r[4]==2) & (r[5]==0) & (r[6]==0))
psi <-
-3*(2*a12*a14^2 + 2*a11*a14*a24 + a11*a12*a44)
if ((r[1]==3) & (r[2]==1) & (r[3]==0) & (r[4]==1) & (r[5]==1) & (r[6]==0))
psi <-
-3*(2*a12*a14*a15 + a11*a15*a24 + a11*a14*a25 + a11*a12*a45)
if ((r[1]==3) & (r[2]==1) & (r[3]==0) & (r[4]==1) & (r[5]==0) & (r[6]==1))
psi <-
-3*(2*a12*a14*a16 + a11*a16*a24 + a11*a14*a26 + a11*a12*a46)
if ((r[1]==3) & (r[2]==1) & (r[3]==0) & (r[4]==0) & (r[5]==2) & (r[6]==0))
psi <-
-3*(2*a12*a15^2 + 2*a11*a15*a25 + a11*a12*a55)
if ((r[1]==3) & (r[2]==1) & (r[3]==0) & (r[4]==0) & (r[5]==1) & (r[6]==1))
psi <-
-3*(2*a12*a15*a16 + a11*a16*a25 + a11*a15*a26 + a11*a12*a56)
if ((r[1]==3) & (r[2]==1) & (r[3]==0) & (r[4]==0) & (r[5]==0) & (r[6]==2))
psi <-
-3*(2*a12*a16^2 + 2*a11*a16*a26 + a11*a12*a66)
if ((r[1]==3) & (r[2]==0) & (r[3]==3) & (r[4]==0) & (r[5]==0) & (r[6]==0))
psi <-
-6*a13^3 - 9*a11*a13*a33
if ((r[1]==3) & (r[2]==0) & (r[3]==2) & (r[4]==1) & (r[5]==0) & (r[6]==0))
psi <-
-3*(2*a13^2*a14 + a11*a14*a33 + 2*a11*a13*a34)
  if ((r[1]==3) & (r[2]==0) & (r[3]==2) & (r[4]==0) & (r[5]==1) & (r[6]==0))
psi <-
-3*(2*a13^2*a15 + a11*a15*a33 + 2*a11*a13*a35)
if ((r[1]==3) & (r[2]==0) & (r[3]==2) & (r[4]==0) & (r[5]==0) & (r[6]==1))
psi <-
-3*(2*a13^2*a16 + a11*a16*a33 + 2*a11*a13*a36)
if ((r[1]==3) & (r[2]==0) & (r[3]==1) & (r[4]==2) & (r[5]==0) & (r[6]==0))
psi <-
-3*(2*a13*a14^2 + 2*a11*a14*a34 + a11*a13*a44)
if ((r[1]==3) & (r[2]==0) & (r[3]==1) & (r[4]==1) & (r[5]==1) & (r[6]==0))
psi <-
-3*(2*a13*a14*a15 + a11*a15*a34 + a11*a14*a35 + a11*a13*a45)
if ((r[1]==3) & (r[2]==0) & (r[3]==1) & (r[4]==1) & (r[5]==0) & (r[6]==1))
psi <-
-3*(2*a13*a14*a16 + a11*a16*a34 + a11*a14*a36 + a11*a13*a46)
if ((r[1]==3) & (r[2]==0) & (r[3]==1) & (r[4]==0) & (r[5]==2) & (r[6]==0))
psi <-
-3*(2*a13*a15^2 + 2*a11*a15*a35 + a11*a13*a55)
if ((r[1]==3) & (r[2]==0) & (r[3]==1) & (r[4]==0) & (r[5]==1) & (r[6]==1))
psi <-
-3*(2*a13*a15*a16 + a11*a16*a35 + a11*a15*a36 + a11*a13*a56)
if ((r[1]==3) & (r[2]==0) & (r[3]==1) & (r[4]==0) & (r[5]==0) & (r[6]==2))
psi <-
-3*(2*a13*a16^2 + 2*a11*a16*a36 + a11*a13*a66)
if ((r[1]==3) & (r[2]==0) & (r[3]==0) & (r[4]==3) & (r[5]==0) & (r[6]==0))
psi <-
-6*a14^3 - 9*a11*a14*a44
if ((r[1]==3) & (r[2]==0) & (r[3]==0) & (r[4]==2) & (r[5]==1) & (r[6]==0))
psi <-
-3*(2*a14^2*a15 + a11*a15*a44 + 2*a11*a14*a45)
if ((r[1]==3) & (r[2]==0) & (r[3]==0) & (r[4]==2) & (r[5]==0) & (r[6]==1))
psi <-
-3*(2*a14^2*a16 + a11*a16*a44 + 2*a11*a14*a46)
if ((r[1]==3) & (r[2]==0) & (r[3]==0) & (r[4]==1) & (r[5]==2) & (r[6]==0))
psi <-
-3*(2*a14*a15^2 + 2*a11*a15*a45 + a11*a14*a55)
if ((r[1]==3) & (r[2]==0) & (r[3]==0) & (r[4]==1) & (r[5]==1) & (r[6]==1))
psi <-
-3*(2*a14*a15*a16 + a11*a16*a45 + a11*a15*a46 + a11*a14*a56)
if ((r[1]==3) & (r[2]==0) & (r[3]==0) & (r[4]==1) & (r[5]==0) & (r[6]==2))
psi <-
-3*(2*a14*a16^2 + 2*a11*a16*a46 + a11*a14*a66)
if ((r[1]==3) & (r[2]==0) & (r[3]==0) & (r[4]==0) & (r[5]==3) & (r[6]==0))
psi <-
-6*a15^3 - 9*a11*a15*a55
if ((r[1]==3) & (r[2]==0) & (r[3]==0) & (r[4]==0) & (r[5]==2) & (r[6]==1))
psi <-
-3*(2*a15^2*a16 + a11*a16*a55 + 2*a11*a15*a56)
if ((r[1]==3) & (r[2]==0) & (r[3]==0) & (r[4]==0) & (r[5]==1) & (r[6]==2))
psi <-
-3*(2*a15*a16^2 + 2*a11*a16*a56 + a11*a15*a66)
if ((r[1]==3) & (r[2]==0) & (r[3]==0) & (r[4]==0) & (r[5]==0) & (r[6]==3))
psi <-
-6*a16^3 - 9*a11*a16*a66
if ((r[1]==2) & (r[2]==4) & (r[3]==0) & (r[4]==0) & (r[5]==0) & (r[6]==0))
psi <-
-3*a22*(4*a12^2 + a11*a22)
if ((r[1]==2) & (r[2]==3) & (r[3]==1) & (r[4]==0) & (r[5]==0) & (r[6]==0))
psi <-
-3*(2*a12*a13*a22 + 2*a12^2*a23 + a11*a22*a23)
if ((r[1]==2) & (r[2]==3) & (r[3]==0) & (r[4]==1) & (r[5]==0) & (r[6]==0))
psi <-
-3*(2*a12*a14*a22 + 2*a12^2*a24 + a11*a22*a24)
if ((r[1]==2) & (r[2]==3) & (r[3]==0) & (r[4]==0) & (r[5]==1) & (r[6]==0))
psi <-
-3*(2*a12*a15*a22 + 2*a12^2*a25 + a11*a22*a25)
if ((r[1]==2) & (r[2]==3) & (r[3]==0) & (r[4]==0) & (r[5]==0) & (r[6]==1))
psi <-
-3*(2*a12*a16*a22 + 2*a12^2*a26 + a11*a22*a26)
if ((r[1]==2) & (r[2]==2) & (r[3]==2) & (r[4]==0) & (r[5]==0) & (r[6]==0))
psi <-
-2*a13^2*a22 - 8*a12*a13*a23 - 2*a11*a23^2 - 2*a12^2*a33 - a11*a22*a33
if ((r[1]==2) & (r[2]==2) & (r[3]==1) & (r[4]==1) & (r[5]==0) & (r[6]==0))
psi <-
-2*a13*a14*a22 - 4*a12*a14*a23 - 4*a12*a13*a24 - 2*a11*a23*a24 - 
 2*a12^2*a34 - a11*a22*a34
if ((r[1]==2) & (r[2]==2) & (r[3]==1) & (r[4]==0) & (r[5]==1) & (r[6]==0))
psi <-
-2*a13*a15*a22 - 4*a12*a15*a23 - 4*a12*a13*a25 - 2*a11*a23*a25 - 
 2*a12^2*a35 - a11*a22*a35
if ((r[1]==2) & (r[2]==2) & (r[3]==1) & (r[4]==0) & (r[5]==0) & (r[6]==1))
psi <-
-2*a13*a16*a22 - 4*a12*a16*a23 - 4*a12*a13*a26 - 2*a11*a23*a26 - 
 2*a12^2*a36 - a11*a22*a36
if ((r[1]==2) & (r[2]==2) & (r[3]==0) & (r[4]==2) & (r[5]==0) & (r[6]==0))
psi <-
-2*a14^2*a22 - 8*a12*a14*a24 - 2*a11*a24^2 - 2*a12^2*a44 - a11*a22*a44
if ((r[1]==2) & (r[2]==2) & (r[3]==0) & (r[4]==1) & (r[5]==1) & (r[6]==0))
psi <-
-2*a14*a15*a22 - 4*a12*a15*a24 - 4*a12*a14*a25 - 2*a11*a24*a25 - 
 2*a12^2*a45 - a11*a22*a45
if ((r[1]==2) & (r[2]==2) & (r[3]==0) & (r[4]==1) & (r[5]==0) & (r[6]==1))
psi <-
-2*a14*a16*a22 - 4*a12*a16*a24 - 4*a12*a14*a26 - 2*a11*a24*a26 - 
 2*a12^2*a46 - a11*a22*a46
if ((r[1]==2) & (r[2]==2) & (r[3]==0) & (r[4]==0) & (r[5]==2) & (r[6]==0))
psi <-
-2*a15^2*a22 - 8*a12*a15*a25 - 2*a11*a25^2 - 2*a12^2*a55 - a11*a22*a55
if ((r[1]==2) & (r[2]==2) & (r[3]==0) & (r[4]==0) & (r[5]==1) & (r[6]==1))
psi <-
-2*a15*a16*a22 - 4*a12*a16*a25 - 4*a12*a15*a26 - 2*a11*a25*a26 - 
 2*a12^2*a56 - a11*a22*a56
if ((r[1]==2) & (r[2]==2) & (r[3]==0) & (r[4]==0) & (r[5]==0) & (r[6]==2))
psi <-
-2*a16^2*a22 - 8*a12*a16*a26 - 2*a11*a26^2 - 2*a12^2*a66 - a11*a22*a66
if ((r[1]==2) & (r[2]==1) & (r[3]==3) & (r[4]==0) & (r[5]==0) & (r[6]==0))
psi <-
-3*(2*a13^2*a23 + 2*a12*a13*a33 + a11*a23*a33)
  if ((r[1]==2) & (r[2]==1) & (r[3]==2) & (r[4]==1) & (r[5]==0) & (r[6]==0))
psi <-
-4*a13*a14*a23 - 2*a13^2*a24 - 2*a12*a14*a33 - a11*a24*a33 - 4*a12*a13*a34 - 
 2*a11*a23*a34
if ((r[1]==2) & (r[2]==1) & (r[3]==2) & (r[4]==0) & (r[5]==1) & (r[6]==0))
psi <-
-4*a13*a15*a23 - 2*a13^2*a25 - 2*a12*a15*a33 - a11*a25*a33 - 4*a12*a13*a35 - 
 2*a11*a23*a35
if ((r[1]==2) & (r[2]==1) & (r[3]==2) & (r[4]==0) & (r[5]==0) & (r[6]==1))
psi <-
-4*a13*a16*a23 - 2*a13^2*a26 - 2*a12*a16*a33 - a11*a26*a33 - 4*a12*a13*a36 - 
 2*a11*a23*a36
if ((r[1]==2) & (r[2]==1) & (r[3]==1) & (r[4]==2) & (r[5]==0) & (r[6]==0))
psi <-
-2*a14^2*a23 - 4*a13*a14*a24 - 4*a12*a14*a34 - 2*a11*a24*a34 - 
 2*a12*a13*a44 - a11*a23*a44
if ((r[1]==2) & (r[2]==1) & (r[3]==1) & (r[4]==1) & (r[5]==1) & (r[6]==0))
psi <-
-2*a14*a15*a23 - 2*a13*a15*a24 - 2*a13*a14*a25 - 2*a12*a15*a34 - 
 a11*a25*a34 - 2*a12*a14*a35 - a11*a24*a35 - 2*a12*a13*a45 - a11*a23*a45
if ((r[1]==2) & (r[2]==1) & (r[3]==1) & (r[4]==1) 
& (r[5]==0) & (r[6]==1))
psi <-
-2*a14*a16*a23 - 2*a13*a16*a24 - 2*a13*a14*a26 - 2*a12*a16*a34 - 
 a11*a26*a34 - 2*a12*a14*a36 - a11*a24*a36 - 2*a12*a13*a46 - a11*a23*a46
if ((r[1]==2) & (r[2]==1) & (r[3]==1) & (r[4]==0) 
& (r[5]==2) & (r[6]==0))
psi <-
-2*a15^2*a23 - 4*a13*a15*a25 - 4*a12*a15*a35 - 2*a11*a25*a35 - 
 2*a12*a13*a55 - a11*a23*a55
if ((r[1]==2) & (r[2]==1) & (r[3]==1) & (r[4]==0) 
& (r[5]==1) & (r[6]==1))
psi <-
-2*a15*a16*a23 - 2*a13*a16*a25 - 2*a13*a15*a26 - 2*a12*a16*a35 - 
 a11*a26*a35 - 2*a12*a15*a36 - a11*a25*a36 - 2*a12*a13*a56 - a11*a23*a56
if ((r[1]==2) & (r[2]==1) & (r[3]==1) & (r[4]==0) 
& (r[5]==0) & (r[6]==2))
psi <-
-2*a16^2*a23 - 4*a13*a16*a26 - 4*a12*a16*a36 - 2*a11*a26*a36 - 
 2*a12*a13*a66 - a11*a23*a66
if ((r[1]==2) & (r[2]==1) & (r[3]==0) & (r[4]==3) 
& (r[5]==0) & (r[6]==0))
psi <-
-3*(2*a14^2*a24 + 2*a12*a14*a44 + a11*a24*a44)
if ((r[1]==2) & (r[2]==1) & (r[3]==0) & (r[4]==2) 
& (r[5]==1) & (r[6]==0))
psi <-
-4*a14*a15*a24 - 2*a14^2*a25 - 2*a12*a15*a44 - a11*a25*a44 - 4*a12*a14*a45 - 
 2*a11*a24*a45
if ((r[1]==2) & (r[2]==1) & (r[3]==0) & (r[4]==2) 
& (r[5]==0) & (r[6]==1))
psi <-
-4*a14*a16*a24 - 2*a14^2*a26 - 2*a12*a16*a44 - a11*a26*a44 - 4*a12*a14*a46 - 
 2*a11*a24*a46
if ((r[1]==2) & (r[2]==1) & (r[3]==0) & (r[4]==1) 
& (r[5]==2) & (r[6]==0))
psi <-
-2*a15^2*a24 - 4*a14*a15*a25 - 4*a12*a15*a45 - 2*a11*a25*a45 - 
 2*a12*a14*a55 - a11*a24*a55
if ((r[1]==2) & (r[2]==1) & (r[3]==0) & (r[4]==1) 
& (r[5]==1) & (r[6]==1))
psi <-
-2*a15*a16*a24 - 2*a14*a16*a25 - 2*a14*a15*a26 - 2*a12*a16*a45 - 
 a11*a26*a45 - 2*a12*a15*a46 - a11*a25*a46 - 2*a12*a14*a56 - a11*a24*a56
if ((r[1]==2) & (r[2]==1) & (r[3]==0) & (r[4]==1) 
& (r[5]==0) & (r[6]==2))
psi <-
-2*a16^2*a24 - 4*a14*a16*a26 - 4*a12*a16*a46 - 2*a11*a26*a46 - 
 2*a12*a14*a66 - a11*a24*a66
if ((r[1]==2) & (r[2]==1) & (r[3]==0) & (r[4]==0) 
& (r[5]==3) & (r[6]==0))
psi <-
-3*(2*a15^2*a25 + 2*a12*a15*a55 + a11*a25*a55)
if ((r[1]==2) & (r[2]==1) & (r[3]==0) & (r[4]==0) 
& (r[5]==2) & (r[6]==1))
psi <-
-4*a15*a16*a25 - 2*a15^2*a26 - 2*a12*a16*a55 - a11*a26*a55 - 4*a12*a15*a56 - 
 2*a11*a25*a56
if ((r[1]==2) & (r[2]==1) & (r[3]==0) & (r[4]==0) 
& (r[5]==1) & (r[6]==2))
psi <-
-2*a16^2*a25 - 4*a15*a16*a26 - 4*a12*a16*a56 - 2*a11*a26*a56 - 
 2*a12*a15*a66 - a11*a25*a66
if ((r[1]==2) & (r[2]==1) & (r[3]==0) & (r[4]==0) 
& (r[5]==0) & (r[6]==3))
psi <-
-3*(2*a16^2*a26 + 2*a12*a16*a66 + a11*a26*a66)
if ((r[1]==2) & (r[2]==0) & (r[3]==4) & (r[4]==0) 
& (r[5]==0) & (r[6]==0))
psi <-
-3*a33*(4*a13^2 + a11*a33)
if ((r[1]==2) & (r[2]==0) & (r[3]==3) & (r[4]==1) 
& (r[5]==0) & (r[6]==0))
psi <-
-3*(2*a13*a14*a33 + 2*a13^2*a34 + a11*a33*a34)
if ((r[1]==2) & (r[2]==0) & (r[3]==3) & (r[4]==0) 
& (r[5]==1) & (r[6]==0))
psi <-
-3*(2*a13*a15*a33 + 2*a13^2*a35 + a11*a33*a35)
if ((r[1]==2) & (r[2]==0) & (r[3]==3) & (r[4]==0) 
& (r[5]==0) & (r[6]==1))
psi <-
-3*(2*a13*a16*a33 + 2*a13^2*a36 + a11*a33*a36)
if ((r[1]==2) & (r[2]==0) & (r[3]==2) & (r[4]==2) 
& (r[5]==0) & (r[6]==0))
psi <-
-2*a14^2*a33 - 8*a13*a14*a34 - 2*a11*a34^2 - 2*a13^2*a44 - a11*a33*a44
if ((r[1]==2) & (r[2]==0) & (r[3]==2) & (r[4]==1) 
& (r[5]==1) & (r[6]==0))
psi <-
-2*a14*a15*a33 - 4*a13*a15*a34 - 4*a13*a14*a35 - 2*a11*a34*a35 - 
 2*a13^2*a45 - a11*a33*a45
  if ((r[1]==2) & (r[2]==0) & (r[3]==2) & (r[4]==1) 
& (r[5]==0) & (r[6]==1))
psi <-
-2*a14*a16*a33 - 4*a13*a16*a34 - 4*a13*a14*a36 - 2*a11*a34*a36 - 
 2*a13^2*a46 - a11*a33*a46
if ((r[1]==2) & (r[2]==0) & (r[3]==2) & (r[4]==0) 
& (r[5]==2) & (r[6]==0))
psi <-
-2*a15^2*a33 - 8*a13*a15*a35 - 2*a11*a35^2 - 2*a13^2*a55 - a11*a33*a55
if ((r[1]==2) & (r[2]==0) & (r[3]==2) & (r[4]==0) 
& (r[5]==1) & (r[6]==1))
psi <-
-2*a15*a16*a33 - 4*a13*a16*a35 - 4*a13*a15*a36 - 2*a11*a35*a36 - 
 2*a13^2*a56 - a11*a33*a56
if ((r[1]==2) & (r[2]==0) & (r[3]==2) & (r[4]==0) 
& (r[5]==0) & (r[6]==2))
psi <-
-2*a16^2*a33 - 8*a13*a16*a36 - 2*a11*a36^2 - 2*a13^2*a66 - a11*a33*a66
if ((r[1]==2) & (r[2]==0) & (r[3]==1) & (r[4]==3) 
& (r[5]==0) & (r[6]==0))
psi <-
-3*(2*a14^2*a34 + 2*a13*a14*a44 + a11*a34*a44)
if ((r[1]==2) & (r[2]==0) & (r[3]==1) & (r[4]==2) 
& (r[5]==1) & (r[6]==0))
psi <-
-4*a14*a15*a34 - 2*a14^2*a35 - 2*a13*a15*a44 - a11*a35*a44 - 4*a13*a14*a45 - 
 2*a11*a34*a45
if ((r[1]==2) & (r[2]==0) & (r[3]==1) & (r[4]==2) 
& (r[5]==0) & (r[6]==1))
psi <-
-4*a14*a16*a34 - 2*a14^2*a36 - 2*a13*a16*a44 - a11*a36*a44 - 4*a13*a14*a46 - 
 2*a11*a34*a46
if ((r[1]==2) & (r[2]==0) & (r[3]==1) & (r[4]==1) 
& (r[5]==2) & (r[6]==0))
psi <-
-2*a15^2*a34 - 4*a14*a15*a35 - 4*a13*a15*a45 - 2*a11*a35*a45 - 
 2*a13*a14*a55 - a11*a34*a55
if ((r[1]==2) & (r[2]==0) & (r[3]==1) & (r[4]==1) 
& (r[5]==1) & (r[6]==1))
psi <-
-2*a15*a16*a34 - 2*a14*a16*a35 - 2*a14*a15*a36 - 2*a13*a16*a45 - 
 a11*a36*a45 - 2*a13*a15*a46 - a11*a35*a46 - 2*a13*a14*a56 - a11*a34*a56
if ((r[1]==2) & (r[2]==0) & (r[3]==1) & (r[4]==1) 
& (r[5]==0) & (r[6]==2))
psi <-
-2*a16^2*a34 - 4*a14*a16*a36 - 4*a13*a16*a46 - 2*a11*a36*a46 - 
 2*a13*a14*a66 - a11*a34*a66
if ((r[1]==2) & (r[2]==0) & (r[3]==1) & (r[4]==0) 
& (r[5]==3) & (r[6]==0))