function[hmax, dmax] = trajectory(vel_inicial, angulo)

  vel_inicial_x = vel_inicial * cos(angulo);
  vel_inicial_y = vel_inicial * sen(angulo);

  gravidade = 9.81;

  hmax = (vel_inicial_y^2)/(2*gravidade);

  tempo_para_hmax = vel_inicial_y/g;
  tempo_trajetoria = tempo_para_hmax * 2;

  dmax = vel_inicial * tempo_trajetoria;

  plot(hmax, dmax)