What am I doing wrong?

Hello!

I am trying to solve an equation but I get this error:

Error using sym>notimplemented (line 2682)
Function 'ge' is not implemented for MuPAD symbolic objects.

Error in sym/ge (line 823)
            notimplemented('ge');

Error in Programa_drag_torque_paper (line 41)
    a=double(s{i}(1)>=r_1)&&double(s{i}(1)<=r_2);

I think that the problem is related with how are determined the variables (sym, double,...) but I don't know how to fix it.

Here is the code:

clear all
clc

r_1=178*10^-3;                      
r_2=218*10^-3;                      
h=0.2*10^-3;                        
N=10;                               
Q_o=5.83*10^-5;                     
rho=820;                            
mu=0.0082 
sigma=0.030
omega=[700:100:2200]*(2*pi/60);     
omega_rpm=[700:100:2200];           
theta=0*(2*pi/60);                  
r_m=(r_1+r_2)/2;                    

Re_h=(rho*omega*r_1*h)/mu;         

m=zeros(16,1);
G_r=m;
f=m;
syms r_o_1
syms r
s=cell(16,1);
r_o_final_1=m;
T_1=m;
t_final_1=m;

for i=1:16

      if Re_h(i)>=500
          f(i)=0.885*Re_h(i)^-0.367;            
      else
          f(i)=0.09;
      end
      G_r(i)=1/12*(1+(0.00069.*Re_h(i).^0.88));    

      s{i}=solve(((rho*(omega(i)^2/2))*(f(i)+(1/4))*(r_o_1^2))-((mu*Q_o*r_o_1)/(2*pi*r_m*h^3*G_r(i)))-((rho*(omega(i)^2/2))*(f(i)+(1/4))*r_1^2)-((2*sigma*cos(theta))/h)+((mu*Q_o*r_1)/(2*pi*r_m*h^3*G_r(i))),r_o_1);

      a=double(s{i}(1)>=r_1)&&double(s{i}(1)<=r_2);
      b=double(s{i}(2)>=r_1)&&double(s{i}(2)<=r_2);

      if a
          r_o_final_1(i)=s{i}(1);
      elseif b
          r_o_final_1(i)=s{i}(2);
      else
          r_o_final_1(i)=r_2;
      end

      T_1(i)=int(((mu*omega(i)*r^3)/h)*(1+(0.0012*Re_h(i)^0.94)),r_1,r_o_final_1(i));        
      T_final_1(i)=2*N*pi*T_1(i);         

end

It would be great if you could find a solution.

Thank you so much for your answers!!!

Answer 1

[Walter Roberson]

 a=double(s{i}(1)) >=r_1&&double(s{i}(1)) <=r_2;  
 b=double(s{i}(2))>=r_1&&double(s{i}(2)) <=r_2;

Notice I moved the )

[Lorena Ortiz]

Thank you Walter Roberson!!

Now it works :)