B=B-(B-atan((z+N*(e^2).*sin(B))./sqrt(x.^2+y.^2)))./(1-N*(e^2).*cos(B).*(x.^2+y.^2)./(x.^2+y.^2+(z+N.*sin(B)*(e^2)))); % 这句有问题，有一些不必要的点（已修改）
