%load AA.mat 
%for i=1:15
%   for j=1:4
%      fprintf('At2(%i,%i) = %s;\n',i,j,char(At2(i,j)));
%   end
%end

At1(1,1) = -u11^2*up31+u11*up11*u31;
At1(1,2) = u11^2*u21;
At1(1,3) = -u11*u31^2-u11^3;
At1(1,4) = u11*u21*u31;
At1(2,1) = -u21*up31*u11+u31*up11*u21;
At1(2,2) = u11*u21^2;
At1(2,3) = -u11^2*u21-u21*u31^2;
At1(2,4) = u31*u21^2;
At1(3,1) = -u31*up31*u11+u31^2*up11;
At1(3,2) = u11*u21*u31;
At1(3,3) = -u31^3-u11^2*u31;
At1(3,4) = u21*u31^2;
At1(4,1) = -u12^2*up32+u12*up12*u32;
At1(4,2) = u12^2*u22;
At1(4,3) = -u12*u32^2-u12^3;
At1(4,4) = u12*u22*u32;
At1(5,1) = -u22*up32*u12+u32*up12*u22;
At1(5,2) = u12*u22^2;
At1(5,3) = -u12^2*u22-u22*u32^2;
At1(5,4) = u32*u22^2;
At1(6,1) = -u32*up32*u12+u32^2*up12;
At1(6,2) = u12*u22*u32;
At1(6,3) = -u32^3-u12^2*u32;
At1(6,4) = u22*u32^2;
At1(7,1) = -u13^2*up33+u13*up13*u33;
At1(7,2) = u13^2*u23;
At1(7,3) = -u13*u33^2-u13^3;
At1(7,4) = u13*u23*u33;
At1(8,1) = -u23*up33*u13+u33*up13*u23;
At1(8,2) = u13*u23^2;
At1(8,3) = -u13^2*u23-u23*u33^2;
At1(8,4) = u33*u23^2;
At1(9,1) = -u33*up33*u13+u33^2*up13;
At1(9,2) = u13*u23*u33;
At1(9,3) = -u33^3-u13^2*u33;
At1(9,4) = u23*u33^2;
At1(10,1) = -u14^2*up34+u14*up14*u34;
At1(10,2) = u14^2*u24;
At1(10,3) = -u14*u34^2-u14^3;
At1(10,4) = u14*u24*u34;
At1(11,1) = -u24*up34*u14+u34*up14*u24;
At1(11,2) = u14*u24^2;
At1(11,3) = -u14^2*u24-u24*u34^2;
At1(11,4) = u34*u24^2;
At1(12,1) = -u34*up34*u14+u34^2*up14;
At1(12,2) = u14*u24*u34;
At1(12,3) = -u34^3-u14^2*u34;
At1(12,4) = u24*u34^2;
At1(13,1) = -u15^2*up35+u15*up15*u35;
At1(13,2) = u15^2*u25;
At1(13,3) = -u15*u35^2-u15^3;
At1(13,4) = u25*u35*u15;
At1(14,1) = -u25*up35*u15+u35*up15*u25;
At1(14,2) = u15*u25^2;
At1(14,3) = -u15^2*u25-u25*u35^2;
At1(14,4) = u35*u25^2;
At1(15,1) = -u35*up35*u15+u35^2*up15;
At1(15,2) = u25*u35*u15;
At1(15,3) = -u35^3-u15^2*u35;
At1(15,4) = u25*u35^2;


