#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <time.h> 


int main(void)
{
        FILE *fp;
        long i;
        float b,N0,t0,N,t,b1,u,deltat;
        srand(time(NULL));
        b=2;
        N0=0;
        t0=0;
        deltat=0.01
        fp = fopen("linearbirth.csv", "w");
        fprintf(fp," %g ,  %g  \n", N0,t0);
        for (i = 1; i < 500; i++)
        {
                N = N0;
                t = t0;
                b1 = b*deltat;
                u = rand()/(1.0+RAND_MAX);
                N0=N+(u<b1?1:0);
                t0=t+deltat;
                fprintf(fp," %g ,  %g  \n", N0,t0);
        }
        fclose(fp);
      return 0;
}