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

//int main(int argc, char** argv)
int main()
{
    int step;
	double m=0,P=0;
	double initialp;
	double kplus=10000,kminus=0.00000002,kn=0.00018;

	int nc=2;
	double deltat=1;
	printf("Please input total protein concentration(uM/ml):\n");
        scanf("%lf",&initialp);
        printf("Please input reaction rates for fiber nucleation, elongation and fragmentation:\n");
 	printf("kn=");
        scanf("%lf",&kn);
	printf("k+=");
    scanf("%lf",&kplus);
	printf("k-=");
    scanf("%lf",&kminus);
//        initialp=atof(argv[1]);
//        kn=atof(argv[2]);
//        kplus=atof(argv[3]);
//        kminus=atof(argv[4]);

        printf("Time(h)    M(t)        P(t)      M(t)/P(t)\n");


	m=initialp;
	for(step=0;step<=3600*100;step++)
	{
		P=P+(kminus*(initialp-m-(2*nc-1)*P)+kn*m*m)*deltat;
		m=m-((2*kplus*m-(nc-1)*nc*kminus)*P+nc*kn*m*m)*deltat;
		if((int)(step/3600)*3600==step) printf("%.2f, %.3e, %.3e, %.3e\n", deltat*step/3600, 1-m/initialp, P/initialp, (initialp-m)/P);
	}
        
//      fitting with experimental data of M(t) vs. t:
 



}