At2(1,1) = -u11*up11*u21+u11^2*up21;
At2(1,2) = u11^2*u31;
At2(1,3) = u11*u21*u31;
At2(1,4) = -u11^3-u11*u21^2;
At2(2,1) = -u21^2*up11+u21*up21*u11;
At2(2,2) = u11*u21*u31;
At2(2,3) = u31*u21^2;
At2(2,4) = -u11^2*u21-u21^3;
At2(3,1) = -u31*up11*u21+u31*up21*u11;
At2(3,2) = u11*u31^2;
At2(3,3) = u21*u31^2;
At2(3,4) = -u31*u21^2-u11^2*u31;
At2(4,1) = -u12*up12*u22+u12^2*up22;
At2(4,2) = u12^2*u32;
At2(4,3) = u12*u22*u32;
At2(4,4) = -u12^3-u12*u22^2;
At2(5,1) = -u22^2*up12+u22*up22*u12;
At2(5,2) = u12*u22*u32;
At2(5,3) = u32*u22^2;
At2(5,4) = -u12^2*u22-u22^3;
At2(6,1) = -u32*up12*u22+u32*up22*u12;
At2(6,2) = u12*u32^2;
At2(6,3) = u22*u32^2;
At2(6,4) = -u32*u22^2-u12^2*u32;
At2(7,1) = -u13*up13*u23+u13^2*up23;
At2(7,2) = u13^2*u33;
At2(7,3) = u13*u23*u33;
At2(7,4) = -u13^3-u13*u23^2;
At2(8,1) = -u23^2*up13+u23*up23*u13;
At2(8,2) = u13*u23*u33;
At2(8,3) = u33*u23^2;
At2(8,4) = -u13^2*u23-u23^3;
At2(9,1) = -u33*up13*u23+u33*up23*u13;
At2(9,2) = u13*u33^2;
At2(9,3) = u23*u33^2;
At2(9,4) = -u33*u23^2-u13^2*u33;
At2(10,1) = -u14*up14*u24+u14^2*up24;
At2(10,2) = u14^2*u34;
At2(10,3) = u14*u24*u34;
At2(10,4) = -u14^3-u14*u24^2;
At2(11,1) = -u24^2*up14+u24*up24*u14;
At2(11,2) = u14*u24*u34;
At2(11,3) = u34*u24^2;
At2(11,4) = -u14^2*u24-u24^3;
At2(12,1) = -u34*up14*u24+u34*up24*u14;
At2(12,2) = u14*u34^2;
At2(12,3) = u24*u34^2;
At2(12,4) = -u34*u24^2-u14^2*u34;
At2(13,1) = -u15*up15*u25+u15^2*up25;
At2(13,2) = u15^2*u35;
At2(13,3) = u25*u35*u15;
At2(13,4) = -u15^3-u15*u25^2;
At2(14,1) = -u25^2*up15+u25*up25*u15;
At2(14,2) = u25*u35*u15;
At2(14,3) = u35*u25^2;
At2(14,4) = -u15^2*u25-u25^3;
At2(15,1) = -u35*up15*u25+u35*up25*u15;
At2(15,2) = u15*u35^2;
At2(15,3) = u25*u35^2;
At2(15,4) = -u35*u25^2-u15^2*u35;



A1(1,1) = -u31*up21*u11+u21*up31*u11;
A1(1,2) = -u11*u31^2-u11*u21^2;
A1(1,3) = u11^2*u21;
A1(1,4) = u11^2*u31;
A1(2,1) = -u21*up21*u31+up31*u21^2;
A1(2,2) = -u21^3-u21*u31^2;
A1(2,3) = u11*u21^2;
A1(2,4) = u11*u21*u31;
A1(3,1) = -up21*u31^2+u31*up31*u21;
A1(3,2) = -u31*u21^2-u31^3;
A1(3,3) = u11*u21*u31;
A1(3,4) = u11*u31^2;
A1(4,1) = -u32*up22*u12+u22*up32*u12;
A1(4,2) = -u12*u32^2-u12*u22^2;
A1(4,3) = u12^2*u22;
A1(4,4) = u12^2*u32;
A1(5,1) = -u22*up22*u32+up32*u22^2;
A1(5,2) = -u22^3-u22*u32^2;
A1(5,3) = u12*u22^2;
A1(5,4) = u12*u22*u32;
A1(6,1) = -up22*u32^2+u32*up32*u22;
A1(6,2) = -u32*u22^2-u32^3;
A1(6,3) = u12*u22*u32;
A1(6,4) = u12*u32^2;
A1(7,1) = -u33*up23*u13+u23*up33*u13;
A1(7,2) = -u13*u33^2-u13*u23^2;
A1(7,3) = u13^2*u23;
A1(7,4) = u13^2*u33;
A1(8,1) = -u23*up23*u33+up33*u23^2;
A1(8,2) = -u23^3-u23*u33^2;
A1(8,3) = u13*u23^2;
A1(8,4) = u13*u23*u33;
A1(9,1) = -up23*u33^2+u33*up33*u23;
A1(9,2) = -u33*u23^2-u33^3;
A1(9,3) = u13*u23*u33;
A1(9,4) = u13*u33^2;
A1(10,1) = -u34*up24*u14+u24*up34*u14;
A1(10,2) = -u14*u34^2-u14*u24^2;
A1(10,3) = u14^2*u24;
A1(10,4) = u14^2*u34;
A1(11,1) = -u24*up24*u34+up34*u24^2;
A1(11,2) = -u24^3-u24*u34^2;
A1(11,3) = u14*u24^2;
A1(11,4) = u14*u24*u34;
A1(12,1) = -up24*u34^2+u34*up34*u24;
A1(12,2) = -u34*u24^2-u34^3;
A1(12,3) = u14*u24*u34;
A1(12,4) = u14*u34^2;
A1(13,1) = -u35*up25*u15+u25*up35*u15;
A1(13,2) = -u15*u35^2-u15*u25^2;
A1(13,3) = u15^2*u25;
A1(13,4) = u15^2*u35;
A1(14,1) = -u25*up25*u35+up35*u25^2;
A1(14,2) = -u25^3-u25*u35^2;
A1(14,3) = u15*u25^2;
A1(14,4) = u25*u35*u15;
A1(15,1) = -up25*u35^2+u35*up35*u25;
A1(15,2) = -u35*u25^2-u35^3;
A1(15,3) = u25*u35*u15;
A1(15,4) = u15*u35^2;