% $$$ for i=1:15
% $$$    for j=1:4
% $$$       fprintf('A(%i,%i) = %s;\n',i,j,char(A(i,j)));
% $$$    end
% $$$ end
A(1,1) = -u31*up21*u11-u11*t2*up11*u21+u11^2*t2*up21-u11^2*t1*up31+u11*t1*up11*u31+u21*up31*u11;
A(1,2) = -u11*u31^2+u11^2*t2*u31+u11^2*t1*u21-u11*u21^2;
A(1,3) = u11*t2*u21*u31-u11*t1*u31^2-u11^3*t1+u21*u11^2;
A(1,4) = u31*u11^2-u11^3*t2-u11*t2*u21^2+u11*t1*u31*u21;
A(2,1) = -u11*t1*up31*u21-u21*up21*u31+up31*u21^2-u21^2*t2*up11+u21*t2*up21*u11+u31*t1*up11*u21;
A(2,2) = u11*t1*u21^2-u21^3-u21*u31^2+u11*t2*u21*u31;
A(2,3) = -u11^2*t1*u21+u11*u21^2+u21^2*t2*u31-u31^2*t1*u21;
A(2,4) = u21*u31*u11-u21*t2*u11^2-u21^3*t2+u31*t1*u21^2;
A(3,1) = -u21*t2*up11*u31-u31*t1*up31*u11+u31^2*t1*up11-up21*u31^2+u31*up31*u21+u11*t2*up21*u31;
A(3,2) = u11*t1*u31*u21-u31*u21^2-u31^3+u11*t2*u31^2;
A(3,3) = u21*t2*u31^2-u31^3*t1-u31*t1*u11^2+u21*u31*u11;
A(3,4) = -u21^2*t2*u31+u31^2*t1*u21+u11*u31^2-u11^2*t2*u31;
A(4,1) = -u32*up22*u12-u12*t2*up12*u22+u12^2*t2*up22-u12^2*t1*up32+u12*t1*up12*u32+u22*up32*u12;
A(4,2) = -u12*u32^2+u12^2*t2*u32+u12^2*t1*u22-u12*u22^2;
A(4,3) = u12*t2*u22*u32-u12*t1*u32^2-u12^3*t1+u22*u12^2;
A(4,4) = u32*u12^2-u12^3*t2-u12*t2*u22^2+u12*t1*u32*u22;
A(5,1) = -u12*t1*up32*u22-u22*up22*u32+up32*u22^2-u22^2*t2*up12+u22*t2*up22*u12+u32*t1*up12*u22;
A(5,2) = u12*t1*u22^2-u22^3-u22*u32^2+u12*t2*u22*u32;
A(5,3) = -u12^2*t1*u22+u12*u22^2+u22^2*t2*u32-u32^2*t1*u22;
A(5,4) = u22*u32*u12-u22*t2*u12^2-u22^3*t2+u32*t1*u22^2;
A(6,1) = -u22*t2*up12*u32-u32*t1*up32*u12+u32^2*t1*up12-up22*u32^2+u32*up32*u22+u12*t2*up22*u32;
A(6,2) = u12*t1*u32*u22-u32*u22^2-u32^3+u12*t2*u32^2;
A(6,3) = u22*t2*u32^2-u32^3*t1-u32*t1*u12^2+u22*u32*u12;
A(6,4) = -u22^2*t2*u32+u32^2*t1*u22+u12*u32^2-u12^2*t2*u32;
A(7,1) = -u33*up23*u13-u13*t2*up13*u23+u13^2*t2*up23-u13^2*t1*up33+u13*t1*up13*u33+u23*up33*u13;
A(7,2) = -u13*u33^2+u13^2*t2*u33+u13^2*t1*u23-u13*u23^2;
A(7,3) = u13*t2*u23*u33-u13*t1*u33^2-u13^3*t1+u23*u13^2;
A(7,4) = u33*u13^2-u13^3*t2-u13*t2*u23^2+u13*t1*u33*u23;
A(8,1) = -u13*t1*up33*u23-u23*up23*u33+up33*u23^2-u23^2*t2*up13+u23*t2*up23*u13+u33*t1*up13*u23;
A(8,2) = u13*t1*u23^2-u23^3-u23*u33^2+u13*t2*u23*u33;
A(8,3) = -u13^2*t1*u23+u13*u23^2+u23^2*t2*u33-u33^2*t1*u23;
A(8,4) = u23*u33*u13-u23*t2*u13^2-u23^3*t2+u33*t1*u23^2;
A(9,1) = -u23*t2*up13*u33-u33*t1*up33*u13+u33^2*t1*up13-up23*u33^2+u33*up33*u23+u13*t2*up23*u33;
A(9,2) = u13*t1*u33*u23-u33*u23^2-u33^3+u13*t2*u33^2;
A(9,3) = u23*t2*u33^2-u33^3*t1-u33*t1*u13^2+u23*u33*u13;
A(9,4) = -u23^2*t2*u33+u33^2*t1*u23+u13*u33^2-u13^2*t2*u33;
A(10,1) = -u34*up24*u14-u14*t2*up14*u24+u14^2*t2*up24-u14^2*t1*up34+u14*t1*up14*u34+u24*up34*u14;
A(10,2) = -u14*u34^2+u14^2*t2*u34+u14^2*t1*u24-u14*u24^2;
A(10,3) = u14*t2*u24*u34-u14*t1*u34^2-u14^3*t1+u24*u14^2;
A(10,4) = u34*u14^2-u14^3*t2-u14*t2*u24^2+u14*t1*u34*u24;
A(11,1) = -u14*t1*up34*u24-u24*up24*u34+up34*u24^2-u24^2*t2*up14+u24*t2*up24*u14+u34*t1*up14*u24;
A(11,2) = u14*t1*u24^2-u24^3-u24*u34^2+u14*t2*u24*u34;
A(11,3) = -u14^2*t1*u24+u14*u24^2+u24^2*t2*u34-u34^2*t1*u24;
A(11,4) = u24*u34*u14-u24*t2*u14^2-u24^3*t2+u34*t1*u24^2;
A(12,1) = -u24*t2*up14*u34-u34*t1*up34*u14+u34^2*t1*up14-up24*u34^2+u34*up34*u24+u14*t2*up24*u34;
A(12,2) = u14*t1*u34*u24-u34*u24^2-u34^3+u14*t2*u34^2;
A(12,3) = u24*t2*u34^2-u34^3*t1-u34*t1*u14^2+u24*u34*u14;
A(12,4) = -u24^2*t2*u34+u34^2*t1*u24+u14*u34^2-u14^2*t2*u34;
A(13,1) = -u35*up25*u15-u15*t2*up15*u25+u15^2*t2*up25-u15^2*t1*up35+u15*t1*up15*u35+u25*up35*u15;
A(13,2) = -u15*u35^2+u15^2*t2*u35+u15^2*t1*u25-u15*u25^2;
A(13,3) = u15*t2*u25*u35-u15*t1*u35^2-u15^3*t1+u25*u15^2;
A(13,4) = u35*u15^2-u15^3*t2-u15*t2*u25^2+u15*t1*u35*u25;
A(14,1) = -u15*t1*up35*u25-u25*up25*u35+up35*u25^2-u25^2*t2*up15+u25*t2*up25*u15+u35*t1*up15*u25;
A(14,2) = u15*t1*u25^2-u25^3-u25*u35^2+u15*t2*u25*u35;
A(14,3) = -u15^2*t1*u25+u15*u25^2+u25^2*t2*u35-u35^2*t1*u25;
A(14,4) = u25*u35*u15-u25*t2*u15^2-u25^3*t2+u35*t1*u25^2;
A(15,1) = -u25*t2*up15*u35-u35*t1*up35*u15+u35^2*t1*up15-up25*u35^2+u35*up35*u25+u15*t2*up25*u35;
A(15,2) = u15*t1*u35*u25-u35*u25^2-u35^3+u15*t2*u35^2;
A(15,3) = u25*t2*u35^2-u35^3*t1-u35*t1*u15^2+u25*u35*u15;
A(15,4) = -u25^2*t2*u35+u35^2*t1*u25+u15*u35^2-u15^2*t2*u35;