MATLAB: What am I doing wrong

doubleequationMATLABsolvesyms

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!!!

Best Answer

 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 )
Related Question