tens314=[[2*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^3*cjz*cjj-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cjj*cja,4*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^3*cjz*cja-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cja^2,2*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^3*cjz*cjy-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cja*cjy,2*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^3*cjz^2-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz^2*cja];[0,0,0,0];[2*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^3*cyj*cya-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyz*cyj*cya,2*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^3*cya^2-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyz*cya^2,2*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^3*cya*cyy-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyz*cya*cyy,4*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^3*cya*cyz-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyz^2*cya];[0,0,0,0]];tens315=[[-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cjj*cjy,2*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^3*cjz*cjy-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cja*cjy,-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cjy^2,-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz^2*cjy];[0,0,0,0];[2*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^3*cyj*cyy-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyz*cyj*cyy,2*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^3*cya*cyy-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyz*cya*cyy,2*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^3*cyy^2-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyz*cyy^2,4*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^3*cyy*cyz-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyz^2*cyy];[0,0,0,0]];tens316=[[-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz^2*cjj,2*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^3*cjz^2-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz^2*cja,-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz^2*cjy,-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz^3];[0,0,0,0];[4*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^3*cyj*cyz-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyz^2*cyj,4*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^3*cya*cyz-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyz^2*cya,4*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^3*cyy*cyz-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyz^2*cyy,6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^3*cyz^2-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyz^3];[0,0,0,0]];tens3(:,:,1,1) =tens31;tens3(:,:,1,2) =tens32;tens3(:,:,1,3) =tens33;tens3(:,:,1,4) =tens34;tens3(:,:,2,1) =tens35;tens3(:,:,2,2) =tens36;tens3(:,:,2,3) =tens37;tens3(:,:,2,4) =tens38;tens3(:,:,3,1) =tens39;tens3(:,:,3,2) =tens310;tens3(:,:,3,3) =tens311;tens3(:,:,3,4) =tens312;tens3(:,:,4,1) =tens313;tens3(:,:,4,2) =tens314;tens3(:,:,4,3) =tens315;tens3(:,:,4,4) =tens316;%---------------------------------------------------------------------------function tens4  = der4(t,kmrgd,b1,b2,ma,mj,my,mz,cjj,cja,cyy,cyz,cjy,cjz,cyj,cya)tens41=[[24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjj^4,-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjj^3+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cja*cjj^3,24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjy*cjj^3,24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cjj^3];[0,0,0,0];[24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyj^4,24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cya*cyj^3,24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyy*cyj^3,-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyj^3+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cyj^3];[0,0,0,0]];tens42=[[-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjj^3+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cja*cjj^3,-12*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjj^2*cja+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cja^2*cjj^2,-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy*cjj^2+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjy*cjj^2*cja,-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cjj^2+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cjj^2*cja];[0,0,0,0];[24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cya*cyj^3,24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cya^2*cyj^2,24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyy*cyj^2*cya,-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyj^2*cya+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cyj^2*cya];[0,0,0,0]];tens43=[[24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjy*cjj^3,-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy*cjj^2+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjy*cjj^2*cja,24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjy^2*cjj^2,24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cjj^2*cjy];[0,0,0,0];[24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyy*cyj^3,24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyy*cyj^2*cya,24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyy^2*cyj^2,-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyj^2*cyy+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cyj^2*cyy];[0,0,0,0]];tens44=[[24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cjj^3,-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cjj^2+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cjj^2*cja,24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cjj^2*cjy,24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz^2*cjj^2];[0,0,0,0];[-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyj^3+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cyj^3,-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyj^2*cya+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cyj^2*cya,-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyj^2*cyy+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cyj^2*cyy,-12*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyj^2*cyz+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz^2*cyj^2];[0,0,0,0]];tens45=[[-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjj^3+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cja*cjj^3,-12*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjj^2*cja+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cja^2*cjj^2,-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy*cjj^2+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjy*cjj^2*cja,-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cjj^2+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cjj^2*cja];[0,0,0,0];[24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cya*cyj^3,24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cya^2*cyj^2,24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyy*cyj^2*cya,-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyj^2*cya+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cyj^2*cya];[0,0,0,0]];tens46=[[-12*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjj^2*cja+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cja^2*cjj^2,-18*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjj*cja^2+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cja^3*cjj,-12*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy*cjj*cja+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjy*cjj*cja^2,-12*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cjj*cja+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cjj*cja^2];[0,0,0,0];[24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cya^2*cyj^2,24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cya^3*cyj,24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyy*cyj*cya^2,-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyj*cya^2+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cyj*cya^2];[0,0,0,0]];tens47=[[-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy*cjj^2+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjy*cjj^2*cja,-12*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy*cjj*cja+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjy*cjj*cja^2,-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy^2*cjj+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjy^2*cjj*cja,-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cjj*cjy+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cjj*cja*cjy];[0,0,0,0];[24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyy*cyj^2*cya,24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyy*cyj*cya^2,24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyy^2*cyj*cya,-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyj*cya*cyy+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cyj*cya*cyy];[0,0,0,0]];tens48=[[-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cjj^2+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cjj^2*cja,-12*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cjj*cja+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cjj*cja^2,-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cjj*cjy+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cjj*cja*cjy,-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz^2*cjj+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz^2*cjj*cja];[0,0,0,0];[-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyj^2*cya+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cyj^2*cya,-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyj*cya^2+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cyj*cya^2,-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyj*cya*cyy+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cyj*cya*cyy,-12*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyj*cya*cyz+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz^2*cyj*cya];[0,0,0,0]];