lesliegower.m

计算动力学系统的分岔图

tens49=[[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]];
tens410=[[-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]];
tens411=[[24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjy^2*cjj^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,24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjy^3*cjj,24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cjj*cjy^2];[0,0,0,0];[24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyy^2*cyj^2,24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyy^2*cyj*cya,24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyy^3*cyj,-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyj*cyy^2+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cyj*cyy^2];[0,0,0,0]];
tens412=[[24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cjj^2*cjy,-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,24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cjj*cjy^2,24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz^2*cjj*cjy];[0,0,0,0];[-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,-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,-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyj*cyy^2+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cyj*cyy^2,-12*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyj*cyy*cyz+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz^2*cyj*cyy];[0,0,0,0]];
tens413=[[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]];
tens414=[[-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]];
tens415=[[24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cjj^2*cjy,-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,24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cjj*cjy^2,24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz^2*cjj*cjy];[0,0,0,0];[-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,-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,-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyj*cyy^2+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cyj*cyy^2,-12*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyj*cyy*cyz+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz^2*cyj*cyy];[0,0,0,0]];
tens416=[[24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz^2*cjj^2,-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,24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz^2*cjj*cjy,24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz^3*cjj];[0,0,0,0];[-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,-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,-12*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyj*cyy*cyz+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz^2*cyj*cyy,-18*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cyj*cyz^2+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz^3*cyj];[0,0,0,0]];
tens417=[[-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]];
tens418=[[-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]];
tens419=[[-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]];
tens420=[[-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]];
tens421=[[-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]];
tens422=[[-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,-24*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cja^3+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cja^4,-18*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy*cja^2+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjy*cja^3,-18*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cja^2+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cja^3];[0,0,0,0];[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*cya^4,24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyy*cya^3,-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya^3+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cya^3];[0,0,0,0]];
tens423=[[-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,-18*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy*cja^2+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjy*cja^3,-12*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy^2*cja+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjy^2*cja^2,-12*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cja*cjy+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cja^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*cya^2,24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyy*cya^3,24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyy^2*cya^2,-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya^2*cyy+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cya^2*cyy];[0,0,0,0]];
tens424=[[-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,-18*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cja^2+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cja^3,-12*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cja*cjy+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cja^2*cjy,-12*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz^2*cja+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz^2*cja^2];[0,0,0,0];[-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*cya^3+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cya^3,-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya^2*cyy+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cya^2*cyy,-12*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya^2*cyz+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz^2*cya^2];[0,0,0,0]];
tens425=[[-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjj^2*cjy+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cja*cjj^2*cjy,-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*cja*cjj,-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cjy*cjj+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cja*cjy*cjj];[0,0,0,0];[24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cya*cyj^2*cyy,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*cya*cyj,-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya*cyy*cyj+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cya*cyy*cyj];[0,0,0,0]];
tens426=[[-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,-18*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy*cja^2+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjy*cja^3,-12*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy^2*cja+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjy^2*cja^2,-12*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cja*cjy+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cja^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*cya^2,24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyy*cya^3,24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyy^2*cya^2,-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya^2*cyy+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cya^2*cyy];[0,0,0,0]];
tens427=[[-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*cja*cjj,-12*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy^2*cja+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjy^2*cja^2,-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjy^3+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjy^3*cja,-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cjy^2+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cja*cjy^2];[0,0,0,0];[24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyy^2*cya*cyj,24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyy^2*cya^2,24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyy^3*cya,-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya*cyy^2+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cya*cyy^2];[0,0,0,0]];
tens428=[[-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cjy*cjj+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cja*cjy*cjj,-12*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cja*cjy+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cja^2*cjy,-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cjy^2+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cja*cjy^2,-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz^2*cjy+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz^2*cja*cjy];[0,0,0,0];[-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya*cyy*cyj+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cya*cyy*cyj,-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya^2*cyy+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cya^2*cyy,-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya*cyy^2+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cya*cyy^2,-12*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya*cyy*cyz+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz^2*cya*cyy];[0,0,0,0]];
tens429=[[-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjj^2*cjz+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cja*cjj^2*cjz,-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*cjy*cjj+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cja*cjy*cjj,-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*cja*cjj];[0,0,0,0];[-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya*cyj^2+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cya*cyj^2*cyz,-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*cya*cyy*cyj+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cya*cyy*cyj,-12*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya*cyz*cyj+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz^2*cya*cyj];[0,0,0,0]];
tens430=[[-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,-18*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cja^2+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cja^3,-12*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cja*cjy+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cja^2*cjy,-12*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz^2*cja+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz^2*cja^2];[0,0,0,0];[-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*cya^3+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cya^3,-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya^2*cyy+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cya^2*cyy,-12*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya^2*cyz+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz^2*cya^2];[0,0,0,0]];
tens431=[[-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cjy*cjj+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cja*cjy*cjj,-12*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cja*cjy+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cja^2*cjy,-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz*cjy^2+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz*cja*cjy^2,-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz^2*cjy+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz^2*cja*cjy];[0,0,0,0];[-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya*cyy*cyj+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cya*cyy*cyj,-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya^2*cyy+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cya^2*cyy,-6*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya*cyy^2+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz*cya*cyy^2,-12*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya*cyy*cyz+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz^2*cya*cyy];[0,0,0,0]];
tens432=[[-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*cja*cjj,-12*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz^2*cja+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz^2*cja^2,-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz^2*cjy+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz^2*cja*cjy,-6*b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^4*cjz^3+24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^5*cjz^3*cja];[0,0,0,0];[-12*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya*cyz*cyj+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz^2*cya*cyj,-12*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya^2*cyz+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz^2*cya^2,-12*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya*cyy*cyz+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz^2*cya*cyy,-18*b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^4*cya*cyz^2+24*b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^5*cyz^3*cya];[0,0,0,0]];
tens433=[[24*b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*k