.
.
.
.
.
.
.
.
.
.
.
%function [bv00,bv0,bv10,bv1,bv20,bv2,bv30,bv3,bv40,bv4,bv50,bv5] = CalculateDriftbCoeffs08A2Dim(wMu0dt,dwMu0dtdw,c,SeriesOrder)
function [b] = CalculateDriftbCoeffs08A2Dim02(wMu0dt,dwMu0dtdw,c,SeriesOrder)
b(1,1)=wMu0dt;
b(1,2)=c(1,2) *dwMu0dtdw(1);
b(1,3)=1/2*(2* c(1,3) *dwMu0dtdw(1)+c(1,2)^2 *dwMu0dtdw(2));
b(1,4)=1/6*(6 *c(1,4) *dwMu0dtdw(1)+6 *c(1,2)* c(1,3) *dwMu0dtdw(2)+c(1,2)^3 *dwMu0dtdw(3));
b(1,5)=1/24*(24* c(1,5) *dwMu0dtdw(1)+12* c(1,3)^2 *dwMu0dtdw(2)+24 *c(1,2)* c(1,4) *dwMu0dtdw(2)+ ...
12* c(1,2)^2* c(1,3) *dwMu0dtdw(3)+c(1,2)^4 *dwMu0dtdw(4));
b(1,6)=1/120*(120* c(1,6) *dwMu0dtdw(1)+120* c(1,3) *c(1,4) *dwMu0dtdw(2)+120* c(1,2)* c(1,5) *dwMu0dtdw(2)+ ...
60 *c(1,2)* c(1,3)^2 *dwMu0dtdw(3)+60* c(1,2)^2* c(1,4) *dwMu0dtdw(3)+ ...
20 *c(1,2)^3* c(1,3) *dwMu0dtdw(4)+c(1,2)^5 *dwMu0dtdw(5));
b(2,1)=c(2,1) *dwMu0dtdw(1);
b(3,1)=1/2*(2* c(3,1) *dwMu0dtdw(1)+c(2,1)^2 *dwMu0dtdw(2));
b(4,1)=1/6*(6 *c(4,1) *dwMu0dtdw(1)+6 *c(2,1)* c(3,1) *dwMu0dtdw(2)+c(2,1)^3 *dwMu0dtdw(3));
b(5,1)=1/24*(24* c(5,1) *dwMu0dtdw(1)+12* c(3,1)^2 *dwMu0dtdw(2)+24 *c(2,1)* c(4,1) *dwMu0dtdw(2)+ ...
12* c(2,1)^2* c(3,1) *dwMu0dtdw(3)+c(2,1)^4 *dwMu0dtdw(4));
b(6,1)=1/120*(120* c(6,1) *dwMu0dtdw(1)+120* c(3,1) *c(4,1) *dwMu0dtdw(2)+120* c(2,1)* c(5,1) *dwMu0dtdw(2)+ ...
60 *c(2,1)* c(3,1)^2 *dwMu0dtdw(3)+60* c(2,1)^2* c(4,1) *dwMu0dtdw(3)+ ...
20 *c(2,1)^3* c(3,1) *dwMu0dtdw(4)+c(2,1)^5 *dwMu0dtdw(5));
b(2,2)=dwMu0dtdw(2).*c(1,2).*c(2,1)+dwMu0dtdw(1).*c(2,2);
b(3,2)=1/2*(dwMu0dtdw(3).*c(1,2).*c(2,1).^2+2.*dwMu0dtdw(2).*c(2,1).*c(2,2)+dwMu0dtdw(2).*c(1,2).*2.*c(3,1)+dwMu0dtdw(1).*2.*c(3,2));
b(4,2)=1/6*(dwMu0dtdw(4).*c(1,2).*c(2,1).^3+3.*dwMu0dtdw(3).*c(2,1).^2.*c(2,2)+3.*dwMu0dtdw(3).*c(1,2).*c(2,1).*2.*c(3,1)+3.*dwMu0dtdw(2).*c(2,2).*2.*c(3,1)+3.*dwMu0dtdw(2).*c(2,1).*2.*c(3,2)+dwMu0dtdw(2).*c(1,2).*6.*c(4,1)+dwMu0dtdw(1).*6.*c(4,2));
b(5,2)=1/24*(dwMu0dtdw(5).*c(1,2).*c(2,1).^4+4.*dwMu0dtdw(4).*c(2,1).^3.*c(2,2)+6.*dwMu0dtdw(4).*c(1,2).*c(2,1).^2.*2.*c(3,1)+12.*dwMu0dtdw(3).*c(2,1).*c(2,2).*2.*c(3,1)+3.*dwMu0dtdw(3).*c(1,2).*2.*c(3,1).^2+6.*dwMu0dtdw(3).*c(2,1).^2.*2.*c(3,2)+6.*dwMu0dtdw(2).*2.*c(3,1).*2.*c(3,2)+4.*dwMu0dtdw(3).*c(1,2).*c(2,1).*6.*c(4,1)+4.*dwMu0dtdw(2).*c(2,2).*6.*c(4,1)+4.*dwMu0dtdw(2).*c(2,1).*6.*c(4,2)+dwMu0dtdw(2).*c(1,2).*24.*c(5,1)+dwMu0dtdw(1).*24.*c(5,2));
b(6,2)=1/120*(dwMu0dtdw(6).*c(1,2).*c(2,1).^5+5.*dwMu0dtdw(5).*c(2,1).^4.*c(2,2)+10.*dwMu0dtdw(5).*c(1,2).*c(2,1).^3.*2.*c(3,1)+30.*dwMu0dtdw(4).*c(2,1).^2.*c(2,2).*2.*c(3,1)+15.*dwMu0dtdw(4).*c(1,2).*c(2,1).*2.*c(3,1).^2+15.*dwMu0dtdw(3).*c(2,2).*2.*c(3,1).^2+10.*dwMu0dtdw(4).*c(2,1).^3.*2.*c(3,2)+30.*dwMu0dtdw(3).*c(2,1).*2.*c(3,1).*2.*c(3,2)+10.*dwMu0dtdw(4).*c(1,2).*c(2,1).^2.*6.*c(4,1)+20.*dwMu0dtdw(3).*c(2,1).*c(2,2).*6.*c(4,1)+10.*dwMu0dtdw(3).*c(1,2).*2.*c(3,1).*6.*c(4,1)+10.*dwMu0dtdw(2).*2.*c(3,2).*6.*c(4,1)+10.*dwMu0dtdw(3).*c(2,1).^2.*6.*c(4,2)+10.*dwMu0dtdw(2).*2.*c(3,1).*6.*c(4,2)+5.*dwMu0dtdw(3).*c(1,2).*c(2,1).*24.*c(5,1)+5.*dwMu0dtdw(2).*c(2,2).*24.*c(5,1)+5.*dwMu0dtdw(2).*c(2,1).*24.*c(5,2)+dwMu0dtdw(2).*c(1,2).*120.*c(6,1)+dwMu0dtdw(1).*120.*c(6,2));
b(2,3)=1/2*((dwMu0dtdw(3).*c(1,2).^2+dwMu0dtdw(2).*2.*c(1,3)).*c(2,1)+2.*dwMu0dtdw(2).*c(1,2).*c(2,2)+dwMu0dtdw(1).*2.*c(2,3));
b(3,3)=1/4*((dwMu0dtdw(4).*c(1,2).^2+dwMu0dtdw(3).*2.*c(1,3)).*c(2,1).^2+4.*dwMu0dtdw(3).*c(1,2).*c(2,1).*c(2,2)+dwMu0dtdw(2).*(2.*c(2,2).^2+2.*c(2,1).*2.*c(2,3))+(dwMu0dtdw(3).*c(1,2).^2+dwMu0dtdw(2).*2.*c(1,3)).*2.*c(3,1)+2.*dwMu0dtdw(2).*c(1,2).*2.*c(3,2)+dwMu0dtdw(1).*4.*c(3,3));
b(4,3)=1/12*((dwMu0dtdw(5).*c(1,2).^2+dwMu0dtdw(4).*2.*c(1,3)).*c(2,1).^3+6.*dwMu0dtdw(4).*c(1,2).*c(2,1).^2.*c(2,2)+dwMu0dtdw(3).*(6.*c(2,1).*c(2,2).^2+3.*c(2,1).^2.*2.*c(2,3))+3.*(dwMu0dtdw(4).*c(1,2).^2+dwMu0dtdw(3).*2.*c(1,3)).*c(2,1).*2.*c(3,1)+6.*dwMu0dtdw(3).*c(1,2).*(c(2,2).*2.*c(3,1)+c(2,1).*2.*c(3,2))+3.*dwMu0dtdw(2).*(2.*c(2,3).*2.*c(3,1)+2.*c(2,2).*2.*c(3,2)+c(2,1).*4.*c(3,3))+(dwMu0dtdw(3).*c(1,2).^2+dwMu0dtdw(2).*2.*c(1,3)).*6.*c(4,1)+2.*dwMu0dtdw(2).*c(1,2).*6.*c(4,2)+dwMu0dtdw(1).*12.*c(4,3));
b(5,3)=1/48*((dwMu0dtdw(6).*c(1,2).^2+dwMu0dtdw(5).*2.*c(1,3)).*c(2,1).^4+8.*dwMu0dtdw(5).*c(1,2).*c(2,1).^3.*c(2,2)+dwMu0dtdw(4).*(12.*c(2,1).^2.*c(2,2).^2+4.*c(2,1).^3.*2.*c(2,3))+6.*(dwMu0dtdw(5).*c(1,2).^2+dwMu0dtdw(4).*2.*c(1,3)).*c(2,1).^2.*2.*c(3,1)+3.*(dwMu0dtdw(4).*c(1,2).^2+dwMu0dtdw(3).*2.*c(1,3)).*2.*c(3,1).^2+12.*dwMu0dtdw(3).*c(1,2).*2.*c(3,1).*2.*c(3,2)+12.*dwMu0dtdw(4).*c(1,2).*(2.*c(2,1).*c(2,2).*2.*c(3,1)+c(2,1).^2.*2.*c(3,2))+6.*dwMu0dtdw(3).*((2.*c(2,2).^2+2.*c(2,1).*2.*c(2,3)).*2.*c(3,1)+4.*c(2,1).*c(2,2).*2.*c(3,2)+c(2,1).^2.*4.*c(3,3))+3.*dwMu0dtdw(2).*(2.*2.*c(3,2).^2+2.*2.*c(3,1).*4.*c(3,3))+4.*(dwMu0dtdw(4).*c(1,2).^2+dwMu0dtdw(3).*2.*c(1,3)).*c(2,1).*6.*c(4,1)+8.*dwMu0dtdw(3).*c(1,2).*(c(2,2).*6.*c(4,1)+c(2,1).*6.*c(4,2))+4.*dwMu0dtdw(2).*(2.*c(2,3).*6.*c(4,1)+2.*c(2,2).*6.*c(4,2)+c(2,1).*12.*c(4,3))+(dwMu0dtdw(3).*c(1,2).^2+dwMu0dtdw(2).*2.*c(1,3)).*24.*c(5,1)+2.*dwMu0dtdw(2).*c(1,2).*24.*c(5,2)+dwMu0dtdw(1).*48.*c(5,3));
b(6,3)=1/240*((dwMu0dtdw(7).*c(1,2).^2+dwMu0dtdw(6).*2.*c(1,3)).*c(2,1).^5+10.*dwMu0dtdw(6).*c(1,2).*c(2,1).^4.*c(2,2)+dwMu0dtdw(5).*(20.*c(2,1).^3.*c(2,2).^2+5.*c(2,1).^4.*2.*c(2,3))+10.*(dwMu0dtdw(6).*c(1,2).^2+dwMu0dtdw(5).*2.*c(1,3)).*c(2,1).^3.*2.*c(3,1)+15.*(dwMu0dtdw(5).*c(1,2).^2+dwMu0dtdw(4).*2.*c(1,3)).*c(2,1).*2.*c(3,1).^2+20.*dwMu0dtdw(5).*c(1,2).*(3.*c(2,1).^2.*c(2,2).*2.*c(3,1)+c(2,1).^3.*2.*c(3,2))+30.*dwMu0dtdw(4).*c(1,2).*(c(2,2).*2.*c(3,1).^2+2.*c(2,1).*2.*c(3,1).*2.*c(3,2))+10.*dwMu0dtdw(4).*((6.*c(2,1).*c(2,2).^2+3.*c(2,1).^2.*2.*c(2,3)).*2.*c(3,1)+6.*c(2,1).^2.*c(2,2).*2.*c(3,2)+c(2,1).^3.*4.*c(3,3))+15.*dwMu0dtdw(3).*(2.*c(2,3).*2.*c(3,1).^2+4.*c(2,2).*2.*c(3,1).*2.*c(3,2)+c(2,1).*(2.*2.*c(3,2).^2+2.*2.*c(3,1).*4.*c(3,3)))+10.*(dwMu0dtdw(5).*c(1,2).^2+dwMu0dtdw(4).*2.*c(1,3)).*c(2,1).^2.*6.*c(4,1)+10.*(dwMu0dtdw(4).*c(1,2).^2+dwMu0dtdw(3).*2.*c(1,3)).*2.*c(3,1).*6.*c(4,1)+20.*dwMu0dtdw(4).*c(1,2).*(2.*c(2,1).*c(2,2).*6.*c(4,1)+c(2,1).^2.*6.*c(4,2))+20.*dwMu0dtdw(3).*c(1,2).*(2.*c(3,2).*6.*c(4,1)+2.*c(3,1).*6.*c(4,2))+10.*dwMu0dtdw(3).*((2.*c(2,2).^2+2.*c(2,1).*2.*c(2,3)).*6.*c(4,1)+4.*c(2,1).*c(2,2).*6.*c(4,2)+c(2,1).^2.*12.*c(4,3))+10.*dwMu0dtdw(2).*(4.*c(3,3).*6.*c(4,1)+2.*2.*c(3,2).*6.*c(4,2)+2.*c(3,1).*12.*c(4,3))+5.*(dwMu0dtdw(4).*c(1,2).^2+dwMu0dtdw(3).*2.*c(1,3)).*c(2,1).*24.*c(5,1)+10.*dwMu0dtdw(3).*c(1,2).*(c(2,2).*24.*c(5,1)+c(2,1).*24.*c(5,2))+5.*dwMu0dtdw(2).*(2.*c(2,3).*24.*c(5,1)+2.*c(2,2).*24.*c(5,2)+c(2,1).*48.*c(5,3))+(dwMu0dtdw(3).*c(1,2).^2+dwMu0dtdw(2).*2.*c(1,3)).*120.*c(6,1)+2.*dwMu0dtdw(2).*c(1,2).*120.*c(6,2)+dwMu0dtdw(1).*240.*c(6,3));
b(2,4)=1/6*((dwMu0dtdw(4).*c(1,2).^3+3.*dwMu0dtdw(3).*c(1,2).*2.*c(1,3)+dwMu0dtdw(2).*6.*c(1,4)).*c(2,1)+3.*(dwMu0dtdw(3).*c(1,2).^2+dwMu0dtdw(2).*2.*c(1,3)).*c(2,2)+3.*dwMu0dtdw(2).*c(1,2).*2.*c(2,3)+dwMu0dtdw(1).*6.*c(2,4));
b(3,4)=1/12*((dwMu0dtdw(5).*c(1,2).^3+3.*dwMu0dtdw(4).*c(1,2).*2.*c(1,3)+dwMu0dtdw(3).*6.*c(1,4)).*c(2,1).^2+6.*(dwMu0dtdw(4).*c(1,2).^2+dwMu0dtdw(3).*2.*c(1,3)).*c(2,1).*c(2,2)+3.*dwMu0dtdw(3).*c(1,2).*(2.*c(2,2).^2+2.*c(2,1).*2.*c(2,3))+dwMu0dtdw(2).*(6.*c(2,2).*2.*c(2,3)+2.*c(2,1).*6.*c(2,4))+(dwMu0dtdw(4).*c(1,2).^3+3.*dwMu0dtdw(3).*c(1,2).*2.*c(1,3)+dwMu0dtdw(2).*6.*c(1,4)).*2.*c(3,1)+3.*(dwMu0dtdw(3).*c(1,2).^2+dwMu0dtdw(2).*2.*c(1,3)).*2.*c(3,2)+3.*dwMu0dtdw(2).*c(1,2).*4.*c(3,3)+dwMu0dtdw(1).*12.*c(3,4));
b(4,4)=1/36*((dwMu0dtdw(6).*c(1,2).^3+3.*dwMu0dtdw(5).*c(1,2).*2.*c(1,3)+dwMu0dtdw(4).*6.*c(1,4)).*c(2,1).^3+9.*(dwMu0dtdw(5).*c(1,2).^2+dwMu0dtdw(4).*2.*c(1,3)).*c(2,1).^2.*c(2,2)+3.*dwMu0dtdw(4).*c(1,2).*(6.*c(2,1).*c(2,2).^2+3.*c(2,1).^2.*2.*c(2,3))+dwMu0dtdw(3).*(6.*c(2,2).^3+18.*c(2,1).*c(2,2).*2.*c(2,3)+3.*c(2,1).^2.*6.*c(2,4))+3.*(dwMu0dtdw(5).*c(1,2).^3+3.*dwMu0dtdw(4).*c(1,2).*2.*c(1,3)+dwMu0dtdw(3).*6.*c(1,4)).*c(2,1).*2.*c(3,1)+9.*(dwMu0dtdw(4).*c(1,2).^2+dwMu0dtdw(3).*2.*c(1,3)).*(c(2,2).*2.*c(3,1)+c(2,1).*2.*c(3,2))+9.*dwMu0dtdw(3).*c(1,2).*(2.*c(2,3).*2.*c(3,1)+2.*c(2,2).*2.*c(3,2)+c(2,1).*4.*c(3,3))+3.*dwMu0dtdw(2).*(6.*c(2,4).*2.*c(3,1)+3.*2.*c(2,3).*2.*c(3,2)+3.*c(2,2).*4.*c(3,3)+c(2,1).*12.*c(3,4))+(dwMu0dtdw(4).*c(1,2).^3+3.*dwMu0dtdw(3).*c(1,2).*2.*c(1,3)+dwMu0dtdw(2).*6.*c(1,4)).*6.*c(4,1)+3.*(dwMu0dtdw(3).*c(1,2).^2+dwMu0dtdw(2).*2.*c(1,3)).*6.*c(4,2)+3.*dwMu0dtdw(2).*c(1,2).*12.*c(4,3)+dwMu0dtdw(1).*36.*c(4,4));
b(5,4)=1/144*((dwMu0dtdw(7).*c(1,2).^3+3.*dwMu0dtdw(6).*c(1,2).*2.*c(1,3)+dwMu0dtdw(5).*6.*c(1,4)).*c(2,1).^4+12.*(dwMu0dtdw(6).*c(1,2).^2+dwMu0dtdw(5).*2.*c(1,3)).*c(2,1).^3.*c(2,2)+3.*dwMu0dtdw(5).*c(1,2).*(12.*c(2,1).^2.*c(2,2).^2+4.*c(2,1).^3.*2.*c(2,3))+dwMu0dtdw(4).*(24.*c(2,1).*c(2,2).^3+36.*c(2,1).^2.*c(2,2).*2.*c(2,3)+4.*c(2,1).^3.*6.*c(2,4))+6.*(dwMu0dtdw(6).*c(1,2).^3+3.*dwMu0dtdw(5).*c(1,2).*2.*c(1,3)+dwMu0dtdw(4).*6.*c(1,4)).*c(2,1).^2.*2.*c(3,1)+3.*(dwMu0dtdw(5).*c(1,2).^3+3.*dwMu0dtdw(4).*c(1,2).*2.*c(1,3)+dwMu0dtdw(3).*6.*c(1,4)).*2.*c(3,1).^2+18.*(dwMu0dtdw(4).*c(1,2).^2+dwMu0dtdw(3).*2.*c(1,3)).*2.*c(3,1).*2.*c(3,2)+18.*(dwMu0dtdw(5).*c(1,2).^2+dwMu0dtdw(4).*2.*c(1,3)).*(2.*c(2,1).*c(2,2).*2.*c(3,1)+c(2,1).^2.*2.*c(3,2))+18.*dwMu0dtdw(4).*c(1,2).*((2.*c(2,2).^2+2.*c(2,1).*2.*c(2,3)).*2.*c(3,1)+4.*c(2,1).*c(2,2).*2.*c(3,2)+c(2,1).^2.*4.*c(3,3))+9.*dwMu0dtdw(3).*c(1,2).*(2.*2.*c(3,2).^2+2.*2.*c(3,1).*4.*c(3,3))+6.*dwMu0dtdw(3).*((6.*c(2,2).*2.*c(2,3)+2.*c(2,1).*6.*c(2,4)).*2.*c(3,1)+3.*(2.*c(2,2).^2+2.*c(2,1).*2.*c(2,3)).*2.*c(3,2)+6.*c(2,1).*c(2,2).*4.*c(3,3)+c(2,1).^2.*12.*c(3,4))+3.*dwMu0dtdw(2).*(6.*2.*c(3,2).*4.*c(3,3)+2.*2.*c(3,1).*12.*c(3,4))+4.*(dwMu0dtdw(5).*c(1,2).^3+3.*dwMu0dtdw(4).*c(1,2).*2.*c(1,3)+dwMu0dtdw(3).*6.*c(1,4)).*c(2,1).*6.*c(4,1)+12.*(dwMu0dtdw(4).*c(1,2).^2+dwMu0dtdw(3).*2.*c(1,3)).*(c(2,2).*6.*c(4,1)+c(2,1).*6.*c(4,2))+12.*dwMu0dtdw(3).*c(1,2).*(2.*c(2,3).*6.*c(4,1)+2.*c(2,2).*6.*c(4,2)+c(2,1).*12.*c(4,3))+4.*dwMu0dtdw(2).*(6.*c(2,4).*6.*c(4,1)+3.*2.*c(2,3).*6.*c(4,2)+3.*c(2,2).*12.*c(4,3)+c(2,1).*36.*c(4,4))+(dwMu0dtdw(4).*c(1,2).^3+3.*dwMu0dtdw(3).*c(1,2).*2.*c(1,3)+dwMu0dtdw(2).*6.*c(1,4)).*24.*c(5,1)+3.*(dwMu0dtdw(3).*c(1,2).^2+dwMu0dtdw(2).*2.*c(1,3)).*24.*c(5,2)+3.*dwMu0dtdw(2).*c(1,2).*48.*c(5,3)+dwMu0dtdw(1).*144.*c(5,4));
b(6,4)=1/720*((dwMu0dtdw(8).*c(1,2).^3+3.*dwMu0dtdw(7).*c(1,2).*2.*c(1,3)+dwMu0dtdw(6).*6.*c(1,4)).*c(2,1).^5+15.*(dwMu0dtdw(7).*c(1,2).^2+dwMu0dtdw(6).*2.*c(1,3)).*c(2,1).^4.*c(2,2)+3.*dwMu0dtdw(6).*c(1,2).*(20.*c(2,1).^3.*c(2,2).^2+5.*c(2,1).^4.*2.*c(2,3))+dwMu0dtdw(5).*(60.*c(2,1).^2.*c(2,2).^3+60.*c(2,1).^3.*c(2,2).*2.*c(2,3)+5.*c(2,1).^4.*6.*c(2,4))+10.*(dwMu0dtdw(7).*c(1,2).^3+3.*dwMu0dtdw(6).*c(1,2).*2.*c(1,3)+dwMu0dtdw(5).*6.*c(1,4)).*c(2,1).^3.*2.*c(3,1)+15.*(dwMu0dtdw(6).*c(1,2).^3+3.*dwMu0dtdw(5).*c(1,2).*2.*c(1,3)+dwMu0dtdw(4).*6.*c(1,4)).*c(2,1).*2.*c(3,1).^2+30.*(dwMu0dtdw(6).*c(1,2).^2+dwMu0dtdw(5).*2.*c(1,3)).*(3.*c(2,1).^2.*c(2,2).*2.*c(3,1)+c(2,1).^3.*2.*c(3,2))+45.*(dwMu0dtdw(5).*c(1,2).^2+dwMu0dtdw(4).*2.*c(1,3)).*(c(2,2).*2.*c(3,1).^2+2.*c(2,1).*2.*c(3,1).*2.*c(3,2))+30.*dwMu0dtdw(5).*c(1,2).*((6.*c(2,1).*c(2,2).^2+3.*c(2,1).^2.*2.*c(2,3)).*2.*c(3,1)+6.*c(2,1).^2.*c(2,2).*2.*c(3,2)+c(2,1).^3.*4.*c(3,3))+45.*dwMu0dtdw(4).*c(1,2).*(2.*c(2,3).*2.*c(3,1).^2+4.*c(2,2).*2.*c(3,1).*2.*c(3,2)+c(2,1).*(2.*2.*c(3,2).^2+2.*2.*c(3,1).*4.*c(3,3)))+10.*dwMu0dtdw(4).*((6.*c(2,2).^3+18.*c(2,1).*c(2,2).*2.*c(2,3)+3.*c(2,1).^2.*6.*c(2,4)).*2.*c(3,1)+3.*(6.*c(2,1).*c(2,2).^2+3.*c(2,1).^2.*2.*c(2,3)).*2.*c(3,2)+9.*c(2,1).^2.*c(2,2).*4.*c(3,3)+c(2,1).^3.*12.*c(3,4))+15.*dwMu0dtdw(3).*(6.*c(2,4).*2.*c(3,1).^2+6.*2.*c(2,3).*2.*c(3,1).*2.*c(3,2)+3.*c(2,2).*(2.*2.*c(3,2).^2+2.*2.*c(3,1).*4.*c(3,3))+c(2,1).*(6.*2.*c(3,2).*4.*c(3,3)+2.*2.*c(3,1).*12.*c(3,4)))+10.*(dwMu0dtdw(6).*c(1,2).^3+3.*dwMu0dtdw(5).*c(1,2).*2.*c(1,3)+dwMu0dtdw(4).*6.*c(1,4)).*c(2,1).^2.*6.*c(4,1)+10.*(dwMu0dtdw(5).*c(1,2).^3+3.*dwMu0dtdw(4).*c(1,2).*2.*c(1,3)+dwMu0dtdw(3).*6.*c(1,4)).*2.*c(3,1).*6.*c(4,1)+30.*(dwMu0dtdw(5).*c(1,2).^2+dwMu0dtdw(4).*2.*c(1,3)).*(2.*c(2,1).*c(2,2).*6.*c(4,1)+c(2,1).^2.*6.*c(4,2))+30.*(dwMu0dtdw(4).*c(1,2).^2+dwMu0dtdw(3).*2.*c(1,3)).*(2.*c(3,2).*6.*c(4,1)+2.*c(3,1).*6.*c(4,2))+30.*dwMu0dtdw(4).*c(1,2).*((2.*c(2,2).^2+2.*c(2,1).*2.*c(2,3)).*6.*c(4,1)+4.*c(2,1).*c(2,2).*6.*c(4,2)+c(2,1).^2.*12.*c(4,3))+30.*dwMu0dtdw(3).*c(1,2).*(4.*c(3,3).*6.*c(4,1)+2.*2.*c(3,2).*6.*c(4,2)+2.*c(3,1).*12.*c(4,3))+10.*dwMu0dtdw(3).*((6.*c(2,2).*2.*c(2,3)+2.*c(2,1).*6.*c(2,4)).*6.*c(4,1)+3.*(2.*c(2,2).^2+2.*c(2,1).*2.*c(2,3)).*6.*c(4,2)+6.*c(2,1).*c(2,2).*12.*c(4,3)+c(2,1).^2.*36.*c(4,4))+10.*dwMu0dtdw(2).*(12.*c(3,4).*6.*c(4,1)+3.*4.*c(3,3).*6.*c(4,2)+3.*2.*c(3,2).*12.*c(4,3)+2.*c(3,1).*36.*c(4,4))+5.*(dwMu0dtdw(5).*c(1,2).^3+3.*dwMu0dtdw(4).*c(1,2).*2.*c(1,3)+dwMu0dtdw(3).*6.*c(1,4)).*c(2,1).*24.*c(5,1)+15.*(dwMu0dtdw(4).*c(1,2).^2+dwMu0dtdw(3).*2.*c(1,3)).*(c(2,2).*24.*c(5,1)+c(2,1).*24.*c(5,2))+15.*dwMu0dtdw(3).*c(1,2).*(2.*c(2,3).*24.*c(5,1)+2.*c(2,2).*24.*c(5,2)+c(2,1).*48.*c(5,3))+5.*dwMu0dtdw(2).*(6.*c(2,4).*24.*c(5,1)+3.*2.*c(2,3).*24.*c(5,2)+3.*c(2,2).*48.*c(5,3)+c(2,1).*144.*c(5,4))+(dwMu0dtdw(4).*c(1,2).^3+3.*dwMu0dtdw(3).*c(1,2).*2.*c(1,3)+dwMu0dtdw(2).*6.*c(1,4)).*120.*c(6,1)+3.*(dwMu0dtdw(3).*c(1,2).^2+dwMu0dtdw(2).*2.*c(1,3)).*120.*c(6,2)+3.*dwMu0dtdw(2).*c(1,2).*240.*c(6,3)+dwMu0dtdw(1).*720.*c(6,4));
b(2,5)=1/24*((dwMu0dtdw(5).*c(1,2).^4+6.*dwMu0dtdw(4).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(3).*2.*c(1,3).^2+4.*dwMu0dtdw(3).*c(1,2).*6.*c(1,4)+dwMu0dtdw(2).*24.*c(1,5)).*c(2,1)+4.*(dwMu0dtdw(4).*c(1,2).^3+3.*dwMu0dtdw(3).*c(1,2).*2.*c(1,3)+dwMu0dtdw(2).*6.*c(1,4)).*c(2,2)+6.*(dwMu0dtdw(3).*c(1,2).^2+dwMu0dtdw(2).*2.*c(1,3)).*2.*c(2,3)+4.*dwMu0dtdw(2).*c(1,2).*6.*c(2,4)+dwMu0dtdw(1).*24.*c(2,5));
b(3,5)=1/48*((dwMu0dtdw(6).*c(1,2).^4+6.*dwMu0dtdw(5).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(4).*2.*c(1,3).^2+4.*dwMu0dtdw(4).*c(1,2).*6.*c(1,4)+dwMu0dtdw(3).*24.*c(1,5)).*c(2,1).^2+8.*(dwMu0dtdw(5).*c(1,2).^3+3.*dwMu0dtdw(4).*c(1,2).*2.*c(1,3)+dwMu0dtdw(3).*6.*c(1,4)).*c(2,1).*c(2,2)+6.*(dwMu0dtdw(4).*c(1,2).^2+dwMu0dtdw(3).*2.*c(1,3)).*(2.*c(2,2).^2+2.*c(2,1).*2.*c(2,3))+4.*dwMu0dtdw(3).*c(1,2).*(6.*c(2,2).*2.*c(2,3)+2.*c(2,1).*6.*c(2,4))+dwMu0dtdw(2).*(6.*2.*c(2,3).^2+8.*c(2,2).*6.*c(2,4)+2.*c(2,1).*24.*c(2,5))+(dwMu0dtdw(5).*c(1,2).^4+6.*dwMu0dtdw(4).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(3).*2.*c(1,3).^2+4.*dwMu0dtdw(3).*c(1,2).*6.*c(1,4)+dwMu0dtdw(2).*24.*c(1,5)).*2.*c(3,1)+4.*(dwMu0dtdw(4).*c(1,2).^3+3.*dwMu0dtdw(3).*c(1,2).*2.*c(1,3)+dwMu0dtdw(2).*6.*c(1,4)).*2.*c(3,2)+6.*(dwMu0dtdw(3).*c(1,2).^2+dwMu0dtdw(2).*2.*c(1,3)).*4.*c(3,3)+4.*dwMu0dtdw(2).*c(1,2).*12.*c(3,4)+dwMu0dtdw(1).*48.*c(3,5));
b(4,5)=1/144*((dwMu0dtdw(7).*c(1,2).^4+6.*dwMu0dtdw(6).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(5).*2.*c(1,3).^2+4.*dwMu0dtdw(5).*c(1,2).*6.*c(1,4)+dwMu0dtdw(4).*24.*c(1,5)).*c(2,1).^3+12.*(dwMu0dtdw(6).*c(1,2).^3+3.*dwMu0dtdw(5).*c(1,2).*2.*c(1,3)+dwMu0dtdw(4).*6.*c(1,4)).*c(2,1).^2.*c(2,2)+6.*(dwMu0dtdw(5).*c(1,2).^2+dwMu0dtdw(4).*2.*c(1,3)).*(6.*c(2,1).*c(2,2).^2+3.*c(2,1).^2.*2.*c(2,3))+4.*dwMu0dtdw(4).*c(1,2).*(6.*c(2,2).^3+18.*c(2,1).*c(2,2).*2.*c(2,3)+3.*c(2,1).^2.*6.*c(2,4))+dwMu0dtdw(3).*(36.*c(2,2).^2.*2.*c(2,3)+18.*c(2,1).*2.*c(2,3).^2+24.*c(2,1).*c(2,2).*6.*c(2,4)+3.*c(2,1).^2.*24.*c(2,5))+3.*(dwMu0dtdw(6).*c(1,2).^4+6.*dwMu0dtdw(5).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(4).*2.*c(1,3).^2+4.*dwMu0dtdw(4).*c(1,2).*6.*c(1,4)+dwMu0dtdw(3).*24.*c(1,5)).*c(2,1).*2.*c(3,1)+12.*(dwMu0dtdw(5).*c(1,2).^3+3.*dwMu0dtdw(4).*c(1,2).*2.*c(1,3)+dwMu0dtdw(3).*6.*c(1,4)).*(c(2,2).*2.*c(3,1)+c(2,1).*2.*c(3,2))+18.*(dwMu0dtdw(4).*c(1,2).^2+dwMu0dtdw(3).*2.*c(1,3)).*(2.*c(2,3).*2.*c(3,1)+2.*c(2,2).*2.*c(3,2)+c(2,1).*4.*c(3,3))+12.*dwMu0dtdw(3).*c(1,2).*(6.*c(2,4).*2.*c(3,1)+3.*2.*c(2,3).*2.*c(3,2)+3.*c(2,2).*4.*c(3,3)+c(2,1).*12.*c(3,4))+3.*dwMu0dtdw(2).*(24.*c(2,5).*2.*c(3,1)+4.*6.*c(2,4).*2.*c(3,2)+6.*2.*c(2,3).*4.*c(3,3)+4.*c(2,2).*12.*c(3,4)+c(2,1).*48.*c(3,5))+(dwMu0dtdw(5).*c(1,2).^4+6.*dwMu0dtdw(4).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(3).*2.*c(1,3).^2+4.*dwMu0dtdw(3).*c(1,2).*6.*c(1,4)+dwMu0dtdw(2).*24.*c(1,5)).*6.*c(4,1)+4.*(dwMu0dtdw(4).*c(1,2).^3+3.*dwMu0dtdw(3).*c(1,2).*2.*c(1,3)+dwMu0dtdw(2).*6.*c(1,4)).*6.*c(4,2)+6.*(dwMu0dtdw(3).*c(1,2).^2+dwMu0dtdw(2).*2.*c(1,3)).*12.*c(4,3)+4.*dwMu0dtdw(2).*c(1,2).*36.*c(4,4)+dwMu0dtdw(1).*144.*c(4,5));
b(5,5)=1/576*((dwMu0dtdw(8).*c(1,2).^4+6.*dwMu0dtdw(7).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(6).*2.*c(1,3).^2+4.*dwMu0dtdw(6).*c(1,2).*6.*c(1,4)+dwMu0dtdw(5).*24.*c(1,5)).*c(2,1).^4+16.*(dwMu0dtdw(7).*c(1,2).^3+3.*dwMu0dtdw(6).*c(1,2).*2.*c(1,3)+dwMu0dtdw(5).*6.*c(1,4)).*c(2,1).^3.*c(2,2)+6.*(dwMu0dtdw(6).*c(1,2).^2+dwMu0dtdw(5).*2.*c(1,3)).*(12.*c(2,1).^2.*c(2,2).^2+4.*c(2,1).^3.*2.*c(2,3))+4.*dwMu0dtdw(5).*c(1,2).*(24.*c(2,1).*c(2,2).^3+36.*c(2,1).^2.*c(2,2).*2.*c(2,3)+4.*c(2,1).^3.*6.*c(2,4))+dwMu0dtdw(4).*(24.*c(2,2).^4+144.*c(2,1).*c(2,2).^2.*2.*c(2,3)+36.*c(2,1).^2.*2.*c(2,3).^2+48.*c(2,1).^2.*c(2,2).*6.*c(2,4)+4.*c(2,1).^3.*24.*c(2,5))+6.*(dwMu0dtdw(7).*c(1,2).^4+6.*dwMu0dtdw(6).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(5).*2.*c(1,3).^2+4.*dwMu0dtdw(5).*c(1,2).*6.*c(1,4)+dwMu0dtdw(4).*24.*c(1,5)).*c(2,1).^2.*2.*c(3,1)+3.*(dwMu0dtdw(6).*c(1,2).^4+6.*dwMu0dtdw(5).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(4).*2.*c(1,3).^2+4.*dwMu0dtdw(4).*c(1,2).*6.*c(1,4)+dwMu0dtdw(3).*24.*c(1,5)).*2.*c(3,1).^2+24.*(dwMu0dtdw(5).*c(1,2).^3+3.*dwMu0dtdw(4).*c(1,2).*2.*c(1,3)+dwMu0dtdw(3).*6.*c(1,4)).*2.*c(3,1).*2.*c(3,2)+24.*(dwMu0dtdw(6).*c(1,2).^3+3.*dwMu0dtdw(5).*c(1,2).*2.*c(1,3)+dwMu0dtdw(4).*6.*c(1,4)).*(2.*c(2,1).*c(2,2).*2.*c(3,1)+c(2,1).^2.*2.*c(3,2))+36.*(dwMu0dtdw(5).*c(1,2).^2+dwMu0dtdw(4).*2.*c(1,3)).*((2.*c(2,2).^2+2.*c(2,1).*2.*c(2,3)).*2.*c(3,1)+4.*c(2,1).*c(2,2).*2.*c(3,2)+c(2,1).^2.*4.*c(3,3))+18.*(dwMu0dtdw(4).*c(1,2).^2+dwMu0dtdw(3).*2.*c(1,3)).*(2.*2.*c(3,2).^2+2.*2.*c(3,1).*4.*c(3,3))+24.*dwMu0dtdw(4).*c(1,2).*((6.*c(2,2).*2.*c(2,3)+2.*c(2,1).*6.*c(2,4)).*2.*c(3,1)+3.*(2.*c(2,2).^2+2.*c(2,1).*2.*c(2,3)).*2.*c(3,2)+6.*c(2,1).*c(2,2).*4.*c(3,3)+c(2,1).^2.*12.*c(3,4))+12.*dwMu0dtdw(3).*c(1,2).*(6.*2.*c(3,2).*4.*c(3,3)+2.*2.*c(3,1).*12.*c(3,4))+6.*dwMu0dtdw(3).*((6.*2.*c(2,3).^2+8.*c(2,2).*6.*c(2,4)+2.*c(2,1).*24.*c(2,5)).*2.*c(3,1)+4.*(6.*c(2,2).*2.*c(2,3)+2.*c(2,1).*6.*c(2,4)).*2.*c(3,2)+6.*(2.*c(2,2).^2+2.*c(2,1).*2.*c(2,3)).*4.*c(3,3)+8.*c(2,1).*c(2,2).*12.*c(3,4)+c(2,1).^2.*48.*c(3,5))+3.*dwMu0dtdw(2).*(6.*4.*c(3,3).^2+8.*2.*c(3,2).*12.*c(3,4)+2.*2.*c(3,1).*48.*c(3,5))+4.*(dwMu0dtdw(6).*c(1,2).^4+6.*dwMu0dtdw(5).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(4).*2.*c(1,3).^2+4.*dwMu0dtdw(4).*c(1,2).*6.*c(1,4)+dwMu0dtdw(3).*24.*c(1,5)).*c(2,1).*6.*c(4,1)+16.*(dwMu0dtdw(5).*c(1,2).^3+3.*dwMu0dtdw(4).*c(1,2).*2.*c(1,3)+dwMu0dtdw(3).*6.*c(1,4)).*(c(2,2).*6.*c(4,1)+c(2,1).*6.*c(4,2))+24.*(dwMu0dtdw(4).*c(1,2).^2+dwMu0dtdw(3).*2.*c(1,3)).*(2.*c(2,3).*6.*c(4,1)+2.*c(2,2).*6.*c(4,2)+c(2,1).*12.*c(4,3))+16.*dwMu0dtdw(3).*c(1,2).*(6.*c(2,4).*6.*c(4,1)+3.*2.*c(2,3).*6.*c(4,2)+3.*c(2,2).*12.*c(4,3)+c(2,1).*36.*c(4,4))+4.*dwMu0dtdw(2).*(24.*c(2,5).*6.*c(4,1)+4.*6.*c(2,4).*6.*c(4,2)+6.*2.*c(2,3).*12.*c(4,3)+4.*c(2,2).*36.*c(4,4)+c(2,1).*144.*c(4,5))+(dwMu0dtdw(5).*c(1,2).^4+6.*dwMu0dtdw(4).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(3).*2.*c(1,3).^2+4.*dwMu0dtdw(3).*c(1,2).*6.*c(1,4)+dwMu0dtdw(2).*24.*c(1,5)).*24.*c(5,1)+4.*(dwMu0dtdw(4).*c(1,2).^3+3.*dwMu0dtdw(3).*c(1,2).*2.*c(1,3)+dwMu0dtdw(2).*6.*c(1,4)).*24.*c(5,2)+6.*(dwMu0dtdw(3).*c(1,2).^2+dwMu0dtdw(2).*2.*c(1,3)).*48.*c(5,3)+4.*dwMu0dtdw(2).*c(1,2).*144.*c(5,4)+dwMu0dtdw(1).*576.*c(5,5));
b(6,5)=1/2880*((dwMu0dtdw(9).*c(1,2).^4+6.*dwMu0dtdw(8).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(7).*2.*c(1,3).^2+4.*dwMu0dtdw(7).*c(1,2).*6.*c(1,4)+dwMu0dtdw(6).*24.*c(1,5)).*c(2,1).^5+20.*(dwMu0dtdw(8).*c(1,2).^3+3.*dwMu0dtdw(7).*c(1,2).*2.*c(1,3)+dwMu0dtdw(6).*6.*c(1,4)).*c(2,1).^4.*c(2,2)+6.*(dwMu0dtdw(7).*c(1,2).^2+dwMu0dtdw(6).*2.*c(1,3)).*(20.*c(2,1).^3.*c(2,2).^2+5.*c(2,1).^4.*2.*c(2,3))+4.*dwMu0dtdw(6).*c(1,2).*(60.*c(2,1).^2.*c(2,2).^3+60.*c(2,1).^3.*c(2,2).*2.*c(2,3)+5.*c(2,1).^4.*6.*c(2,4))+dwMu0dtdw(5).*(120.*c(2,1).*c(2,2).^4+360.*c(2,1).^2.*c(2,2).^2.*2.*c(2,3)+60.*c(2,1).^3.*2.*c(2,3).^2+80.*c(2,1).^3.*c(2,2).*6.*c(2,4)+5.*c(2,1).^4.*24.*c(2,5))+10.*(dwMu0dtdw(8).*c(1,2).^4+6.*dwMu0dtdw(7).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(6).*2.*c(1,3).^2+4.*dwMu0dtdw(6).*c(1,2).*6.*c(1,4)+dwMu0dtdw(5).*24.*c(1,5)).*c(2,1).^3.*2.*c(3,1)+15.*(dwMu0dtdw(7).*c(1,2).^4+6.*dwMu0dtdw(6).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(5).*2.*c(1,3).^2+4.*dwMu0dtdw(5).*c(1,2).*6.*c(1,4)+dwMu0dtdw(4).*24.*c(1,5)).*c(2,1).*2.*c(3,1).^2+40.*(dwMu0dtdw(7).*c(1,2).^3+3.*dwMu0dtdw(6).*c(1,2).*2.*c(1,3)+dwMu0dtdw(5).*6.*c(1,4)).*(3.*c(2,1).^2.*c(2,2).*2.*c(3,1)+c(2,1).^3.*2.*c(3,2))+60.*(dwMu0dtdw(6).*c(1,2).^3+3.*dwMu0dtdw(5).*c(1,2).*2.*c(1,3)+dwMu0dtdw(4).*6.*c(1,4)).*(c(2,2).*2.*c(3,1).^2+2.*c(2,1).*2.*c(3,1).*2.*c(3,2))+60.*(dwMu0dtdw(6).*c(1,2).^2+dwMu0dtdw(5).*2.*c(1,3)).*((6.*c(2,1).*c(2,2).^2+3.*c(2,1).^2.*2.*c(2,3)).*2.*c(3,1)+6.*c(2,1).^2.*c(2,2).*2.*c(3,2)+c(2,1).^3.*4.*c(3,3))+90.*(dwMu0dtdw(5).*c(1,2).^2+dwMu0dtdw(4).*2.*c(1,3)).*(2.*c(2,3).*2.*c(3,1).^2+4.*c(2,2).*2.*c(3,1).*2.*c(3,2)+c(2,1).*(2.*2.*c(3,2).^2+2.*2.*c(3,1).*4.*c(3,3)))+40.*dwMu0dtdw(5).*c(1,2).*((6.*c(2,2).^3+18.*c(2,1).*c(2,2).*2.*c(2,3)+3.*c(2,1).^2.*6.*c(2,4)).*2.*c(3,1)+3.*(6.*c(2,1).*c(2,2).^2+3.*c(2,1).^2.*2.*c(2,3)).*2.*c(3,2)+9.*c(2,1).^2.*c(2,2).*4.*c(3,3)+c(2,1).^3.*12.*c(3,4))+60.*dwMu0dtdw(4).*c(1,2).*(6.*c(2,4).*2.*c(3,1).^2+6.*2.*c(2,3).*2.*c(3,1).*2.*c(3,2)+3.*c(2,2).*(2.*2.*c(3,2).^2+2.*2.*c(3,1).*4.*c(3,3))+c(2,1).*(6.*2.*c(3,2).*4.*c(3,3)+2.*2.*c(3,1).*12.*c(3,4)))+10.*dwMu0dtdw(4).*((36.*c(2,2).^2.*2.*c(2,3)+18.*c(2,1).*2.*c(2,3).^2+24.*c(2,1).*c(2,2).*6.*c(2,4)+3.*c(2,1).^2.*24.*c(2,5)).*2.*c(3,1)+4.*(6.*c(2,2).^3+18.*c(2,1).*c(2,2).*2.*c(2,3)+3.*c(2,1).^2.*6.*c(2,4)).*2.*c(3,2)+6.*(6.*c(2,1).*c(2,2).^2+3.*c(2,1).^2.*2.*c(2,3)).*4.*c(3,3)+12.*c(2,1).^2.*c(2,2).*12.*c(3,4)+c(2,1).^3.*48.*c(3,5))+15.*dwMu0dtdw(3).*(24.*c(2,5).*2.*c(3,1).^2+8.*6.*c(2,4).*2.*c(3,1).*2.*c(3,2)+6.*2.*c(2,3).*(2.*2.*c(3,2).^2+2.*2.*c(3,1).*4.*c(3,3))+4.*c(2,2).*(6.*2.*c(3,2).*4.*c(3,3)+2.*2.*c(3,1).*12.*c(3,4))+c(2,1).*(6.*4.*c(3,3).^2+8.*2.*c(3,2).*12.*c(3,4)+2.*2.*c(3,1).*48.*c(3,5)))+10.*(dwMu0dtdw(7).*c(1,2).^4+6.*dwMu0dtdw(6).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(5).*2.*c(1,3).^2+4.*dwMu0dtdw(5).*c(1,2).*6.*c(1,4)+dwMu0dtdw(4).*24.*c(1,5)).*c(2,1).^2.*6.*c(4,1)+10.*(dwMu0dtdw(6).*c(1,2).^4+6.*dwMu0dtdw(5).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(4).*2.*c(1,3).^2+4.*dwMu0dtdw(4).*c(1,2).*6.*c(1,4)+dwMu0dtdw(3).*24.*c(1,5)).*2.*c(3,1).*6.*c(4,1)+40.*(dwMu0dtdw(6).*c(1,2).^3+3.*dwMu0dtdw(5).*c(1,2).*2.*c(1,3)+dwMu0dtdw(4).*6.*c(1,4)).*(2.*c(2,1).*c(2,2).*6.*c(4,1)+c(2,1).^2.*6.*c(4,2))+40.*(dwMu0dtdw(5).*c(1,2).^3+3.*dwMu0dtdw(4).*c(1,2).*2.*c(1,3)+dwMu0dtdw(3).*6.*c(1,4)).*(2.*c(3,2).*6.*c(4,1)+2.*c(3,1).*6.*c(4,2))+60.*(dwMu0dtdw(5).*c(1,2).^2+dwMu0dtdw(4).*2.*c(1,3)).*((2.*c(2,2).^2+2.*c(2,1).*2.*c(2,3)).*6.*c(4,1)+4.*c(2,1).*c(2,2).*6.*c(4,2)+c(2,1).^2.*12.*c(4,3))+60.*(dwMu0dtdw(4).*c(1,2).^2+dwMu0dtdw(3).*2.*c(1,3)).*(4.*c(3,3).*6.*c(4,1)+2.*2.*c(3,2).*6.*c(4,2)+2.*c(3,1).*12.*c(4,3))+40.*dwMu0dtdw(4).*c(1,2).*((6.*c(2,2).*2.*c(2,3)+2.*c(2,1).*6.*c(2,4)).*6.*c(4,1)+3.*(2.*c(2,2).^2+2.*c(2,1).*2.*c(2,3)).*6.*c(4,2)+6.*c(2,1).*c(2,2).*12.*c(4,3)+c(2,1).^2.*36.*c(4,4))+40.*dwMu0dtdw(3).*c(1,2).*(12.*c(3,4).*6.*c(4,1)+3.*4.*c(3,3).*6.*c(4,2)+3.*2.*c(3,2).*12.*c(4,3)+2.*c(3,1).*36.*c(4,4))+10.*dwMu0dtdw(3).*((6.*2.*c(2,3).^2+8.*c(2,2).*6.*c(2,4)+2.*c(2,1).*24.*c(2,5)).*6.*c(4,1)+4.*(6.*c(2,2).*2.*c(2,3)+2.*c(2,1).*6.*c(2,4)).*6.*c(4,2)+6.*(2.*c(2,2).^2+2.*c(2,1).*2.*c(2,3)).*12.*c(4,3)+8.*c(2,1).*c(2,2).*36.*c(4,4)+c(2,1).^2.*144.*c(4,5))+10.*dwMu0dtdw(2).*(48.*c(3,5).*6.*c(4,1)+4.*12.*c(3,4).*6.*c(4,2)+6.*4.*c(3,3).*12.*c(4,3)+4.*2.*c(3,2).*36.*c(4,4)+2.*c(3,1).*144.*c(4,5))+5.*(dwMu0dtdw(6).*c(1,2).^4+6.*dwMu0dtdw(5).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(4).*2.*c(1,3).^2+4.*dwMu0dtdw(4).*c(1,2).*6.*c(1,4)+dwMu0dtdw(3).*24.*c(1,5)).*c(2,1).*24.*c(5,1)+20.*(dwMu0dtdw(5).*c(1,2).^3+3.*dwMu0dtdw(4).*c(1,2).*2.*c(1,3)+dwMu0dtdw(3).*6.*c(1,4)).*(c(2,2).*24.*c(5,1)+c(2,1).*24.*c(5,2))+30.*(dwMu0dtdw(4).*c(1,2).^2+dwMu0dtdw(3).*2.*c(1,3)).*(2.*c(2,3).*24.*c(5,1)+2.*c(2,2).*24.*c(5,2)+c(2,1).*48.*c(5,3))+20.*dwMu0dtdw(3).*c(1,2).*(6.*c(2,4).*24.*c(5,1)+3.*2.*c(2,3).*24.*c(5,2)+3.*c(2,2).*48.*c(5,3)+c(2,1).*144.*c(5,4))+5.*dwMu0dtdw(2).*(24.*c(2,5).*24.*c(5,1)+4.*6.*c(2,4).*24.*c(5,2)+6.*2.*c(2,3).*48.*c(5,3)+4.*c(2,2).*144.*c(5,4)+c(2,1).*576.*c(5,5))+(dwMu0dtdw(5).*c(1,2).^4+6.*dwMu0dtdw(4).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(3).*2.*c(1,3).^2+4.*dwMu0dtdw(3).*c(1,2).*6.*c(1,4)+dwMu0dtdw(2).*24.*c(1,5)).*120.*c(6,1)+4.*(dwMu0dtdw(4).*c(1,2).^3+3.*dwMu0dtdw(3).*c(1,2).*2.*c(1,3)+dwMu0dtdw(2).*6.*c(1,4)).*120.*c(6,2)+6.*(dwMu0dtdw(3).*c(1,2).^2+dwMu0dtdw(2).*2.*c(1,3)).*240.*c(6,3)+4.*dwMu0dtdw(2).*c(1,2).*720.*c(6,4)+dwMu0dtdw(1).*2880.*c(6,5));
b(2,6)=1/120*((dwMu0dtdw(6).*c(1,2).^5+10.*dwMu0dtdw(5).*c(1,2).^3.*2.*c(1,3)+15.*dwMu0dtdw(4).*c(1,2).*2.*c(1,3).^2+10.*dwMu0dtdw(4).*c(1,2).^2.*6.*c(1,4)+10.*dwMu0dtdw(3).*2.*c(1,3).*6.*c(1,4)+5.*dwMu0dtdw(3).*c(1,2).*24.*c(1,5)+dwMu0dtdw(2).*120.*c(1,6)).*c(2,1)+5.*(dwMu0dtdw(5).*c(1,2).^4+6.*dwMu0dtdw(4).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(3).*2.*c(1,3).^2+4.*dwMu0dtdw(3).*c(1,2).*6.*c(1,4)+dwMu0dtdw(2).*24.*c(1,5)).*c(2,2)+10.*(dwMu0dtdw(4).*c(1,2).^3+3.*dwMu0dtdw(3).*c(1,2).*2.*c(1,3)+dwMu0dtdw(2).*6.*c(1,4)).*2.*c(2,3)+10.*(dwMu0dtdw(3).*c(1,2).^2+dwMu0dtdw(2).*2.*c(1,3)).*6.*c(2,4)+5.*dwMu0dtdw(2).*c(1,2).*24.*c(2,5)+dwMu0dtdw(1).*120.*c(2,6));
b(3,6)=1/240*((dwMu0dtdw(7).*c(1,2).^5+10.*dwMu0dtdw(6).*c(1,2).^3.*2.*c(1,3)+15.*dwMu0dtdw(5).*c(1,2).*2.*c(1,3).^2+10.*dwMu0dtdw(5).*c(1,2).^2.*6.*c(1,4)+10.*dwMu0dtdw(4).*2.*c(1,3).*6.*c(1,4)+5.*dwMu0dtdw(4).*c(1,2).*24.*c(1,5)+dwMu0dtdw(3).*120.*c(1,6)).*c(2,1).^2+10.*(dwMu0dtdw(6).*c(1,2).^4+6.*dwMu0dtdw(5).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(4).*2.*c(1,3).^2+4.*dwMu0dtdw(4).*c(1,2).*6.*c(1,4)+dwMu0dtdw(3).*24.*c(1,5)).*c(2,1).*c(2,2)+10.*(dwMu0dtdw(5).*c(1,2).^3+3.*dwMu0dtdw(4).*c(1,2).*2.*c(1,3)+dwMu0dtdw(3).*6.*c(1,4)).*(2.*c(2,2).^2+2.*c(2,1).*2.*c(2,3))+10.*(dwMu0dtdw(4).*c(1,2).^2+dwMu0dtdw(3).*2.*c(1,3)).*(6.*c(2,2).*2.*c(2,3)+2.*c(2,1).*6.*c(2,4))+5.*dwMu0dtdw(3).*c(1,2).*(6.*2.*c(2,3).^2+8.*c(2,2).*6.*c(2,4)+2.*c(2,1).*24.*c(2,5))+dwMu0dtdw(2).*(20.*2.*c(2,3).*6.*c(2,4)+10.*c(2,2).*24.*c(2,5)+2.*c(2,1).*120.*c(2,6))+(dwMu0dtdw(6).*c(1,2).^5+10.*dwMu0dtdw(5).*c(1,2).^3.*2.*c(1,3)+15.*dwMu0dtdw(4).*c(1,2).*2.*c(1,3).^2+10.*dwMu0dtdw(4).*c(1,2).^2.*6.*c(1,4)+10.*dwMu0dtdw(3).*2.*c(1,3).*6.*c(1,4)+5.*dwMu0dtdw(3).*c(1,2).*24.*c(1,5)+dwMu0dtdw(2).*120.*c(1,6)).*2.*c(3,1)+5.*(dwMu0dtdw(5).*c(1,2).^4+6.*dwMu0dtdw(4).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(3).*2.*c(1,3).^2+4.*dwMu0dtdw(3).*c(1,2).*6.*c(1,4)+dwMu0dtdw(2).*24.*c(1,5)).*2.*c(3,2)+10.*(dwMu0dtdw(4).*c(1,2).^3+3.*dwMu0dtdw(3).*c(1,2).*2.*c(1,3)+dwMu0dtdw(2).*6.*c(1,4)).*4.*c(3,3)+10.*(dwMu0dtdw(3).*c(1,2).^2+dwMu0dtdw(2).*2.*c(1,3)).*12.*c(3,4)+5.*dwMu0dtdw(2).*c(1,2).*48.*c(3,5)+dwMu0dtdw(1).*240.*c(3,6));
b(4,6)=1/720*((dwMu0dtdw(8).*c(1,2).^5+10.*dwMu0dtdw(7).*c(1,2).^3.*2.*c(1,3)+15.*dwMu0dtdw(6).*c(1,2).*2.*c(1,3).^2+10.*dwMu0dtdw(6).*c(1,2).^2.*6.*c(1,4)+10.*dwMu0dtdw(5).*2.*c(1,3).*6.*c(1,4)+5.*dwMu0dtdw(5).*c(1,2).*24.*c(1,5)+dwMu0dtdw(4).*120.*c(1,6)).*c(2,1).^3+15.*(dwMu0dtdw(7).*c(1,2).^4+6.*dwMu0dtdw(6).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(5).*2.*c(1,3).^2+4.*dwMu0dtdw(5).*c(1,2).*6.*c(1,4)+dwMu0dtdw(4).*24.*c(1,5)).*c(2,1).^2.*c(2,2)+10.*(dwMu0dtdw(6).*c(1,2).^3+3.*dwMu0dtdw(5).*c(1,2).*2.*c(1,3)+dwMu0dtdw(4).*6.*c(1,4)).*(6.*c(2,1).*c(2,2).^2+3.*c(2,1).^2.*2.*c(2,3))+10.*(dwMu0dtdw(5).*c(1,2).^2+dwMu0dtdw(4).*2.*c(1,3)).*(6.*c(2,2).^3+18.*c(2,1).*c(2,2).*2.*c(2,3)+3.*c(2,1).^2.*6.*c(2,4))+5.*dwMu0dtdw(4).*c(1,2).*(36.*c(2,2).^2.*2.*c(2,3)+18.*c(2,1).*2.*c(2,3).^2+24.*c(2,1).*c(2,2).*6.*c(2,4)+3.*c(2,1).^2.*24.*c(2,5))+dwMu0dtdw(3).*(90.*c(2,2).*2.*c(2,3).^2+60.*c(2,2).^2.*6.*c(2,4)+60.*c(2,1).*2.*c(2,3).*6.*c(2,4)+30.*c(2,1).*c(2,2).*24.*c(2,5)+3.*c(2,1).^2.*120.*c(2,6))+3.*(dwMu0dtdw(7).*c(1,2).^5+10.*dwMu0dtdw(6).*c(1,2).^3.*2.*c(1,3)+15.*dwMu0dtdw(5).*c(1,2).*2.*c(1,3).^2+10.*dwMu0dtdw(5).*c(1,2).^2.*6.*c(1,4)+10.*dwMu0dtdw(4).*2.*c(1,3).*6.*c(1,4)+5.*dwMu0dtdw(4).*c(1,2).*24.*c(1,5)+dwMu0dtdw(3).*120.*c(1,6)).*c(2,1).*2.*c(3,1)+15.*(dwMu0dtdw(6).*c(1,2).^4+6.*dwMu0dtdw(5).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(4).*2.*c(1,3).^2+4.*dwMu0dtdw(4).*c(1,2).*6.*c(1,4)+dwMu0dtdw(3).*24.*c(1,5)).*(c(2,2).*2.*c(3,1)+c(2,1).*2.*c(3,2))+30.*(dwMu0dtdw(5).*c(1,2).^3+3.*dwMu0dtdw(4).*c(1,2).*2.*c(1,3)+dwMu0dtdw(3).*6.*c(1,4)).*(2.*c(2,3).*2.*c(3,1)+2.*c(2,2).*2.*c(3,2)+c(2,1).*4.*c(3,3))+30.*(dwMu0dtdw(4).*c(1,2).^2+dwMu0dtdw(3).*2.*c(1,3)).*(6.*c(2,4).*2.*c(3,1)+3.*2.*c(2,3).*2.*c(3,2)+3.*c(2,2).*4.*c(3,3)+c(2,1).*12.*c(3,4))+15.*dwMu0dtdw(3).*c(1,2).*(24.*c(2,5).*2.*c(3,1)+4.*6.*c(2,4).*2.*c(3,2)+6.*2.*c(2,3).*4.*c(3,3)+4.*c(2,2).*12.*c(3,4)+c(2,1).*48.*c(3,5))+3.*dwMu0dtdw(2).*(120.*c(2,6).*2.*c(3,1)+5.*24.*c(2,5).*2.*c(3,2)+10.*6.*c(2,4).*4.*c(3,3)+10.*2.*c(2,3).*12.*c(3,4)+5.*c(2,2).*48.*c(3,5)+c(2,1).*240.*c(3,6))+(dwMu0dtdw(6).*c(1,2).^5+10.*dwMu0dtdw(5).*c(1,2).^3.*2.*c(1,3)+15.*dwMu0dtdw(4).*c(1,2).*2.*c(1,3).^2+10.*dwMu0dtdw(4).*c(1,2).^2.*6.*c(1,4)+10.*dwMu0dtdw(3).*2.*c(1,3).*6.*c(1,4)+5.*dwMu0dtdw(3).*c(1,2).*24.*c(1,5)+dwMu0dtdw(2).*120.*c(1,6)).*6.*c(4,1)+5.*(dwMu0dtdw(5).*c(1,2).^4+6.*dwMu0dtdw(4).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(3).*2.*c(1,3).^2+4.*dwMu0dtdw(3).*c(1,2).*6.*c(1,4)+dwMu0dtdw(2).*24.*c(1,5)).*6.*c(4,2)+10.*(dwMu0dtdw(4).*c(1,2).^3+3.*dwMu0dtdw(3).*c(1,2).*2.*c(1,3)+dwMu0dtdw(2).*6.*c(1,4)).*12.*c(4,3)+10.*(dwMu0dtdw(3).*c(1,2).^2+dwMu0dtdw(2).*2.*c(1,3)).*36.*c(4,4)+5.*dwMu0dtdw(2).*c(1,2).*144.*c(4,5)+dwMu0dtdw(1).*720.*c(4,6));
b(5,6)=1/2880*((dwMu0dtdw(9).*c(1,2).^5+10.*dwMu0dtdw(8).*c(1,2).^3.*2.*c(1,3)+15.*dwMu0dtdw(7).*c(1,2).*2.*c(1,3).^2+10.*dwMu0dtdw(7).*c(1,2).^2.*6.*c(1,4)+10.*dwMu0dtdw(6).*2.*c(1,3).*6.*c(1,4)+5.*dwMu0dtdw(6).*c(1,2).*24.*c(1,5)+dwMu0dtdw(5).*120.*c(1,6)).*c(2,1).^4+20.*(dwMu0dtdw(8).*c(1,2).^4+6.*dwMu0dtdw(7).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(6).*2.*c(1,3).^2+4.*dwMu0dtdw(6).*c(1,2).*6.*c(1,4)+dwMu0dtdw(5).*24.*c(1,5)).*c(2,1).^3.*c(2,2)+10.*(dwMu0dtdw(7).*c(1,2).^3+3.*dwMu0dtdw(6).*c(1,2).*2.*c(1,3)+dwMu0dtdw(5).*6.*c(1,4)).*(12.*c(2,1).^2.*c(2,2).^2+4.*c(2,1).^3.*2.*c(2,3))+10.*(dwMu0dtdw(6).*c(1,2).^2+dwMu0dtdw(5).*2.*c(1,3)).*(24.*c(2,1).*c(2,2).^3+36.*c(2,1).^2.*c(2,2).*2.*c(2,3)+4.*c(2,1).^3.*6.*c(2,4))+5.*dwMu0dtdw(5).*c(1,2).*(24.*c(2,2).^4+144.*c(2,1).*c(2,2).^2.*2.*c(2,3)+36.*c(2,1).^2.*2.*c(2,3).^2+48.*c(2,1).^2.*c(2,2).*6.*c(2,4)+4.*c(2,1).^3.*24.*c(2,5))+dwMu0dtdw(4).*(240.*c(2,2).^3.*2.*c(2,3)+360.*c(2,1).*c(2,2).*2.*c(2,3).^2+240.*c(2,1).*c(2,2).^2.*6.*c(2,4)+120.*c(2,1).^2.*2.*c(2,3).*6.*c(2,4)+60.*c(2,1).^2.*c(2,2).*24.*c(2,5)+4.*c(2,1).^3.*120.*c(2,6))+6.*(dwMu0dtdw(8).*c(1,2).^5+10.*dwMu0dtdw(7).*c(1,2).^3.*2.*c(1,3)+15.*dwMu0dtdw(6).*c(1,2).*2.*c(1,3).^2+10.*dwMu0dtdw(6).*c(1,2).^2.*6.*c(1,4)+10.*dwMu0dtdw(5).*2.*c(1,3).*6.*c(1,4)+5.*dwMu0dtdw(5).*c(1,2).*24.*c(1,5)+dwMu0dtdw(4).*120.*c(1,6)).*c(2,1).^2.*2.*c(3,1)+3.*(dwMu0dtdw(7).*c(1,2).^5+10.*dwMu0dtdw(6).*c(1,2).^3.*2.*c(1,3)+15.*dwMu0dtdw(5).*c(1,2).*2.*c(1,3).^2+10.*dwMu0dtdw(5).*c(1,2).^2.*6.*c(1,4)+10.*dwMu0dtdw(4).*2.*c(1,3).*6.*c(1,4)+5.*dwMu0dtdw(4).*c(1,2).*24.*c(1,5)+dwMu0dtdw(3).*120.*c(1,6)).*2.*c(3,1).^2+30.*(dwMu0dtdw(6).*c(1,2).^4+6.*dwMu0dtdw(5).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(4).*2.*c(1,3).^2+4.*dwMu0dtdw(4).*c(1,2).*6.*c(1,4)+dwMu0dtdw(3).*24.*c(1,5)).*2.*c(3,1).*2.*c(3,2)+30.*(dwMu0dtdw(7).*c(1,2).^4+6.*dwMu0dtdw(6).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(5).*2.*c(1,3).^2+4.*dwMu0dtdw(5).*c(1,2).*6.*c(1,4)+dwMu0dtdw(4).*24.*c(1,5)).*(2.*c(2,1).*c(2,2).*2.*c(3,1)+c(2,1).^2.*2.*c(3,2))+60.*(dwMu0dtdw(6).*c(1,2).^3+3.*dwMu0dtdw(5).*c(1,2).*2.*c(1,3)+dwMu0dtdw(4).*6.*c(1,4)).*((2.*c(2,2).^2+2.*c(2,1).*2.*c(2,3)).*2.*c(3,1)+4.*c(2,1).*c(2,2).*2.*c(3,2)+c(2,1).^2.*4.*c(3,3))+30.*(dwMu0dtdw(5).*c(1,2).^3+3.*dwMu0dtdw(4).*c(1,2).*2.*c(1,3)+dwMu0dtdw(3).*6.*c(1,4)).*(2.*2.*c(3,2).^2+2.*2.*c(3,1).*4.*c(3,3))+60.*(dwMu0dtdw(5).*c(1,2).^2+dwMu0dtdw(4).*2.*c(1,3)).*((6.*c(2,2).*2.*c(2,3)+2.*c(2,1).*6.*c(2,4)).*2.*c(3,1)+3.*(2.*c(2,2).^2+2.*c(2,1).*2.*c(2,3)).*2.*c(3,2)+6.*c(2,1).*c(2,2).*4.*c(3,3)+c(2,1).^2.*12.*c(3,4))+30.*(dwMu0dtdw(4).*c(1,2).^2+dwMu0dtdw(3).*2.*c(1,3)).*(6.*2.*c(3,2).*4.*c(3,3)+2.*2.*c(3,1).*12.*c(3,4))+30.*dwMu0dtdw(4).*c(1,2).*((6.*2.*c(2,3).^2+8.*c(2,2).*6.*c(2,4)+2.*c(2,1).*24.*c(2,5)).*2.*c(3,1)+4.*(6.*c(2,2).*2.*c(2,3)+2.*c(2,1).*6.*c(2,4)).*2.*c(3,2)+6.*(2.*c(2,2).^2+2.*c(2,1).*2.*c(2,3)).*4.*c(3,3)+8.*c(2,1).*c(2,2).*12.*c(3,4)+c(2,1).^2.*48.*c(3,5))+15.*dwMu0dtdw(3).*c(1,2).*(6.*4.*c(3,3).^2+8.*2.*c(3,2).*12.*c(3,4)+2.*2.*c(3,1).*48.*c(3,5))+6.*dwMu0dtdw(3).*((20.*2.*c(2,3).*6.*c(2,4)+10.*c(2,2).*24.*c(2,5)+2.*c(2,1).*120.*c(2,6)).*2.*c(3,1)+5.*(6.*2.*c(2,3).^2+8.*c(2,2).*6.*c(2,4)+2.*c(2,1).*24.*c(2,5)).*2.*c(3,2)+10.*(6.*c(2,2).*2.*c(2,3)+2.*c(2,1).*6.*c(2,4)).*4.*c(3,3)+10.*(2.*c(2,2).^2+2.*c(2,1).*2.*c(2,3)).*12.*c(3,4)+10.*c(2,1).*c(2,2).*48.*c(3,5)+c(2,1).^2.*240.*c(3,6))+3.*dwMu0dtdw(2).*(20.*4.*c(3,3).*12.*c(3,4)+10.*2.*c(3,2).*48.*c(3,5)+2.*2.*c(3,1).*240.*c(3,6))+4.*(dwMu0dtdw(7).*c(1,2).^5+10.*dwMu0dtdw(6).*c(1,2).^3.*2.*c(1,3)+15.*dwMu0dtdw(5).*c(1,2).*2.*c(1,3).^2+10.*dwMu0dtdw(5).*c(1,2).^2.*6.*c(1,4)+10.*dwMu0dtdw(4).*2.*c(1,3).*6.*c(1,4)+5.*dwMu0dtdw(4).*c(1,2).*24.*c(1,5)+dwMu0dtdw(3).*120.*c(1,6)).*c(2,1).*6.*c(4,1)+20.*(dwMu0dtdw(6).*c(1,2).^4+6.*dwMu0dtdw(5).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(4).*2.*c(1,3).^2+4.*dwMu0dtdw(4).*c(1,2).*6.*c(1,4)+dwMu0dtdw(3).*24.*c(1,5)).*(c(2,2).*6.*c(4,1)+c(2,1).*6.*c(4,2))+40.*(dwMu0dtdw(5).*c(1,2).^3+3.*dwMu0dtdw(4).*c(1,2).*2.*c(1,3)+dwMu0dtdw(3).*6.*c(1,4)).*(2.*c(2,3).*6.*c(4,1)+2.*c(2,2).*6.*c(4,2)+c(2,1).*12.*c(4,3))+40.*(dwMu0dtdw(4).*c(1,2).^2+dwMu0dtdw(3).*2.*c(1,3)).*(6.*c(2,4).*6.*c(4,1)+3.*2.*c(2,3).*6.*c(4,2)+3.*c(2,2).*12.*c(4,3)+c(2,1).*36.*c(4,4))+20.*dwMu0dtdw(3).*c(1,2).*(24.*c(2,5).*6.*c(4,1)+4.*6.*c(2,4).*6.*c(4,2)+6.*2.*c(2,3).*12.*c(4,3)+4.*c(2,2).*36.*c(4,4)+c(2,1).*144.*c(4,5))+4.*dwMu0dtdw(2).*(120.*c(2,6).*6.*c(4,1)+5.*24.*c(2,5).*6.*c(4,2)+10.*6.*c(2,4).*12.*c(4,3)+10.*2.*c(2,3).*36.*c(4,4)+5.*c(2,2).*144.*c(4,5)+c(2,1).*720.*c(4,6))+(dwMu0dtdw(6).*c(1,2).^5+10.*dwMu0dtdw(5).*c(1,2).^3.*2.*c(1,3)+15.*dwMu0dtdw(4).*c(1,2).*2.*c(1,3).^2+10.*dwMu0dtdw(4).*c(1,2).^2.*6.*c(1,4)+10.*dwMu0dtdw(3).*2.*c(1,3).*6.*c(1,4)+5.*dwMu0dtdw(3).*c(1,2).*24.*c(1,5)+dwMu0dtdw(2).*120.*c(1,6)).*24.*c(5,1)+5.*(dwMu0dtdw(5).*c(1,2).^4+6.*dwMu0dtdw(4).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(3).*2.*c(1,3).^2+4.*dwMu0dtdw(3).*c(1,2).*6.*c(1,4)+dwMu0dtdw(2).*24.*c(1,5)).*24.*c(5,2)+10.*(dwMu0dtdw(4).*c(1,2).^3+3.*dwMu0dtdw(3).*c(1,2).*2.*c(1,3)+dwMu0dtdw(2).*6.*c(1,4)).*48.*c(5,3)+10.*(dwMu0dtdw(3).*c(1,2).^2+dwMu0dtdw(2).*2.*c(1,3)).*144.*c(5,4)+5.*dwMu0dtdw(2).*c(1,2).*576.*c(5,5)+dwMu0dtdw(1).*2880.*c(5,6));
b(6,6)=1/14400*((dwMu0dtdw(10).*c(1,2).^5+10.*dwMu0dtdw(9).*c(1,2).^3.*2.*c(1,3)+15.*dwMu0dtdw(8).*c(1,2).*2.*c(1,3).^2+10.*dwMu0dtdw(8).*c(1,2).^2.*6.*c(1,4)+10.*dwMu0dtdw(7).*2.*c(1,3).*6.*c(1,4)+5.*dwMu0dtdw(7).*c(1,2).*24.*c(1,5)+dwMu0dtdw(6).*120.*c(1,6)).*c(2,1).^5+25.*(dwMu0dtdw(9).*c(1,2).^4+6.*dwMu0dtdw(8).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(7).*2.*c(1,3).^2+4.*dwMu0dtdw(7).*c(1,2).*6.*c(1,4)+dwMu0dtdw(6).*24.*c(1,5)).*c(2,1).^4.*c(2,2)+10.*(dwMu0dtdw(8).*c(1,2).^3+3.*dwMu0dtdw(7).*c(1,2).*2.*c(1,3)+dwMu0dtdw(6).*6.*c(1,4)).*(20.*c(2,1).^3.*c(2,2).^2+5.*c(2,1).^4.*2.*c(2,3))+10.*(dwMu0dtdw(7).*c(1,2).^2+dwMu0dtdw(6).*2.*c(1,3)).*(60.*c(2,1).^2.*c(2,2).^3+60.*c(2,1).^3.*c(2,2).*2.*c(2,3)+5.*c(2,1).^4.*6.*c(2,4))+5.*dwMu0dtdw(6).*c(1,2).*(120.*c(2,1).*c(2,2).^4+360.*c(2,1).^2.*c(2,2).^2.*2.*c(2,3)+60.*c(2,1).^3.*2.*c(2,3).^2+80.*c(2,1).^3.*c(2,2).*6.*c(2,4)+5.*c(2,1).^4.*24.*c(2,5))+dwMu0dtdw(5).*(120.*c(2,2).^5+1200.*c(2,1).*c(2,2).^3.*2.*c(2,3)+900.*c(2,1).^2.*c(2,2).*2.*c(2,3).^2+600.*c(2,1).^2.*c(2,2).^2.*6.*c(2,4)+200.*c(2,1).^3.*2.*c(2,3).*6.*c(2,4)+100.*c(2,1).^3.*c(2,2).*24.*c(2,5)+5.*c(2,1).^4.*120.*c(2,6))+10.*(dwMu0dtdw(9).*c(1,2).^5+10.*dwMu0dtdw(8).*c(1,2).^3.*2.*c(1,3)+15.*dwMu0dtdw(7).*c(1,2).*2.*c(1,3).^2+10.*dwMu0dtdw(7).*c(1,2).^2.*6.*c(1,4)+10.*dwMu0dtdw(6).*2.*c(1,3).*6.*c(1,4)+5.*dwMu0dtdw(6).*c(1,2).*24.*c(1,5)+dwMu0dtdw(5).*120.*c(1,6)).*c(2,1).^3.*2.*c(3,1)+15.*(dwMu0dtdw(8).*c(1,2).^5+10.*dwMu0dtdw(7).*c(1,2).^3.*2.*c(1,3)+15.*dwMu0dtdw(6).*c(1,2).*2.*c(1,3).^2+10.*dwMu0dtdw(6).*c(1,2).^2.*6.*c(1,4)+10.*dwMu0dtdw(5).*2.*c(1,3).*6.*c(1,4)+5.*dwMu0dtdw(5).*c(1,2).*24.*c(1,5)+dwMu0dtdw(4).*120.*c(1,6)).*c(2,1).*2.*c(3,1).^2+50.*(dwMu0dtdw(8).*c(1,2).^4+6.*dwMu0dtdw(7).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(6).*2.*c(1,3).^2+4.*dwMu0dtdw(6).*c(1,2).*6.*c(1,4)+dwMu0dtdw(5).*24.*c(1,5)).*(3.*c(2,1).^2.*c(2,2).*2.*c(3,1)+c(2,1).^3.*2.*c(3,2))+75.*(dwMu0dtdw(7).*c(1,2).^4+6.*dwMu0dtdw(6).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(5).*2.*c(1,3).^2+4.*dwMu0dtdw(5).*c(1,2).*6.*c(1,4)+dwMu0dtdw(4).*24.*c(1,5)).*(c(2,2).*2.*c(3,1).^2+2.*c(2,1).*2.*c(3,1).*2.*c(3,2))+100.*(dwMu0dtdw(7).*c(1,2).^3+3.*dwMu0dtdw(6).*c(1,2).*2.*c(1,3)+dwMu0dtdw(5).*6.*c(1,4)).*((6.*c(2,1).*c(2,2).^2+3.*c(2,1).^2.*2.*c(2,3)).*2.*c(3,1)+6.*c(2,1).^2.*c(2,2).*2.*c(3,2)+c(2,1).^3.*4.*c(3,3))+150.*(dwMu0dtdw(6).*c(1,2).^3+3.*dwMu0dtdw(5).*c(1,2).*2.*c(1,3)+dwMu0dtdw(4).*6.*c(1,4)).*(2.*c(2,3).*2.*c(3,1).^2+4.*c(2,2).*2.*c(3,1).*2.*c(3,2)+c(2,1).*(2.*2.*c(3,2).^2+2.*2.*c(3,1).*4.*c(3,3)))+100.*(dwMu0dtdw(6).*c(1,2).^2+dwMu0dtdw(5).*2.*c(1,3)).*((6.*c(2,2).^3+18.*c(2,1).*c(2,2).*2.*c(2,3)+3.*c(2,1).^2.*6.*c(2,4)).*2.*c(3,1)+3.*(6.*c(2,1).*c(2,2).^2+3.*c(2,1).^2.*2.*c(2,3)).*2.*c(3,2)+9.*c(2,1).^2.*c(2,2).*4.*c(3,3)+c(2,1).^3.*12.*c(3,4))+150.*(dwMu0dtdw(5).*c(1,2).^2+dwMu0dtdw(4).*2.*c(1,3)).*(6.*c(2,4).*2.*c(3,1).^2+6.*2.*c(2,3).*2.*c(3,1).*2.*c(3,2)+3.*c(2,2).*(2.*2.*c(3,2).^2+2.*2.*c(3,1).*4.*c(3,3))+c(2,1).*(6.*2.*c(3,2).*4.*c(3,3)+2.*2.*c(3,1).*12.*c(3,4)))+50.*dwMu0dtdw(5).*c(1,2).*((36.*c(2,2).^2.*2.*c(2,3)+18.*c(2,1).*2.*c(2,3).^2+24.*c(2,1).*c(2,2).*6.*c(2,4)+3.*c(2,1).^2.*24.*c(2,5)).*2.*c(3,1)+4.*(6.*c(2,2).^3+18.*c(2,1).*c(2,2).*2.*c(2,3)+3.*c(2,1).^2.*6.*c(2,4)).*2.*c(3,2)+6.*(6.*c(2,1).*c(2,2).^2+3.*c(2,1).^2.*2.*c(2,3)).*4.*c(3,3)+12.*c(2,1).^2.*c(2,2).*12.*c(3,4)+c(2,1).^3.*48.*c(3,5))+75.*dwMu0dtdw(4).*c(1,2).*(24.*c(2,5).*2.*c(3,1).^2+8.*6.*c(2,4).*2.*c(3,1).*2.*c(3,2)+6.*2.*c(2,3).*(2.*2.*c(3,2).^2+2.*2.*c(3,1).*4.*c(3,3))+4.*c(2,2).*(6.*2.*c(3,2).*4.*c(3,3)+2.*2.*c(3,1).*12.*c(3,4))+c(2,1).*(6.*4.*c(3,3).^2+8.*2.*c(3,2).*12.*c(3,4)+2.*2.*c(3,1).*48.*c(3,5)))+10.*dwMu0dtdw(4).*((90.*c(2,2).*2.*c(2,3).^2+60.*c(2,2).^2.*6.*c(2,4)+60.*c(2,1).*2.*c(2,3).*6.*c(2,4)+30.*c(2,1).*c(2,2).*24.*c(2,5)+3.*c(2,1).^2.*120.*c(2,6)).*2.*c(3,1)+5.*(36.*c(2,2).^2.*2.*c(2,3)+18.*c(2,1).*2.*c(2,3).^2+24.*c(2,1).*c(2,2).*6.*c(2,4)+3.*c(2,1).^2.*24.*c(2,5)).*2.*c(3,2)+10.*(6.*c(2,2).^3+18.*c(2,1).*c(2,2).*2.*c(2,3)+3.*c(2,1).^2.*6.*c(2,4)).*4.*c(3,3)+10.*(6.*c(2,1).*c(2,2).^2+3.*c(2,1).^2.*2.*c(2,3)).*12.*c(3,4)+15.*c(2,1).^2.*c(2,2).*48.*c(3,5)+c(2,1).^3.*240.*c(3,6))+15.*dwMu0dtdw(3).*(120.*c(2,6).*2.*c(3,1).^2+10.*24.*c(2,5).*2.*c(3,1).*2.*c(3,2)+10.*6.*c(2,4).*(2.*2.*c(3,2).^2+2.*2.*c(3,1).*4.*c(3,3))+10.*2.*c(2,3).*(6.*2.*c(3,2).*4.*c(3,3)+2.*2.*c(3,1).*12.*c(3,4))+5.*c(2,2).*(6.*4.*c(3,3).^2+8.*2.*c(3,2).*12.*c(3,4)+2.*2.*c(3,1).*48.*c(3,5))+c(2,1).*(20.*4.*c(3,3).*12.*c(3,4)+10.*2.*c(3,2).*48.*c(3,5)+2.*2.*c(3,1).*240.*c(3,6)))+10.*(dwMu0dtdw(8).*c(1,2).^5+10.*dwMu0dtdw(7).*c(1,2).^3.*2.*c(1,3)+15.*dwMu0dtdw(6).*c(1,2).*2.*c(1,3).^2+10.*dwMu0dtdw(6).*c(1,2).^2.*6.*c(1,4)+10.*dwMu0dtdw(5).*2.*c(1,3).*6.*c(1,4)+5.*dwMu0dtdw(5).*c(1,2).*24.*c(1,5)+dwMu0dtdw(4).*120.*c(1,6)).*c(2,1).^2.*6.*c(4,1)+10.*(dwMu0dtdw(7).*c(1,2).^5+10.*dwMu0dtdw(6).*c(1,2).^3.*2.*c(1,3)+15.*dwMu0dtdw(5).*c(1,2).*2.*c(1,3).^2+10.*dwMu0dtdw(5).*c(1,2).^2.*6.*c(1,4)+10.*dwMu0dtdw(4).*2.*c(1,3).*6.*c(1,4)+5.*dwMu0dtdw(4).*c(1,2).*24.*c(1,5)+dwMu0dtdw(3).*120.*c(1,6)).*2.*c(3,1).*6.*c(4,1)+50.*(dwMu0dtdw(7).*c(1,2).^4+6.*dwMu0dtdw(6).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(5).*2.*c(1,3).^2+4.*dwMu0dtdw(5).*c(1,2).*6.*c(1,4)+dwMu0dtdw(4).*24.*c(1,5)).*(2.*c(2,1).*c(2,2).*6.*c(4,1)+c(2,1).^2.*6.*c(4,2))+50.*(dwMu0dtdw(6).*c(1,2).^4+6.*dwMu0dtdw(5).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(4).*2.*c(1,3).^2+4.*dwMu0dtdw(4).*c(1,2).*6.*c(1,4)+dwMu0dtdw(3).*24.*c(1,5)).*(2.*c(3,2).*6.*c(4,1)+2.*c(3,1).*6.*c(4,2))+100.*(dwMu0dtdw(6).*c(1,2).^3+3.*dwMu0dtdw(5).*c(1,2).*2.*c(1,3)+dwMu0dtdw(4).*6.*c(1,4)).*((2.*c(2,2).^2+2.*c(2,1).*2.*c(2,3)).*6.*c(4,1)+4.*c(2,1).*c(2,2).*6.*c(4,2)+c(2,1).^2.*12.*c(4,3))+100.*(dwMu0dtdw(5).*c(1,2).^3+3.*dwMu0dtdw(4).*c(1,2).*2.*c(1,3)+dwMu0dtdw(3).*6.*c(1,4)).*(4.*c(3,3).*6.*c(4,1)+2.*2.*c(3,2).*6.*c(4,2)+2.*c(3,1).*12.*c(4,3))+100.*(dwMu0dtdw(5).*c(1,2).^2+dwMu0dtdw(4).*2.*c(1,3)).*((6.*c(2,2).*2.*c(2,3)+2.*c(2,1).*6.*c(2,4)).*6.*c(4,1)+3.*(2.*c(2,2).^2+2.*c(2,1).*2.*c(2,3)).*6.*c(4,2)+6.*c(2,1).*c(2,2).*12.*c(4,3)+c(2,1).^2.*36.*c(4,4))+100.*(dwMu0dtdw(4).*c(1,2).^2+dwMu0dtdw(3).*2.*c(1,3)).*(12.*c(3,4).*6.*c(4,1)+3.*4.*c(3,3).*6.*c(4,2)+3.*2.*c(3,2).*12.*c(4,3)+2.*c(3,1).*36.*c(4,4))+50.*dwMu0dtdw(4).*c(1,2).*((6.*2.*c(2,3).^2+8.*c(2,2).*6.*c(2,4)+2.*c(2,1).*24.*c(2,5)).*6.*c(4,1)+4.*(6.*c(2,2).*2.*c(2,3)+2.*c(2,1).*6.*c(2,4)).*6.*c(4,2)+6.*(2.*c(2,2).^2+2.*c(2,1).*2.*c(2,3)).*12.*c(4,3)+8.*c(2,1).*c(2,2).*36.*c(4,4)+c(2,1).^2.*144.*c(4,5))+50.*dwMu0dtdw(3).*c(1,2).*(48.*c(3,5).*6.*c(4,1)+4.*12.*c(3,4).*6.*c(4,2)+6.*4.*c(3,3).*12.*c(4,3)+4.*2.*c(3,2).*36.*c(4,4)+2.*c(3,1).*144.*c(4,5))+10.*dwMu0dtdw(3).*((20.*2.*c(2,3).*6.*c(2,4)+10.*c(2,2).*24.*c(2,5)+2.*c(2,1).*120.*c(2,6)).*6.*c(4,1)+5.*(6.*2.*c(2,3).^2+8.*c(2,2).*6.*c(2,4)+2.*c(2,1).*24.*c(2,5)).*6.*c(4,2)+10.*(6.*c(2,2).*2.*c(2,3)+2.*c(2,1).*6.*c(2,4)).*12.*c(4,3)+10.*(2.*c(2,2).^2+2.*c(2,1).*2.*c(2,3)).*36.*c(4,4)+10.*c(2,1).*c(2,2).*144.*c(4,5)+c(2,1).^2.*720.*c(4,6))+10.*dwMu0dtdw(2).*(240.*c(3,6).*6.*c(4,1)+5.*48.*c(3,5).*6.*c(4,2)+10.*12.*c(3,4).*12.*c(4,3)+10.*4.*c(3,3).*36.*c(4,4)+5.*2.*c(3,2).*144.*c(4,5)+2.*c(3,1).*720.*c(4,6))+5.*(dwMu0dtdw(7).*c(1,2).^5+10.*dwMu0dtdw(6).*c(1,2).^3.*2.*c(1,3)+15.*dwMu0dtdw(5).*c(1,2).*2.*c(1,3).^2+10.*dwMu0dtdw(5).*c(1,2).^2.*6.*c(1,4)+10.*dwMu0dtdw(4).*2.*c(1,3).*6.*c(1,4)+5.*dwMu0dtdw(4).*c(1,2).*24.*c(1,5)+dwMu0dtdw(3).*120.*c(1,6)).*c(2,1).*24.*c(5,1)+25.*(dwMu0dtdw(6).*c(1,2).^4+6.*dwMu0dtdw(5).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(4).*2.*c(1,3).^2+4.*dwMu0dtdw(4).*c(1,2).*6.*c(1,4)+dwMu0dtdw(3).*24.*c(1,5)).*(c(2,2).*24.*c(5,1)+c(2,1).*24.*c(5,2))+50.*(dwMu0dtdw(5).*c(1,2).^3+3.*dwMu0dtdw(4).*c(1,2).*2.*c(1,3)+dwMu0dtdw(3).*6.*c(1,4)).*(2.*c(2,3).*24.*c(5,1)+2.*c(2,2).*24.*c(5,2)+c(2,1).*48.*c(5,3))+50.*(dwMu0dtdw(4).*c(1,2).^2+dwMu0dtdw(3).*2.*c(1,3)).*(6.*c(2,4).*24.*c(5,1)+3.*2.*c(2,3).*24.*c(5,2)+3.*c(2,2).*48.*c(5,3)+c(2,1).*144.*c(5,4))+25.*dwMu0dtdw(3).*c(1,2).*(24.*c(2,5).*24.*c(5,1)+4.*6.*c(2,4).*24.*c(5,2)+6.*2.*c(2,3).*48.*c(5,3)+4.*c(2,2).*144.*c(5,4)+c(2,1).*576.*c(5,5))+5.*dwMu0dtdw(2).*(120.*c(2,6).*24.*c(5,1)+5.*24.*c(2,5).*24.*c(5,2)+10.*6.*c(2,4).*48.*c(5,3)+10.*2.*c(2,3).*144.*c(5,4)+5.*c(2,2).*576.*c(5,5)+c(2,1).*2880.*c(5,6))+(dwMu0dtdw(6).*c(1,2).^5+10.*dwMu0dtdw(5).*c(1,2).^3.*2.*c(1,3)+15.*dwMu0dtdw(4).*c(1,2).*2.*c(1,3).^2+10.*dwMu0dtdw(4).*c(1,2).^2.*6.*c(1,4)+10.*dwMu0dtdw(3).*2.*c(1,3).*6.*c(1,4)+5.*dwMu0dtdw(3).*c(1,2).*24.*c(1,5)+dwMu0dtdw(2).*120.*c(1,6)).*120.*c(6,1)+5.*(dwMu0dtdw(5).*c(1,2).^4+6.*dwMu0dtdw(4).*c(1,2).^2.*2.*c(1,3)+3.*dwMu0dtdw(3).*2.*c(1,3).^2+4.*dwMu0dtdw(3).*c(1,2).*6.*c(1,4)+dwMu0dtdw(2).*24.*c(1,5)).*120.*c(6,2)+10.*(dwMu0dtdw(4).*c(1,2).^3+3.*dwMu0dtdw(3).*c(1,2).*2.*c(1,3)+dwMu0dtdw(2).*6.*c(1,4)).*240.*c(6,3)+10.*(dwMu0dtdw(3).*c(1,2).^2+dwMu0dtdw(2).*2.*c(1,3)).*720.*c(6,4)+5.*dwMu0dtdw(2).*c(1,2).*2880.*c(6,5)+dwMu0dtdw(1).*14400.*c(6,6));
% bv00=b(1,1);
% bv0(1:5)=b(1,2:6);
% bv10=b(2,1);
% bv1(1:5)=b(2,2:6);
% bv20=b(3,1);
% bv2(1:5)=b(3,2:6);
% bv30=b(4,1);
% bv3(1:5)=b(4,2:6);
% bv40=b(5,1);
% bv4(1:5)=b(5,2:6);
% bv50=b(6,1);
% bv5(1:5)=b(6,2:6);
%
end
.
.
.
.
.
.
function [c] =SeriesProduct2D(a,b)
c(1,1)=a(1,1).* b(1,1);
c(2,1)=(a(2,1).* b(1,1)+a(1,1).* b(2,1));% Z
c(3,1)=(a(3,1).* b(1,1)+a(2,1).* b(2,1)+a(1,1).* b(3,1));% Z^2
c(4,1)=(a(4,1).* b(1,1)+a(3,1).* b(2,1)+a(2,1).* b(3,1)+a(1,1).* b(4,1));% Z^3
c(5,1)=(a(5,1).* b(1,1)+a(4,1).* b(2,1)+a(3,1).* b(3,1)+a(2,1).* b(4,1)+a(1,1).* b(5,1));% Z^4
c(6,1)=(a(6,1).* b(1,1)+a(5,1).* b(2,1)+a(4,1).* b(3,1)+a(3,1).* b(4,1)+a(2,1).* b(5,1)+a(1,1).* b(6,1));% Z^5
%+W
c(1,2)=(a(1,2).* b(1,1)+a(1,1).* b(1,2));
c(2,2)=(a(2,2).* b(1,1)+a(2,1).* b(1,2)+a(1,2) .*b(2,1)+a(1,1).* b(2,2));% Z
c(3,2)=(a(3,2).* b(1,1)+a(3,1).* b(1,2)+a(2,2).* b(2,1)+a(2,1).* b(2,2)+a(1,2).* b(3,1)+a(1,1).* b(3,2));% Z^2
c(4,2)=(a(4,2).* b(1,1)+a(4,1).* b(1,2)+a(3,2).* b(2,1)+a(3,1).* b(2,2)+a(2,2).* b(3,1)+a(2,1).* b(3,2)+a(1,2).* b(4,1)+a(1,1).* b(4,2));% Z^3
c(5,2)=(a(5,2).* b(1,1)+a(5,1).* b(1,2)+a(4,2).* b(2,1)+a(4,1).* b(2,2)+a(3,2).* b(3,1)+a(3,1).* b(3,2)+a(2,2).* b(4,1)+a(2,1).* b(4,2)+a(1,2).* b(5,1)+a(1,1).* b(5,2));% Z^4
c(6,2)=(a(6,2).* b(1,1)+a(6,1).* b(1,2)+a(5,2).* b(2,1)+a(5,1).* b(2,2)+a(4,2).* b(3,1)+a(4,1).* b(3,2)+a(3,2).* b(4,1)+a(3,1).* b(4,2)+a(2,2).* b(5,1)+a(2,1).* b(5,2)+a(1,2).* b(6,1)+a(1,1).* b(6,2));% Z^5
;%+W^2
c(1,3)=(a(1,3).* b(1,1)+a(1,2).* b(1,2)+a(1,1).* b(1,3));;%
c(2,3)=(a(2,3).* b(1,1)+a(2,2).* b(1,2)+a(2,1).* b(1,3)+a(1,3).* b(2,1)+a(1,2).* b(2,2)+a(1,1).* b(2,3));% Z
c(3,3)=(a(3,3).* b(1,1)+a(3,2).* b(1,2)+a(3,1).* b(1,3)+a(2,3).* b(2,1)+a(2,2).* b(2,2)+a(2,1).* b(2,3)+a(1,3).* b(3,1)+a(1,2).* b(3,2)+a(1,1).* b(3,3));% Z^2
c(4,3)=(a(4,3).* b(1,1)+a(4,2).* b(1,2)+a(4,1).* b(1,3)+a(3,3).* b(2,1)+a(3,2).* b(2,2)+a(3,1).* b(2,3)+a(2,3).* b(3,1)+a(2,2).* b(3,2)+a(2,1).* b(3,3)+a(1,3).* b(4,1)+a(1,2).* b(4,2)+a(1,1).* b(4,3));% Z^3
c(5,3)=(a(5,3).* b(1,1)+a(5,2).* b(1,2)+a(5,1).* b(1,3)+a(4,3).* b(2,1)+a(4,2).* b(2,2)+a(4,1).* b(2,3)+a(3,3).* b(3,1)+a(3,2).* b(3,2)+a(3,1).* b(3,3)+a(2,3).* b(4,1)+a(2,2).* b(4,2)+a(2,1).* b(4,3)+a(1,3).* b(5,1)+a(1,2).* b(5,2)+a(1,1).* b(5,3));% Z^4
c(6,3)=(a(6,3).* b(1,1)+a(6,2).* b(1,2)+a(6,1).* b(1,3)+a(5,3).* b(2,1)+a(5,2).* b(2,2)+a(5,1).* b(2,3)+a(4,3).* b(3,1)+a(4,2).* b(3,2)+a(4,1).* b(3,3)+a(3,3).* b(4,1)+a(3,2).* b(4,2)+a(3,1).* b(4,3)+a(2,3).* b(5,1)+a(2,2).* b(5,2)+a(2,1).* b(5,3)+a(1,3).* b(6,1)+a(1,2).* b(6,2)+a(1,1).* b(6,3));% Z^5
%;%+W^3
c(1,4)=(a(1,4).* b(1,1)+a(1,3).* b(1,2)+a(1,2).* b(1,3)+a(1,1).* b(1,4));;%
c(2,4)=(a(2,4).* b(1,1)+a(2,3).* b(1,2)+a(2,2).* b(1,3)+a(2,1).* b(1,4)+a(1,4).* b(2,1)+a(1,3).* b(2,2)+a(1,2).* b(2,3)+a(1,1).* b(2,4));% Z
c(3,4)=(a(3,4).* b(1,1)+a(3,3).* b(1,2)+a(3,2).* b(1,3)+a(3,1).* b(1,4)+a(2,4).* b(2,1)+a(2,3).* b(2,2)+a(2,2).* b(2,3)+a(2,1).* b(2,4)+a(1,4).* b(3,1)+a(1,3).* b(3,2)+a(1,2).* b(3,3)+a(1,1).* b(3,4));% Z^2
c(4,4)=(a(4,4).* b(1,1)+a(4,3).* b(1,2)+a(4,2).* b(1,3)+a(4,1).* b(1,4)+a(3,4).* b(2,1)+a(3,3).* b(2,2)+a(3,2).* b(2,3)+a(3,1).* b(2,4)+a(2,4).* b(3,1)+a(2,3).* b(3,2)+a(2,2).* b(3,3)+a(2,1).* b(3,4)+a(1,4).* b(4,1)+a(1,3).* b(4,2)+a(1,2).* b(4,3)+a(1,1).* b(4,4));% Z^3
c(5,4)=(a(5,4).* b(1,1)+a(5,3).* b(1,2)+a(5,2).* b(1,3)+a(5,1).* b(1,4)+a(4,4).* b(2,1)+a(4,3).* b(2,2)+a(4,2).* b(2,3)+a(4,1).* b(2,4)+a(3,4).* b(3,1)+a(3,3).* b(3,2)+a(3,2).* b(3,3)+a(3,1).* b(3,4)+a(2,4).* b(4,1)+a(2,3).* b(4,2)+a(2,2).* b(4,3)+a(2,1).* b(4,4)+a(1,4).* b(5,1)+a(1,3).* b(5,2)+a(1,2) .*b(5,3)+a(1,1).* b(5,4));% Z^4
c(6,4)=(a(6,4).* b(1,1)+a(6,3).* b(1,2)+a(6,2).* b(1,3)+a(6,1).* b(1,4)+a(5,4).* b(2,1)+a(5,3).* b(2,2)+a(5,2).* b(2,3)+a(5,1).* b(2,4)+a(4,4).* b(3,1)+a(4,3).* b(3,2)+a(4,2).* b(3,3)+a(4,1).* b(3,4)+a(3,4).* b(4,1)+a(3,3).* b(4,2)+a(3,2).* b(4,3)+a(3,1).* b(4,4)+a(2,4).* b(5,1)+a(2,3).* b(5,2)+a(2,2).* b(5,3)+a(2,1).* b(5,4)+a(1,4).* b(6,1)+a(1,3).* b(6,2)+a(1,2).* b(6,3)+a(1,1).* b(6,4));% Z^5
%;%+W^4
c(1,5)=(a(1,5).* b(1,1)+a(1,4).* b(1,2)+a(1,3).* b(1,3)+a(1,2).* b(1,4)+a(1,1).* b(1,5));
c(2,5)=(a(2,5).* b(1,1)+a(2,4).* b(1,2)+a(2,3).* b(1,3)+a(2,2).* b(1,4)+a(2,1).* b(1,5)+a(1,5).* b(2,1)+a(1,4).* b(2,2)+a(1,3).* b(2,3)+a(1,2).* b(2,4)+a(1,1).* b(2,5));% Z
c(3,5)=(a(3,5).* b(1,1)+a(3,4).* b(1,2)+a(3,3).* b(1,3)+a(3,2).* b(1,4)+a(3,1).* b(1,5)+a(2,5).* b(2,1)+a(2,4).* b(2,2)+a(2,3).* b(2,3)+a(2,2).* b(2,4)+a(2,1).* b(2,5)+a(1,5).* b(3,1)+a(1,4).* b(3,2)+a(1,3).* b(3,3)+a(1,2).* b(3,4)+a(1,1).* b(3,5));% Z^2
c(4,5)=(a(4,5).* b(1,1)+a(4,4).* b(1,2)+a(4,3).* b(1,3)+a(4,2).* b(1,4)+a(4,1).* b(1,5)+a(3,5).* b(2,1)+a(3,4).* b(2,2)+a(3,3).* b(2,3)+a(3,2).* b(2,4)+a(3,1).* b(2,5)+a(2,5).* b(3,1)+a(2,4).* b(3,2)+a(2,3).* b(3,3)+a(2,2).* b(3,4)+a(2,1).* b(3,5)+a(1,5).* b(4,1)+a(1,4).* b(4,2)+a(1,3).* b(4,3)+a(1,2).* b(4,4)+a(1,1).* b(4,5));% Z^3
c(5,5)=(a(5,5).* b(1,1)+a(5,4).* b(1,2)+a(5,3).* b(1,3)+a(5,2).* b(1,4)+a(5,1).* b(1,5)+a(4,5).* b(2,1)+a(4,4).* b(2,2)+a(4,3).* b(2,3)+a(4,2).* b(2,4)+a(4,1).* b(2,5)+a(3,5).* b(3,1)+a(3,4).* b(3,2)+a(3,3).* b(3,3)+a(3,2).* b(3,4)+a(3,1).* b(3,5)+a(2,5).* b(4,1)+a(2,4).* b(4,2)+a(2,3).* b(4,3)+a(2,2).* b(4,4)+a(2,1).* b(4,5)+a(1,5).* b(5,1)+a(1,4).* b(5,2)+a(1,3).* b(5,3)+a(1,2).* b(5,4)+a(1,1).* b(5,5));% Z^4
c(6,5)=(a(6,5).* b(1,1)+a(6,4).* b(1,2)+a(6,3).* b(1,3)+a(6,2).* b(1,4)+a(6,1).* b(1,5)+a(5,5).* b(2,1)+a(5,4).* b(2,2)+a(5,3).* b(2,3)+a(5,2).* b(2,4)+a(5,1).* b(2,5)+a(4,5).* b(3,1)+a(4,4).* b(3,2)+a(4,3).* b(3,3)+a(4,2).* b(3,4)+a(4,1).* b(3,5)+a(3,5).* b(4,1)+a(3,4).* b(4,2)+a(3,3).* b(4,3)+a(3,2).* b(4,4)+a(3,1).* b(4,5)+a(2,5).* b(5,1)+a(2,4).* b(5,2)+a(2,3).* b(5,3)+a(2,2).* b(5,4)+a(2,1).* b(5,5)+a(1,5).* b(6,1)+a(1,4).* b(6,2)+a(1,3).* b(6,3)+a(1,2).* b(6,4)+a(1,1).* b(6,5));% Z^5
%;%W^5
c(1,6)=(a(1,6).* b(1,1)+a(1,5).* b(1,2)+a(1,4).* b(1,3)+a(1,3).* b(1,4)+a(1,2).* b(1,5)+a(1,1).* b(1,6));
c(2,6)=(a(2,6).* b(1,1)+a(2,5).* b(1,2)+a(2,4).* b(1,3)+a(2,3).* b(1,4)+a(2,2).* b(1,5)+a(2,1).* b(1,6)+a(1,6).* b(2,1)+a(1,5).* b(2,2)+a(1,4).* b(2,3)+a(1,3).* b(2,4)+a(1,2).* b(2,5)+a(1,1).* b(2,6));% Z
c(3,6)=(a(3,6).* b(1,1)+a(3,5).* b(1,2)+a(3,4).* b(1,3)+a(3,3).* b(1,4)+a(3,2).* b(1,5)+a(3,1).* b(1,6)+a(2,6).* b(2,1)+a(2,5).* b(2,2)+a(2,4).* b(2,3)+a(2,3).* b(2,4)+a(2,2).* b(2,5)+a(2,1).* b(2,6)+a(1,6).* b(3,1)+a(1,5).* b(3,2)+a(1,4).* b(3,3)+a(1,3).* b(3,4)+a(1,2).* b(3,5)+a(1,1).* b(3,6));% Z^2
c(4,6)=(a(4,6).* b(1,1)+a(4,5).* b(1,2)+a(4,4).* b(1,3)+a(4,3).* b(1,4)+a(4,2).* b(1,5)+a(4,1).* b(1,6)+a(3,6).* b(2,1)+a(3,5).* b(2,2)+a(3,4).* b(2,3)+a(3,3).* b(2,4)+a(3,2).* b(2,5)+a(3,1).* b(2,6)+a(2,6).* b(3,1)+a(2,5).* b(3,2)+a(2,4).* b(3,3)+a(2,3).* b(3,4)+a(2,2).* b(3,5)+a(2,1).* b(3,6)+a(1,6).* b(4,1)+a(1,5).* b(4,2)+a(1,4).* b(4,3)+a(1,3).* b(4,4)+a(1,2).* b(4,5)+a(1,1).* b(4,6));% Z^3
c(5,6)=(a(5,6).* b(1,1)+a(5,5).* b(1,2)+a(5,4).* b(1,3)+a(5,3).* b(1,4)+a(5,2).* b(1,5)+a(5,1).* b(1,6)+a(4,6).* b(2,1)+a(4,5).* b(2,2)+a(4,4).* b(2,3)+a(4,3).* b(2,4)+a(4,2).* b(2,5)+a(4,1).* b(2,6)+a(3,6).* b(3,1)+a(3,5).* b(3,2)+a(3,4).* b(3,3)+a(3,3).* b(3,4)+a(3,2).* b(3,5)+a(3,1).* b(3,6)+a(2,6).* b(4,1)+a(2,5).* b(4,2)+a(2,4).* b(4,3)+a(2,3).* b(4,4)+a(2,2).* b(4,5)+a(2,1).* b(4,6)+a(1,6).* b(5,1)+a(1,5).* b(5,2)+a(1,4).* b(5,3)+a(1,3).* b(5,4)+a(1,2).* b(5,5)+a(1,1).* b(5,6));% Z^4
c(6,6)=(a(6,6).* b(1,1)+a(6,5).* b(1,2)+a(6,4).* b(1,3)+a(6,3).* b(1,4)+a(6,2).* b(1,5)+a(6,1).* b(1,6)+a(5,6).* b(2,1)+a(5,5).* b(2,2)+a(5,4).* b(2,3)+a(5,3).* b(2,4)+a(5,2).* b(2,5)+a(5,1).* b(2,6)+a(4,6).* b(3,1)+a(4,5).* b(3,2)+a(4,4).* b(3,3)+a(4,3).* b(3,4)+a(4,2).* b(3,5)+a(4,1).* b(3,6)+a(3,6).* b(4,1)+a(3,5).* b(4,2)+a(3,4).* b(4,3)+a(3,3).* b(4,4)+a(3,2).* b(4,5)+a(3,1).* b(4,6)+a(2,6).* b(5,1)+a(2,5).* b(5,2)+a(2,4).* b(5,3)+a(2,3).* b(5,4)+a(2,2).* b(5,5)+a(2,1).* b(5,6)+a(1,6).* b(6,1)+a(1,5).* b(6,2)+a(1,4).* b(6,3)+a(1,3).* b(6,4)+a(1,2).* b(6,5)+a(1,1).* b(6,6));% Z^5
end
.
.
.