Untitled

g = 9.81;
rho = 1.225;
S = 11.92;
m = 1532.55;
C_W = 0.021;
k = 1.075 / (8.5 * pi);
m_r = 0.55;
C_A = 0.2;

f = @(V) (m * V) ./ (-m_r * m * g - (C_W - m_r * C_A)*(rho/2) * V.^2 * S);

result = integral(f,39.223,0);

disp(result)