import std.stdio;
import std.conv;
import std.string;
import std.typecons;
import std.algorithm;
import std.array;
import std.range;
import std.math;
import std.regex : regex;
import std.container;

void main()
{
	while (1) {
		auto x = readln.chomp.to!int;
		auto h = readln.chomp.to!int;
		if (x + h == 0) break;
		writefln("%.8f", 2*x*sqrt((x/2.0)^^2+h^^2)+x^^2);
	}
}