hessp(:,:,1) =hessp1;hessp(:,:,2) =hessp2;hessp(:,:,3) =hessp3;hessp(:,:,4) =hessp4;hessp(:,:,5) =hessp5;hessp(:,:,6) =hessp6;hessp(:,:,7) =hessp7;hessp(:,:,8) =hessp8;hessp(:,:,9) =hessp9;hessp(:,:,10) =hessp10;hessp(:,:,11) =hessp11;hessp(:,:,12) =hessp12;hessp(:,:,13) =hessp13;hessp(:,:,14) =hessp14;%---------------------------------------------------------------------------function tens3  = der3(t,kmrgd,b1,b2,ma,mj,my,mz,cjj,cja,cyy,cyz,cjy,cjz,cyj,cya)tens31=[[-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjj^3,2*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^3*cjj^2-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cja*cjj^2,-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy*cjj^2,-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cjj^2];[0,0,0,0];[-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyj^3,-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya*cyj^2,-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyy*cyj^2,2*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^3*cyj^2-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyz*cyj^2];[0,0,0,0]];tens32=[[2*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^3*cjj^2-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cja*cjj^2,4*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^3*cjj*cja-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cja^2*cjj,2*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^3*cjy*cjj-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy*cjj*cja,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];[0,0,0,0];[-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya*cyj^2,-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya^2*cyj,-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyy*cyj*cya,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];[0,0,0,0]];tens33=[[-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy*cjj^2,2*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^3*cjy*cjj-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy*cjj*cja,-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy^2*cjj,-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cjj*cjy];[0,0,0,0];[-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyy*cyj^2,-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyy*cyj*cya,-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyy^2*cyj,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];[0,0,0,0]];tens34=[[-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cjj^2,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,-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cjj*cjy,-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz^2*cjj];[0,0,0,0];[2*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^3*cyj^2-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyz*cyj^2,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*cyj*cyy-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyz*cyj*cyy,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];[0,0,0,0]];tens35=[[2*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^3*cjj^2-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cja*cjj^2,4*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^3*cjj*cja-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cja^2*cjj,2*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^3*cjy*cjj-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy*cjj*cja,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];[0,0,0,0];[-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya*cyj^2,-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya^2*cyj,-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyy*cyj*cya,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];[0,0,0,0]];tens36=[[4*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^3*cjj*cja-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cja^2*cjj,6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^3*cja^2-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cja^3,4*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^3*cjy*cja-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy*cja^2,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];[0,0,0,0];[-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya^2*cyj,-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya^3,-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyy*cya^2,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];[0,0,0,0]];tens37=[[2*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^3*cjy*cjj-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy*cjj*cja,4*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^3*cjy*cja-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy*cja^2,2*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^3*cjy^2-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy^2*cja,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];[0,0,0,0];[-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyy*cyj*cya,-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyy*cya^2,-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyy^2*cya,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];[0,0,0,0]];tens38=[[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]];tens39=[[-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy*cjj^2,2*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^3*cjy*cjj-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy*cjj*cja,-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy^2*cjj,-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cjj*cjy];[0,0,0,0];[-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyy*cyj^2,-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyy*cyj*cya,-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyy^2*cyj,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];[0,0,0,0]];tens310=[[2*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^3*cjy*cjj-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy*cjj*cja,4*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^3*cjy*cja-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy*cja^2,2*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^3*cjy^2-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy^2*cja,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];[0,0,0,0];[-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyy*cyj*cya,-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyy*cya^2,-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyy^2*cya,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];[0,0,0,0]];tens311=[[-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy^2*cjj,2*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^3*cjy^2-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy^2*cja,-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy^3,-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cjy^2];[0,0,0,0];[-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyy^2*cyj,-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyy^2*cya,-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyy^3,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];[0,0,0,0]];tens312=[[-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]];tens313=[[-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cjj^2,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,-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cjj*cjy,-6*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz^2*cjj];[0,0,0,0];[2*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^3*cyj^2-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyz*cyj^2,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*cyj*cyy-6*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyz*cyj*cyy,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];[0,0,0,0]];