 function [x,y,z]=path(t,gravity,Vo, theta, phi)
 
global rho_air
%define things specific to the drag
cross_section=1e-6; %m
drag_coefficient=.47;
density_of_ice=.9%kg/m^3
mass_of_material=density_of_ice*(4/3)*pi*(cross_section.^3);


     
%need to find the velocity of the material at every time step
V=Vo*(sind(theta))-(gravity*(t))%integrate to get v at this ;


%define drag force
        Drag_force=.5*rho_air*cross_section*drag_coefficient*((V.^2))  ;  
        drag_accel=Drag_force/mass_of_material;

     




        
        
  %assume gravity only in the z direction      
        x = Vo*cosd(theta)*t ;
        y=Vo*cosd(phi)*t; 
  
        z = Vo*(sind(theta)*t)-(0.5*gravity*(t.^2))%-(0.5*drag_accel.*(t.^2));
       
     
        hold on
  
end