matlab求雅可比矩阵特征值

根据已知条件，可以得到雅可比矩阵。目的是想求雅可比矩阵的特征值的实部最大值。遇到的问题：雅可比矩阵是5行5列矩阵，但在j2=eval(j1) 这里出现了4行5列矩阵，这是为什么呢。请大家指教

format long
syms x1 x2 x3 x4 x5 C B   
g1=(x2/(x2+0.51))*(1-x2/23)*(1-x3/46)*(1-x4/12)*(1-5*x5/3);
g2= 1.67+C/0.06+3.291*(1/(x2+7.8));
g3=-2.623+C*80.6+23.97*(22/(x2+13.6));
g4=-0.19+C*30.9+5.4*(x2/(5*x2+7.78));
g5=g2*(0.019/(0.09+C*x2)+8.8/(5.85+12*C*x2));
%变量常微分方程组f
f=[(g1-C)*x1;
  C*(B-x2)-g2*x1;
  g3*x2-C*x3;
  g4*x1-C*x4;
  g5*x1-C*x5;
  ];
v=[x1 x2 x3 x4 x5];
j=jacobian(f,v)       %对常微分方程组f求雅克比矩阵
load tr     %读取存储数据  
[m0,n0]=size(T1);    %T1为tr里的数据
for n=1:n0
    n;
P6=double(T1(:,n))       
j1=subs(j,[x1,x2,x3,x4,x5,C,B]',P6)        
j2=eval(j1)                                                                                                                      
e=eig(j2)                                  
maxreale=max(real(e))                   
end

出现的错误

错误使用 symengine
A square matrix is expected.

出错 sym/privUnaryOp (line 943)
            Csym = mupadmex(op,args{1}.s,varargin{:});

出错 sym/eig (line 51)
    V = privUnaryOp(A, 'symobj::eigenvalues');

出错 wdxfx1 (line 29)
e=eig(j2)           

这个不会

tr 是什么？