Массив для расчета по формулам

uses GraphABC;
 
const
  g = - 9.81;
  s = 10; // масштаб
 
var
  t, dt, x, y, vx, vy: real;
 
begin
 
  // НУ
  t := 0;
  dt := 0.001;
  x := 0;
  y := 0;
  vx := 10;
  vy := 20;
  
  // Интегрирование равнений движения методом Эйлера
  repeat
    x += vx * dt;
    y += vy * dt;
    vy += g * dt;
    t += dt;
    circle(round(x*s) + 100, -round(y*s) + 300, 3);
  until y < -0.1;
  
end.

 uses GraphABC;
 
const
  g = - 9.81;
  s = 20; // масштаб
  n = 19;
  w = 29;
  
type
  mas1 = array[10..n] of integer;
  mas2 = array[20..w] of integer;
 
var
  t, dt, x, y: real;
  vx: mas1;
  vy: mas2;
  i: byte;
  
begin
  t := 0;
  dt := 0.001;
  x := 0;
  y := 0;  
  writeln('10, 11, 12, 13, 14, 15, 16, 17, 18, 19');
  for i:=1 to n do readln(vx[i]);
  writeln(('10, 11, 12, 13, 14, 15, 16, 17, 18, 19:'));
  for i:=1 to n do write(vx[i]:19);
  writeln('20, 21, 22, 23, 24, 25, 26, 27, 28, 29');
  for i:=1 to w do readln(vy[i]);
  writeln(('20, 21, 22, 23, 24, 25, 26, 27, 28, 29:'));
  for i:=1 to w do write(vy[i]:29);
  
repeat
    x += vx * dt;
    y += vy * dt;
    vy += g * dt;
    t += dt;
    circle(round(x*s) + 100, -round(y*s) + 300, 3);
  until y < -0.1;
  
end.

uses GraphABC;
const
  g = - 9.81;
  s = 10; // масштаб
  n=10;
type skor=record
          vx,vy:real;
          end;
var v:array[1..n]of skor;
    t, dt, x, y: real;
begin
setwindowsize(800,400);
centerwindow;  
t := 0;
dt := 0.001;
randomize;
for var i:=1 to n do
 begin 
  v[i].vx := 5+random(20);
  v[i].vy := 5+random(20);
 end; 
// Интегрирование равнений движения методом Эйлера
for var i:=1 to n do
  begin
   x:=0;
   y:=0;
   setpencolor(clrandom);
   repeat
    x += v[i].vx * dt;
    y += v[i].vy * dt;
    v[i].vy += g * dt;
    t += dt;
    circle(round(x*s) + 100, -round(y*s) + 300, 3);
  until y < -0.1;
 end; 
  
end.