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

int main(void)
{
	int n=0,m;
	float ans[10000];
	for(;;)
	{
		float a,b,c,d;
		scanf("%f %f",&a,&b);
		if(a==0&&b==0)
		{
			break;
		}
		c=a/2*a/2+b*b;
		d=sqrt(c);
		ans[n]=a*d/2*4+a*a;
		n++;
	}
	for(m=0;m<n;m++)
	{
		printf("%f\n",ans[m]);
	}
	return 0;
}