grad(exp( -i *magnitude (x,y) ) /magnitude (x,y))