%Datos a IngresarQ= 0.5; % caudal en (m3/s)L=1200; % longitud de la tubera en (m)T=15; % temperatura en (oC)ks=0.0015; % rugusidad del tubo (m2/s)hf=81.67; % prdida (m)Tol=1E-8; % toleranciaImax=12; % numero de iteraciones mximas%Clculo de la Viscosidad en funcin de lamperaturanu=( 1.14-0.031*(T-15)+0.00068*(T-15)^2)*1E-6;%Clculo de kk=9.806*pi^2/(8*L*Q^2);%Datos iniciales para la iteracind=0.1; %valor semilla del dimetroError=1; %error inicialcont=0; %contador inicial%Proceso Iterativowhile(Error>=Tol);D=d;%clculo del nmero de ReynoldRe=4*Q/(pi*D*nu);%clculo del dimetrod=D-((-2*log10( ks/(3.7*D)+ 2.51/(Re*sqrt(k*hf*D^5))))^(-2) - k*hf*D^5)/...((log(10)^2*((10*ks)/(37*D^2) + (251*D^4*hf*k)/(40*Re*(D^5*hf*k)^(3/2))))/...(2*log(251/(100*Re*(D^5*hf*k)^(1/2)) + (10*ks)/(37*D))^3*...(251/(100*Re*(D^5*hf*k)^(1/2)) + (10*ks)/(37*D))) - 5*D^4*hf*k);Error=abs(D-d);cont=cont+1;fprintf('Iteracin nmero %g\n', cont)fprintf('Valor del Dimetro %4.8f\n\n', d)if(cont>=Imax);breakendend