global DeltaCav DeltaFinE Theta
  Theta=0;DeltaCav=0;DeltaFinE=0;
   y0=[60 0 0 Theta*pi/180 0 0];%給函數賦初值
   tspan=[0 10];
     [t,y]=ode45('projectiledes',tspan,y0,[],DeltaCav,DeltaFinE,Theta);
        
        subplot(2,3,5),plot(y(:,5),y(:,6))
      % plot(t,y(:,6))  
       


%--------------------------------------------------------------------------