求救！这个迭代算法怎么解？？

[rmn]

h=z/sin(B)-N(1-E^2)
b=arctan{z(N+H)/sqrt[(x^2+y^2)N(1-E^2)+H]}

x,y,z为已知变量 （double）
N=A/sqrt[(1-E^2sin(B)sin(B)]
A=6378140,E^2=0.006694385

求解h和b。 （double）

可用迭代算法解，初值该赋多少？能贴出其代码吗？
谢谢您！！