clear

t = [-10:10];
alfa = %pi/8;
H = 2;
L = 1;                     //dlugosc liny
G = 10;
X = [1:21];
Y = [1:21];
Z = [1:21];
r = 10;                    //duze R z obrazka
r2 = r + sin(alfa) * L     //male r
V0 = sqrt(2*r2 * G / sin(2*alfa));  //Vo liczona ze wzoru na zasieg(zasieg ma byc 2r)
h = H - cos(alfa) * L

for i=[1:21] do
  X(i) = cos(alfa) * V0 * t(i) - cos(alfa) * V0 * t(1) - r;
  Y(i) = sin(alfa) * V0 * t(i) - 1/2 * G * t(i)^2 - sin(alfa) * V0 * t(1) + 1/2 * G * t(1)^2 + h;
  Z(i) = 0;
end;

param3d(X, Y, Z);

for i=[1:21] do
  X(i) = r * cos(alfa*i);
  Y(i) = r * sin(alfa*i);
  Z(i) = h;
end;

param3d(X, Z, Y);