Question

Please don't attempt if you can not write a MATLAB PROGRAM.
I don't need numerical physics solution, I need MATLAB code.
Read the question carefully

2) Falling Drops Imagine a solid sphere of diameter, d, falling through the air. After its release, the sphere accelerates until it reaches its terminal velocity. The terminal velocity is obtained when the drag force due to air friction, Fa, balances the weight of the sphere, W. The weight of the sphere is simply its mass, m, times acceleration due to gravity (g- 9.8 m/s2), and the drag force can be computed as where Cd is a dimensionless empirical drag coefficient, p is the density of the air, v is the veloci the sphere, and A π*d2/4 is the frontal area of the falling sphere. ty of 1 of 4

Answer 1

clc
clear
close all

R = 287;
P = 101300;
m = 0.5;
g = 9.8;
W = m*g;
d = 15e-2; %m
A = pi*d^2/4;
b1 = 2.156954e-14; b2 = -5.332634033e-11; b3 = 7.477905983e-8; b4 = 2.527878788e-7;

mew = @(T) b1*T^3 + b2*T^2 + b3*T + b4;
rhow = @(T) P/(R*T);
Re = @(v,T) rhow(T)*v*d/mew(T);
cd = @(v,T) 24/Re(v,T) + 6/(1+sqrt(Re(v,T)));

i = 1;

initial_T = 273;
for T = (-60:60)+ 273.15;
fv = @(v) cd(v,T)*0.5*rhow(T)*v^2*A - W
v(i) = fzero(fv,initial_T)
i = i+1;
end

T = (-60:60)+ 273.15;
plot(T,v)