from numpy import sqrt, cos, sin
import ellipsoide
import ausgaben as aus
import winkelumrechnungen as wu


re = ellipsoide.Ellipsoid.init_name("Bessel")
eta = lambda phi: sqrt(re.e_ ** 2 * cos(phi) ** 2)

A = wu.gms2rad([327,0,0])
phi = wu.gms2rad([35,0,0])
h = 3500
dA = eta(phi)**2 * h * sin(A)*cos(A) / re.N(phi)
print(aus.gms("dA", dA, 10))