f(z):=z^2-%i*z=abs(z - %i);

u(x,y):=f(x+y*%i);

realpart(u(x,y));

solve(%,y);

realpart(u(0,1));

is(realpart(u(0,1)));

is(realpart(u(1,0)));

imagpart(u(0,1));

is(imagpart(u(0,1)));

imagepart(u(1,0));

false$