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

#define G0 40
#define GI 10
#define DT 0.5
#define ENDTIM 10

int main() {

  double time=0;
  int i=0;
  double beta[] = {0.001, 0.01, 0.1,0.2,0.5,1.0, 10.0, 100.0, 200.0};
  FILE *fp=NULL;

  fp=fopen("output.k","w");

  if( fp==NULL ) {
    fprintf(stderr," Could not open file to write !");
    exit(-1);
  }

  fprintf(fp,"*KEYWORD\n");
  fprintf(fp,"*DEFINE_TABLE\n1000\n");

  for(i=0; i<9; i++) {
     fprintf(fp,"%20.4f\n", beta[i]);
  }

  for( i=0;i<9;i++) {
    fprintf(fp, "*DEFINE_CURVE_TITLE\n");
    fprintf(fp, "beta = %10.4f, G0=%10.4f, GI=%10.4f\n", beta[i], G0, GI);
    fprintf(fp, "%10d\n", i+1);
    for(time=0; time< ENDTIM; time+=DT) {
       fprintf(fp, "%20.4f, %20.4f\n", time, GI+(G0-GI)*exp(-beta[i]*time));
    }
  }

  fprintf(fp,"*END\n");

  fclose(fp);

}