// 99-Mthod for Uncertain Differential Equation
// Written by Microsoft Visual C++
// Copyright by UTLab @ Tsinghua University
// http://orsc.edu.cn/UTLab

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


int main( void )
{
   int i, j;
   double x, alpha, phiinverse, step=0.001, t;
  
  
   FILE *fp;
   fp=fopen("data.tex","w");

   alpha=0.01;
   for(j=1; j<=99; j++) {
	   alpha=0.01*j;
	   phiinverse=sqrt(3)/3.14159*log(alpha/(1-alpha));
	   x=1; //initial state X_0
	   t=1; //terminal time X_t
	   for(i=1;i<=t/step; i++){
		   x=x+x*step+x*phiinverse*step;
	   }	          
	   fprintf(fp,"%1.3f %1.3f ",x,alpha);
   }
  
   fclose(fp);
   return 1;
}