Matlab - Problema función "ge"

Hola, estoy haciendo un programa M-File que resuelve un ejercicio que estoy haciendo. El problema es que me da un error matlab y no sé que significa porque intenté solucionarlo pero nada. Os pongo todo el programa (no hace falta entender de que va):

clear;
syms a h k
v=12; %Velocidad
Cp=0.4; %Coeficiente de potencia
n=200; % Revoluciones por minuto
p=3; %Número de palas
P=10000; %Potencia
D=sqrt(8*P/(1.23*Cp*pi*v^3)); % Diámetro
r= D/2; % Radio
alfa=5.5; % Ángulo óptimo
Cl=1.0347; % Coeficiente de sustentación óptimo
Eficiencia=95.5429 ; % Eficiencia maxima
E=1/Eficiencia ; %Epsilon , ángulo de planeo
MATRIZ=[];
for i=0.25:0.25:r
L=2*pi*n*i/(60*v); %Lambda, velocidad específica
%%Inicio del proceso de linealización y resolucion del sistema de
%%ecuaciones no lineales
d=[3 3];
X= [a h k];
f(1)=2*a*L^2*(1+h)*(1+a-(1+h)*L*tan(E))/((1+h)*L+(1+a)*tan(E))+k+k*((1+a)*(1+a-(1+h)*L*tan(E))/((1+h)^2*L^2+(1+h)*(1+a)*L*tan(E)));
f(2)=k+(1+tan(E)^2)*(L^2*(1-a^2)+k*((1-a^2)/((1+h)^2*L)))/(1+tan(E)*(1+a)/((1+h)*L))^2;
f(3)=-h+1+(1-a)*(1+a-(1+h)*L*tan(E))/(L^2*(1+h)+(1+a)*L*tan(E));
b=[0 0 0];
TOLERANCIA=0.00000000000001;
X0=[1 1 1]; % Es el punto inicial que se escoge para iterar
X0=X0';
Jf=jacobian(f,X);
t=1;
j=1;
while (t>=TOLERANCIA) & (j<100)
A=double(subs(Jf,X,X0));
B=b'-subs(f,X,X0)'+A*X0;
S=AB;
X0=S;
d=b-subs(f,X,X0);
t=
t=norma(b-subs(f,X,X0));
j=j+1;
end
a1=X0(1); % Primera fila de X0, es decir, a
h1=X0(2); % Segunda fila de X0, es decir, h
fi=atan((1+a)/((1+h)*L)); %Ángulo de inclinación en radianes
fi_grados=fi*180/pi; %Ángulo de inclinación en grados
cuerda=8*pi*r*cos(E)*sin(fi)^2*(1-a)/(3*Cl*cos(fi-E)*(1+a));
MATRIZ=[MATRIZ;[i,L,a1,h1,fi_grados,cuerda]];
end
disp(' radio , lambda , a , h , ángulo de inclinación , cuerda ')
display(MATRIZ)


El error que me sale en matlab es:

??? Function 'ge' is not defined for values of class 'sym'.

Error in ==> ge at 18
[varargout{1:nargout}] = builtin('ge', varargin{:});

Error in ==> metodonoideal at 31
while (t>=TOLERANCIA) & (j<100)

Muchas gracias

Sólo eliminando la t de la linea de enmedio en:
d=b-subs(f,X,X0);
t=
t=norma(b-subs(f,X,X0));
Y definiendo la S como S=A*B elprograma corrió sin dificultad, obtenéndose:
foromatlab
radio , lambda , a , h , ángulo de inclinación , cuerda

MATRIZ =

