#include
#include
#define pi 3.1415926

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

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

	/* Initial Condition */
	b = 0.0;	


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

	/* Iteration */
	for(i = 0; i < iteration; i++){
		V =(atan(sqrt(b/(1-b))) + pi/4) / (sqrt(1-b));
		b = b + db;

		/* Output of the Results */
		
		fprintf(fp,"%6.5lf\t%6.5lf\n",b,V);
		}
	
	fclose(fp);
}