Задачи по моделированию - Pascal - Ответ 3015716

uses crt;
var m,h,v,g,p,mu,c,sf,a,s,dt,t,r,k1,k2:extended;
f:text;
begin
assign(f,'c:\bomba.txt');
rewrite(f);
write(f);
m:=300;
t:=0;
s:=0;
g:=9.8;
v:=0;
a:=g;
r:=1;
mu:=1.002;
c:=0.4;
sf:=pi*r*r;
p:=1000;
k1:=6*pi*mu*r;
k2:=0.5*c*sf*p;
dt:=1;
repeat
t:=t+dt;
s:=s+v*dt+(a*dt*dt)/2;
v:=v+a*dt;
a:=(m*g-k1*v-(k2*v*v))/m;
write(f,t:10:3);
writeln(f,v:10:3);
until  t=200;
close(f);
readln;
end.

var t, x, y, al, ald, dt, ax, ay, m, v, vx, vy, g, pi, mu, r, c, po, sp, k1, k2: extended;
f:text;
begin
assign (f,'f.txt');
rewrite (f);
 
m:=20;
v:=70;
dt:=1;
ald:=0.52;
g:=9.78;
pi:=3.14;
k1:=3.8;
k2:=0.032;
al:=0;
 
 
repeat
al:=al+ald;
vx:=v*cos(180/(al*pi));
vy:=v*sin(180/(al*pi));
t:=0;
x:=0;
y:=1;
ax:=0;
ay:=-g;
 
       repeat
       t:=t+dt;
       x:=x+vx*dt+ax*(dt*dt)/2;
       y:=y+vy*dt+ay*(dt*dt)/2;
       vx:=vx+ax*dt;
       vy:=vy+ay*dt;
       ax:=-(k1*vx-k2*(vx*vx))/m;
       ay:=(m*g-k1*vy-k2*(vy*vy))/m;
       until x=200;
 
If (y>=15) and (y<=17) then
   begin
   write (f,t:10:3);
   write (f,x:10:3);
   write (f,y:10:3);
   writeln (f,al:10:3);
   end;
 
until y>17;
close (f);
end.