[ 1/4, 5*pi/36, -5197115256501860*2^(874), 5028341832536819*2^(872), 180*atan(36/5*(1+a)/(1+h)/pi)/pi, 13610498874333103125/474026950524928*(1+a)/(1+h)^2/pi^2/(25+1296*(1+a)^2/(1+h)^2/pi^2)*(1-a)/cos(atan(36/5*(1+a)/(1+h)/pi)-3016763947417461/288230376151711744)]
[ 1/2, 5*pi/18, NaN, NaN, 180*atan(18/5*(1+a)/(1+h)/pi)/pi, 13610498874333103125/1896107802099712*(1+a)/(1+h)^2/pi^2/(25+324*(1+a)^2/(1+h)^2/pi^2)*(1-a)/cos(atan(18/5*(1+a)/(1+h)/pi)-3016763947417461/288230376151711744)]
[ 3/4, 5*pi/12, 7398517467446094*2^(457), -4858821356317753*2^(457), 180*atan(12/5*(1+a)/(1+h)/pi)/pi, 1512277652703678125/474026950524928*(1+a)/(1+h)^2/pi^2/(25+144*(1+a)^2/(1+h)^2/pi^2)*(1-a)/cos(atan(12/5*(1+a)/(1+h)/pi)-3016763947417461/288230376151711744)]
[ 1, 5*pi/9, 4537009152876595*2^(556), -7850639180538394*2^(555), 180*atan(9/5*(1+a)/(1+h)/pi)/pi, 13610498874333103125/7584431208398848*(1+a)/(1+h)^2/pi^2/(25+81*(1+a)^2/(1+h)^2/pi^2)*(1-a)/cos(atan(9/5*(1+a)/(1+h)/pi)-3016763947417461/288230376151711744)]
[ 5/4, 25*pi/36, 8032490316826180*2^(632), -8592629075634317*2^(632), 180*atan(36/25*(1+a)/(1+h)/pi)/pi, 13610498874333103125/474026950524928*(1+a)/(1+h)^2/pi^2/(625+1296*(1+a)^2/(1+h)^2/pi^2)*(1-a)/cos(atan(36/25*(1+a)/(1+h)/pi)-3016763947417461/288230376151711744)]
[ 3/2, 5*pi/6, 6171635246381203*2^(696), -7832670380295430*2^(696), 180*atan(6/5*(1+a)/(1+h)/pi)/pi, 1512277652703678125/1896107802099712*(1+a)/(1+h)^2/pi^2/(25+36*(1+a)^2/(1+h)^2/pi^2)*(1-a)/cos(atan(6/5*(1+a)/(1+h)/pi)-3016763947417461/288230376151711744)]
[ 7/4, 35*pi/36, 8752271225214568*2^(750), -6400251204213341*2^(751), 180*atan(36/35*(1+a)/(1+h)/pi)/pi, 13610498874333103125/474026950524928*(1+a)/(1+h)^2/pi^2/(1225+1296*(1+a)^2/(1+h)^2/pi^2)*(1-a)/cos(atan(36/35*(1+a)/(1+h)/pi)-3016763947417461/288230376151711744)]
[ 2, 10*pi/9, 7785460692745122*2^(798), -6420327906467066*2^(799), 180*atan(9/10*(1+a)/(1+h)/pi)/pi, 13610498874333103125/7584431208398848*(1+a)/(1+h)^2/pi^2/(100+81*(1+a)^2/(1+h)^2/pi^2)*(1-a)/cos(atan(9/10*(1+a)/(1+h)/pi)-3016763947417461/288230376151711744)]
[ 9/4, 5*pi/4, 6543274246746450*2^(841), -5984037750588861*2^(842), 180*atan(4/5*(1+a)/(1+h)/pi)/pi, 1512277652703678125/4266242554724352*(1+a)/(1+h)^2/pi^2/(25+16*(1+a)^2/(1+h)^2/pi^2)*(1-a)/cos(atan(4/5*(1+a)/(1+h)/pi)-3016763947417461/288230376151711744)]
[

Muchas gracias, perdona la torpeza que no me había fijado.

No te preocupes, es común que cosas así se nos pasen.

Hola, después de haber corregido el error, el progrma en matlab me funcionaba, pero a veces me vuelve a aparecer. Lo que hago es apagar matlab y volver a encender. No es la mejor forma, pero es que no sé porque me sale el error que te había comentado.

Un saludo.

Lo he corrido cinco veces consecutivas y en todas ha funcionado sin problemas, sé que eso no es una prueba contundente de su buen funcionamiento, de hecho es más contundente la falla que comentas ´porque esa sí demuestra, de manera definitiva, que existe algún error en el programa que hay que corregir.
¿Siempre te manda el mismo código de error?

Pues sí, y siempre me sale cuando inicio por primera vez el programa, cuando apago el matlab y lo vuelvo a encender ya no me sale.

El error es este:

??? Function 'ge' is not defined for values of class 'sym'.

Error in ==> ge at 18
[varargout{1:nargout}] = builtin('ge', varargin{:});

Error in ==> metodonoideal at 32
while (t>=TOLERANCIA) & (j<100)