Matlab - Sist. de ecuaciones no lineales

Hola, estoy absolutamente entrampado con la solución de un sistema de ecuaciones no lineales. Yo sé que tiene solución (9 ecs., 9 incógnitas), pues lo vie en un libro y quiero replicarlo para aplicar la técnica a un sistema más grande.

Éste es el programa que hice, con ayuda que ya pedí por acá:

syms b10 b20 b12 c11 c12 c21 c22 st sr
B=[1 b12 ; 0 1];
B0=[b10; b20];
C1=[c11 c12 ; c21 c22];
E=[st 0 ; 0 sr];
A0=inv(B)*B0
A1=inv(B)*C1
e=inv(B)*E*(inv(B))'

a10=A0(1)
a20=A0(2)

a11=A1(1,1)
a12=A1(1,2)
a21=A1(2,1)
a22=A1(2,2)

ett=e(1,1)
err=e(2,2)
etr=e(2,1)

[b10,b20,b12,c11,c12,c21,c22,st,sr]=solve('a10 =b10-b12*b20','a20 =b20','a11 =c11-b12*c21','a12 =c12-b12*c22','a21 =c21','a22 =c22','ett =st+b12^2*sr','err =sr','etr =-b12*sr')

Necesito encontrar las soluciones para los elementos de las matrices B, B0, C1 y E, para después aplicar la solución a valores numéricos de A0, A1 y e.

Agradeceré cualquier pista que me puedan dar. Muchos saludos.

Pablo