Algoritmo de Euler para ecuación diferencial de primer orden

c:
/* Euler  algorithm for first order differential equation */
/* dy/dt = -y(t), y(0)=1, 0<=t<=4, start h=0.5*/
#include <stdio.h>
 
#define dist 0.5                /* stepsize in t */
#define MAX 4.0                 /* max for t */
 
FILE *output;                   /* internal filename */

main()
{
double t, y;
int j;

output=fopen("euler.dat", "w"); /* external filename */

y=1;                            /* initial condition */
fprintf(output, "0\t%f\n", y);

for (j=1;dist*j<=MAX;j++)       /* the time loop */
{
   t=j*dist;
   y-=dist*y;
 
   fprintf (output, "%f\t%f\n", t, y);
}

fclose(output);
}
¿Ya le viste algún error? Dejanos tu correción ;-)

Antes de comentar: Gran parte de los ejercicios propuestos no tienen librerías debido a que Wordpress elimina los tags HTML. Si sabes/tienes/conoces las librerías que hacen falta, déjalo en los comentarios.

Otro punto antes de comentar, Si vas a sugerir un segmento de código en algún lenguaje debes hacerlo así:

De esta manera el código sale coloreado.

Otro punto importante para muchos que sienten que se les ignora: Todos los comentarios los reviso y en su debido momento los apruebo, pero ojo con el último párrafo.

Para poner los símbolos de las librerías sin que desaparezcan, debes escribir los símbolos de mayor > y menor qué < con su entidad HTML correspondiente, así como el símbolo de &

Mis Algoritmos se reserva el derecho de alterar, publicar o no los comentarios así como cambiar estas reglas de uso.

Si estas de acuerdo, adelante puedes comentar :)

P.D. No le hago tareas a nadie, mejor hagan la mía :P