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

int main(int argc, char *argv[]) 
{
	double a, b, c, d, x1,y1, x2,y2;
	
	a= 1.0;
	b= 3.0;
	c= 1.0;
	
	d= b*b - 4.0*a*c;
	if (d > 0.0) {
		x1= (-b + sqrt(d))/(2.0*a);
		x2= (-b - sqrt(d))/(2.0*a);
		
		y1= a*x1*x1 + b*x1 + c;
		y2= a*x2*x2 + b*x2 + c;
		
		printf("x1= %.6lf, y1= %.26le\n", x1, y1);
		printf("x2= %.6lf, y2= %.26le\n", x2, y2);
	}
	else if (d < 0.0) {
		//
		
	}
	else {// d == 0.0
	   //
	   
	}
	
	return 0;
}// end of main()