# 
restart;
phi:=a->(a^2+2)/(2*a);
phi(1);
phi(3/2);
phi(17/12);
a1:=phi(577/408);
a1^2-2;

#