Hi Alexander,
Have a look at the following code. It shows one of the ways by which you can solve your question.
D2 = 11.4;
H2 = 25.2;
g = 9.81;
A2 = 1:1:89;
V = zeros([89 1]);
syms V2
for i = 1:89
eqn = D2-(((V2)^2*(sind(2*A2(i))))/(2*g))*(1+(1+((2*g*H2)/((V2)^2*(sind(A2(i)))^2)))^0.5)==0;
S = solve(eqn,V2);
if S(1)>=0
V(i) = S(1);
else
V(i) = S(2);
end
end
plot(A2,V)
Output:
For each value of A2, two values of V2 are obtained, out of which I have selected the positive values for plotting. You can also select the negative values.
For additional information on the 'solve' function, refer to the following link :
Regards,
Divija
