//Solution 11-5
WD=get_absolute_file_path('11_06_solution.sce');
datafile=WD+filesep()+'11_06_example.sci';
clc;
exec(datafile)
//unit conversions
D = D / 100; //from [cm] to [m]
V = V / 3.6 //from [km/h] to [m/s]
P = P * 1.01325 * 10**5; //from [atm] to [Pa]
T = T + 273; //from [C] to [K]
rho_air = P / (R * T); //from ideal gas equation
A = %pi / 4 * D**2; //frontal area of ball
omega = 2 * %pi * N / 60; //angular velocity of ball [rad/s]
nd_rotation = omega * D /(2 * V); //non dimensional rate of rotation
//from figure 11-53 lift coefficient coefficient corresponding to nd_rotation is
C_L = 0.21;
F_L = C_L * A * rho_air * V**2 / 2; //drag force
W = m * g;
if W > F_L then
    printf("Ball will drop under the combined effect of gravity and lift due to spinning with net force of %1.3f - %1.3f = %1.3f N", W, F_L, W - F_L);
else
    printf("Ball will rise under under the combined effect of gravity and lift due to spinning with net force of %1.3f - %1.3f = %1.3f N", F_L, W, F_L - W);
end