#include
#include
#define pi 3.1415926

main(){
	/* Declaration */
	int i;
	double b, db, V,vg;
	int iteration;
	FILE* fp;

	/* Parameters */
	db = 0.001;
	iteration = 1000;
	

	/* Initial Condition */
	b = 0.0;	


	/* Output of Initial Value */
	fp = fopen("HE11_vg.txt","w");

	/* Iteration */
	for(i = 0; i < iteration; i++){
		V =(atan(sqrt(b/(1-b))) + pi/4) / (sqrt(1-b));
		vg = (0.01 * (1 - b)*sqrt(b) *V /(sqrt(b) * V +1) +1);
        b = b + db;
        
		/* Output of the Results */
		
		fprintf(fp,"%6.5lf\t%6.5lf\t%2.10lf\n",b,V,vg);
		}
	
	fclose(fp);
}