#include
#include

double fun(double t,double y);
double lambda;

void main()
{double t0,tk,y0,dt,bl,ya,t,y,yze,ypom,ble;
FILE *plik2=fopen("wynik2.xls","w");

printf("podaj lambda\n");
scanf("%lf",&lambda);
printf("podaj t0,y0,tk\n");
scanf("%lf%lf%lf",&t0,&y0,&tk);
printf("podaj dt\n");
scanf("%lf",&dt);
//wp
t=t0;
//y=y0;
yze=y0;
//całkujemy
while(t<=tk)
{
//y=y+dt*fun(t,y);
ypom=yze+0.5*dt*fun(t,yze);
yze=yze+dt*fun(t+0.5*dt,ypom);
ya=y0*exp(lambda*(t-t0));
t=t+dt;
//bl=fabs(ya-y);
ble=fabs(ya-yze);
printf("t=%lf bl=%lf\n",t,ble);
fprintf(plik2,"%lf %lf\n",t,ble);
}
}
double fun(double t,double y)
{double f;
f=lambda*y;
return f;
}
double rzut(double t,double v)
{double f;
f=g-c/m*v*v;
return f;
}