Calculo Del Roi

clear all;
clc;
% Este programa es para calcular el NRe y utilizaremos el ciclo if then
% else
display ('Este programa calcula el NRe en el Sistema Internacional');
display (' Ydetermina el regimen del flujo ');
display (' NRe = (diametro*densidad*velocidad/viscocidad ');
display(' ');
res=input ('Deseacorrer el programa? (s para si, cualquier letra para no) ', 's');
while (res=='s'|res=='S')
pp=input ('Para realizar el calculo en el Sistema Internacional presiona 1, para elSistema Ingles presiona 2: ');
if pp==1
diametro=input ('Cual es el valor del diametro? (m): ');
densidad=input ('Cual es el valor de la densidad? (kg/m^3): ');velocidad=input ('Cual es el valor de la velocidad? (m/s): ');
viscocidad=input ('Cual es el valor de la viscocidad? (kg/m.s): ');
Re=diametro*densidad*velocidad/viscocidad;fprintf ('El numero de Reynold es:%3.2f y es un valor adimensional ', Re);
else
if pp==2
diametro=input ('Cual es el valor del diametro? (ft): ');densidad=input ('Cual es el valor de la densidad? (lb/ft^3): ');
velocidad=input ('Cual es el valor de la velocidad? (ft/s): ');
viscocidad=input ('Cual es el valor de laviscocidad? (lb/ft.s): ');
Re=diametro*densidad*velocidad/viscocidad;
fprintf ('El numero de Reynold es:%3.2f y es un valor adimensional ', Re);
if(Re<2100)display('El regimen del flujo es laminar');
else
if(Re>4100)
display('El regimen del flujo es turbulento');
elsedisplay('El regimen del flujo es de transición');
end
end
res=input('Desea volver a correr el programa? (s para sí, cualquier letra para no) ', 's');
end
end
end