Eular法 のプログラム

#include < stdio.h >
main()
{
double t, x, dt, dx, ta, tb;
double f(double ,double );
int i, n;
x = 0.0 ;
t = 0.0 ;
ta = 0.0 ;
tb = 100 ;
n = 1000;
dt = (tb - ta) / n ;
for ( i = 0; i < n ; i++){
dx =f(x, t) *dt ;
x = x+ dx ;
t = t + dt ;
printf(" %e %e \n", t, x);
}
}
double f( double x, double t){
double y, z, alpha, g;
alpha = 0.1;
g = - 9.8;
z = -alpha * x + g ;
return z;
}