//Karyn Czapnik HW #7
#include <iostream>
#include <math.h>
#include <iomanip>

   using namespace std;

   double theNewVelocity (double dragCoefficient, double velocity);
   double yposition (double velocity, double angle, double time);
   double angle = 40;
   double pi = 3.14;
   double rad = angle*(pi/180);

    int main ()
   {
      double drag, velocity, position, time;
      cout << "Enter the coefficient of drag (between 1e-4 and 1e-6):";
      cin >> drag;
      velocity = 1200;
      int spaces = int(yposition (velocity, rad, time));
      cout << setw(spaces) << "* \n";
      cout << spaces << endl;
      double newVelocity = theNewVelocity (drag, newVelocity);

      for (time=2; time<50; time=time+2)
      {
         if (time==2){
            spaces = int (yposition (newVelocity, rad, time));
            cout << setw(spaces) << "*\n";
            cout << spaces << endl;
         }
         else {
            newVelocity = theNewVelocity (drag, newVelocity);
            spaces = int (yposition (newVelocity, rad, time));
            cout << setw(spaces) << "* \n";
            cout << spaces << endl;
         }
      }

   }

    double theNewVelocity (double dragCoefficient, double velocity)
   {
      double velocity1 = dragCoefficient*(velocity*velocity);
      double velocitynow = velocity-velocity1;
      return velocitynow;
   }

    double yposition (double velocity, double angle, double time)
   {
      double a = 9.81;
      double yposition = (velocity*time*sin(rad))-(.5*a*(time*time));
      return yposition;
   }