can somebody please help me with matlab.
clear all;
clc;

den_b = 8550;
den_f = 1100;
visc_f = 0.1;
D_b = 0.0064;
g = 9.81;
u = 0.1;
step = 0.01;
n = 100;
tol = 10^-4;

Re = @(u)(D_b*abs(u)*den_f)/visc_f
bouyancy = (den_f-den_b)*(pi/6).*D_b.^3.*g
Drag = @(u)((2.25./Re(u).^0.31)+(0.358.*Re(u).^0.06)).^3.45
CirArea =(pi/4).*D_b^2
f =@(u) bouyancy-(Drag(u).*CirArea*den_f.*u*abs(u)/2)

for ii = 1:n
u(ii+1) = u(ii)-f(u(ii))/finitediff(u(ii))

Question

can somebody please help me with matlab.
clear all;
clc;

den_b = 8550;
den_f = 1100;
visc_f = 0.1;
D_b = 0.0064;
g = 9.81;
u = 0.1;
step = 0.01;
n = 100;
tol  = 10^-4;

Re = @(u)(D_b*abs(u)*den_f)/visc_f
bouyancy = (den_f-den_b)*(pi/6).*D_b.^3.*g
Drag = @(u)((2.25./Re(u).^0.31)+(0.358.*Re(u).^0.06)).^3.45
CirArea =(pi/4).*D_b^2
f =@(u) bouyancy-(Drag(u).*CirArea*den_f.*u*abs(u)/2)

for ii = 1:n
    u(ii+1) = u(ii)-f(u(ii))/finitediff(u(ii))<tol
end

uarray = 0.1:0.1:100;
yarray = f(uarray);
plot(urray,yarray)

anonymous · Answer

what wrong

anonymous · Answer

@TIMDAVE

anonymous · Answer

this code is not working. in my command window it prints out the functions instead of evaluating it.

mathmate · Answer

@Lion.k
If this is the same problem, it would be nice if you could close this, or refer to the other post.