"Born approximation 3" x = (x1,x2,x3) y = (y1,y2,y3) check(dot(x-y,x-y) == dot(x,x) - 2 dot(x,y) + dot(y,y)) "ok" Run