examen de polo

modelo=input('introduce modelo de coginete   ');
a=xlsread('tablas.xlsx');
b=xlsread('tablas.xlsx',2);
[r,c]=find(a==modelo);
co=a(r,c+1)
Pa=input('introduce el valor de la carga axial  ');
Pr=input('introduce el valor de la carga radial  ');
e=Pa/co
valor=Pa/Pr

e=input('introduce el error    ');

if valor>e
   [r,c]=find(b==e);
   x=b(r,c+1)
   y=b(r,c+2)
end
if valor<=e
   [r,c]=find(b==e) ;
    x=b(r,c+3)
   y=b(r,c+4)
end

    pzero=x*Pr+y*Pa
    tao=(a(r,c+1)/pzero)^3
    Nb=input('introduce la velocidad rotacional    ');

    t10=(10E6*tao)/(60*Nb)