EDO decaimento radioativo

/ método de euler para EDOs:
// f(x+dx) = f(x) + dxF(f,x) (+ condição inicial)

// decaimento radioativo: dN/dt=-(lambda)*N
// solução exata: N(t)=N(0)exp(-(lambda)*t)

#include <stdio.h>

#define dt 0.01
#define lambda 0.05

int main(void)
{
  int i=0;
  double N,t;

  N=10000;

  while (N>0.1)
    {
      N=N-dt*lambda*N;
      t=++i*dt; // note que não foi usado i++! i++: soma 1 ao i e dá o resultado de t, ++i: dá o resultado de t e no final do loop soma 1 ao i
      printf("%f\t %f\n",t,N);
    }

  return 0;
}