Crear cuenta
syms f Re
a=(1./(sqrt(f)))+(0.4)-(1.74*log(Re*(sqrt(f))));
b=solve(a,Re) % b =exp(50/(87*f^(1/2)) + 20/87)/f^(1/2)
quiero reemplazar f en b
f=10^4:10^4:10^6;
n=length(f); %tamaño del vector
fprintf('Re f \n')
disp('=======================')
for i=1 :n
fprintf(' % 3.3f %1.6f\n',Re(i) ,f(i))
end
disp('==================')
syms f Re
a=(1./(sqrt(f)))+(0.4)-(1.74*log(Re*(sqrt(f))));
b=solve(a,Re) % b =exp(50/(87*f^(1/2)) + 20/87)/f^(1/2)
%quiero reemplazar f en b
B=vectorize(inline(b));
f=10^4:10^4:10^6;
A=B(f)
n=length(f); %tamaño del vector
fprintf('Re f \n')
disp('=======================')
for i=1 :n
fprintf(' % 3.3f %1.6f\n',Re(i) ,f(i))
end
disp('==================')
syms f Re
a=(1./(sqrt(f)))+(0.4)-(1.74*log(Re*(sqrt(f))));
b=solve(a,Re) % b =exp(50/(87*f^(1/2)) + 20/87)/f^(1/2)
%quiero reemplazar f en b
B=vectorize(inline(b));
f=10^4:10^4:10^6;
A=B(f)
n=length(f); %tamaño del vector
Re=sym(rand(1,n))
fprintf('Re f \n')
disp('=======================')
for i=1 :n
fprintf(' % 3.3f %1.6f\n',double(Re(i)) ,f(i))
end
disp('==================')
>> simnbolicas
Warning: The solutions are valid under the following conditions: -pi <= (50*sin(angle(f)/2))/(87*abs(f)^(1/2)) & (50*sin(angle(f)/2))/(87*abs(f)^(1/2)) < pi. To include parameters
and conditions in the solution, specify the 'ReturnConditions' option.
> In solve>warnIfParams (line 517)
In solve (line 360)
In simnbolicas (line 3)
Re f
=======================
0.764 10000.000000
0.559 20000.000000
0.184 30000.000000
0.498 40000.000000
0.518 50000.000000
0.994 60000.000000
0.855 70000.000000
0.962 80000.000000
0.679 90000.000000
0.404 100000.000000
0.935 110000.000000
0.479 120000.000000
0.232 130000.000000
0.396 140000.000000
0.705 150000.000000
0.559 160000.000000
0.757 170000.000000
0.995 180000.000000
0.962 190000.000000
0.535 200000.000000