evaluate sqrt(2*a/sqrt((sqrt(a²-b²)-sqrt(b²-y²)a/b)²+y²)-1)*p/b with p=1,a=2,b=1,y=1