I'm aware that there are a few similar questions already answered, but I could not find what I was looking for in any of then, so please bear with me :)
For a school project I need to find an angle ($\theta$) of maximum range ($R$) and its dependency on initial velocity $\theta(v_0)$, considering quadratic drag force and wind blowing in the direction opposite of the throw.
Drag force in $x$ direction:$F_{\textrm{drag},x}=-Kv^2\cos(\alpha)=-Kv_x(v_x^2+v_x^2)^{1 /2}$
Drag force in $y$ direction:$F_{\textrm{drag},y}=-Kv^2\sin(\alpha)=-Kv_y(v_x^2 +v_x^2)^{1/2}$
Wind force: $F_\textrm{wind}=F$
In x direction: $F_{w,x}=F\cos(\theta)$
In y direction: $F_{w,y}=F\sin(\theta)$
So the equations should look like this:
\begin{align}mv_y' &= -kv_x\cdot (v_x^2 + v_y^2)^{1/2}- F\cdot \cos\theta \\ mv_x' &= -kv_y\cdot (v_x^2 + v_y^2)^{1/2}- F\cdot \sin\theta- mg \end{align} Now I'm not sure how to solve this system of two non-linear differential equations... I'm figuring it should be done numerically, but I'm not very familiar with programming (I've only programmed a little bit in Python, but not this kind of stuff). Thank you for any answer!