clear all
close all
syms pi000 pi001 pi002 pi003 pi010 pi011 pi012 pi013 pi020 pi021 pi022 pi023 pi100 pi101 pi102 pi103 pi110 pi111 pi112 pi113 pi120 pi121 pi122 pi123 a1 a2 a3 b11 b21 b22 b31 b32 b33
eq1=(a1+a2+a3)*pi000-b11*pi100-b21*pi010-b31*pi001-b22*pi020-b32*pi002-b33*pi003;
eq2=(a1+a2+b31)*pi001-b11*pi101-b21*pi011-b22*pi021;
eq3=(a1+a2+b32)*pi002-b11*pi102-b21*pi012-b22*pi022;
eq4=(a1+a2+b33)*pi003-b11*pi103-b21*pi013-b22*pi023-a3*pi000;
eq5=(a1+a2+a3+b21)*pi010-b11*pi110-b31*pi011-b32*pi012-b33*pi013;
eq6=(a1+b21+b31)*pi011-b11*pi111;
eq7=(a1+b21+b32)*pi012-b11*pi112-a2*pi010;
eq8=(a1+b21+b33)*pi013-b11*pi113-a3*pi010;
eq9=(a1+a2+a3+b22)*pi020-b11*pi120-b31*pi021-b32*pi022-b33*pi023-a2*pi000;
eq10=(a1+b31+b22