function [x,y] = parabolicTrajectory(x0,y0,v0x,v0y,t)
% Computes a parabolic trajectory under the Earth's gravity.
% Input:
%   x0,y0   : Initial horizontal and vertical position, in meters
%   v0x,v0y : Initial velocity, in m/s
%   t       : vector of times, in s
% Output:
%   x,y    : Vectors of x and y coordinates at times t, in meters
% Usage: 
%   [x,y] = parabolicTrajectory(x0,y0,v0x,v0y,t)
% J.M.Soler. Sep.2010
g=9.81;                  % Gravitational acceleration (m/s^2)
x=x0+v0x*t;              % x values at times t, in m
y=y0+v0y*t-0.5*g*t.^2;   % y values at times t, in m
end