Pasted

#include <stdio.h>
#include <math.h>
 
#define dt 1e-2
 
#define x0 0
#define t0 0
#define tf 10
 
#define N ((tf-t0)/dt)
 
double dxdt(double x, double t)
{
    return -x*x*x + sin(t);
}
 
int main(void)
{
    double x = x0, xn, t = t0, tn;
    int i;
 
    FILE *h = fopen("./euler.dat","w+");
 
    for (i = 0; i < N; ++i)
    {
        xn  = x + dxdt(x,t)*dt;
        tn = t + dt;
        fprintf(h, "%d %f %f\n", i, t, x);
        x = xn;
        t = tn;
    }
 
    fclose(h);
    return 0;
}