Под каким углом к горизонту следует бросить камень, чтобы дальность полета была наибольшей

program model1;
uses graphABC;
const g=9.8; pi=3.14;
var x,y,k,m,t,v0,vx,vy,h,dt,alpha:real;
    xscreen,yscreen:integer;
procedure input(var k,m,alpha,v0:real);
begin
writeln('k'); readln(k);
writeln('m'); readln(m);
writeln('alpha'); readln(alpha);
writeln('v0'); readln(v0);
end;
procedure inic(var x,y,t,vx,vy:real; var xscreen,yscreen:integer);
begin
x:=0;y:=0;t:=0;
vx:=v0*cos(alpha*pi/180);
vy:=v0*sin(alpha*pi/180);
xscreen:=0; yscreen:=450;
moveto(xscreen,yscreen);
dt:=0.01;
end;
function f1(t,x,y,vx,vy:real):real;
begin
f1:=vx;
end;
function f2(t,x,y,vx,vy:real):real;
begin
f2:=vy;
end;
function f3(t,x,y,vx,vy:real):real;
begin
f3:=-k*vx/m;
end;
function f4(t,x,y,vx,vy:real):real;
begin
f4:=-g-k*vy/m;
end;
procedure PK4(var x,h:real;var y1,y2,y3,y4:real);
var k11,k12,k13,k14,k21,k22,k23,k24,k31,k32,k33,k34,k41,k42,k43,k44:real;
begin
k11:=h*f1(x,y1,y2,y3,y4);
k12:=h*f2(x,y1,y2,y3,y4);
k13:=h*f3(x,y1,y2,y3,y4);
k14:=h*f4(x,y1,y2,y3,y4);
k21:=h*f1(x+h/2,y1+k11/2,y2+k12/2,y3+k13/2,y4+k14/2);
k22:=h*f2(x+h/2,y1+k11/2,y2+k12/2,y3+k13/2,y4+k14/2);
k23:=h*f3(x+h/2,y1+k11/2,y2+k12/2,y3+k13/2,y4+k14/2);
k24:=h*f4(x+h/2,y1+k11/2,y2+k12/2,y3+k13/2,y4+k14/2);
k31:=h*f1(x+h/2,y1+k21/2,y2+k22/2,y3+k23/2,y4+k24/2);
k32:=h*f2(x+h/2,y1+k21/2,y2+k22/2,y3+k23/2,y4+k24/2);
k33:=h*f3(x+h/2,y1+k21/2,y2+k22/2,y3+k23/2,y4+k24/2);
k34:=h*f4(x+h/2,y1+k21/2,y2+k22/2,y3+k23/2,y4+k24/2);
k41:=h*f1(x+h,y1+k31,y2+k32,y3+k33,y4+k34);
k42:=h*f2(x+h,y1+k31,y2+k32,y3+k33,y4+k34);
k43:=h*f3(x+h,y1+k31,y2+k32,y3+k33,y4+k34);
k44:=h*f4(x+h,y1+k31,y2+k32,y3+k33,y4+k34);
y1:=y1+1/6*(k11+2*k21+2*k31+k41);
y2:=y2+1/6*(k12+2*k22+2*k32+k42);
y3:=y3+1/6*(k13+2*k23+2*k33+k43);
y4:=y4+1/6*(k14+2*k24+2*k34+k44);
x:=x+h;
end;
begin
input(k,m,alpha,v0);
inic(x,y,t,vx,vy,xscreen,yscreen);
repeat
PK4(t,dt,x,y,vx,vy);
xscreen:=trunc(x*50);
yscreen:=trunc(450-y*50);
lineto(xscreen,yscreen);
t:=t+dt;
until y<=0;
end.