// Zadanie 3c
#include
#include
#define MAXN 10
double ro,g,m,D,Cd,a,pi,V0;
double tp=0;
void fun(double t, double *Y, double *F);
void vrk4( double x0, double y0[], double h, int n, void (*fun)(double, double*, double*), double y1[] );
void main()
{
double Y[4],Y1[4],t,h,p;
double L,H;
int i,n,j;
FILE *pl;
h=0.01;
n=4;
pi=3.14;
g=9.81;
V0=100;
a=3.14/9;
m=10;
D=0.1;
ro=1;
Cd=0.4;

pl=fopen("dane.txt","wt");
p=10;

for(j=0;j<=40;j++)
{
a=p*pi/180;
L=-1.;
H=0.;
Y[0]=0;
Y[1]=V0*cos(a);
Y[2]=0;
Y[3]=V0*sin(a);
t=tp;

for (i=0;i<=5000;i++)
{
t=tp+i*h;

vrk4(t,Y,h,n,fun,Y1);
// fprintf(pl," %lf %lf\n ", Y[0],Y[2]);
// printf("t=%lf, x=%lg, z=%lg\n", t, Y[0], Y[2]);

if(Y[2]<0&&L<0)L=Y[0];
if(Y[2]>H)H=Y[2];
Y[0]=Y1[0];
Y[1]=Y1[1];
Y[2]=Y1[2];
Y[3]=Y1[3];
// if(Y[2]<0)break; //zatrzymaj po upadku
}
printf("p=%lf stopni = %lf radianow, L=%lf, Hmax=%lf\n",p,a,L,H);
p=p+1.;
}
fclose(pl);
}