f := (x,y) -> sqrt(x-y);