Algoritmo para calcular el azimut y dip de los ejes de un elipsoide de influencia

ALGORITMO PARA CALCULAR EL AZIMUT Y DIP DE LOS EJES DE UN ELIPSOIDE DE INFLUENCIA

DATOS DE ENTRADA
Alfa = azimut del eje principal a1
Beta= dip del eje principal a1
Teta= ladeo del ejesecundario a2
a1 a3
a2teta alfa Nortebeta

ALGORITMO (SE PROGRAMO EN VB 6.0)

Dim alfa As Double
Dim beta As Double
Dim teta AsDouble
Dim w, z, m, n, p, r, s, t, x, y, pi As Double
Private Sub form_load()

Text1.Text = ""
Text2.Text = ""
Text3.Text = ""
Text4.Text = ""
Text5.Text = ""
Text6.Text = ""
Text7.Text = ""End Sub
Private Sub calcular_Click()
alfa = Val(Text1.Text)
beta = Val(Text2.Text)
teta = Val(Text7.Text)
pi = 3.14159265358979
m = 90 + alfa + (Atn((Sin(beta * pi / 180)) * (Tan(teta * pi /180)))) * 180 / pi
n = 270 + alfa + (Atn((Sin(beta * pi / 180)) * (Tan(teta * pi / 180)))) * 180 / pi
x = (Cos(beta * pi / 180)) * (Sin(teta * pi / 180))
If x < 1 And x > -1 Then
p = (Atn(x /Sqr(-x * x + 1))) * 180 / pi
Else
If x = 1 Then
p = 90
End If
If x = -1 Then
p = -90
End If
End If
y = ((Cos(beta * pi / 180)) * (Cos(teta * pi / 180)))
If y < 1 And y > -1 Then
t =(Atn(y / Sqr(-y * y + 1))) * 180 / pi
Else
If y = 1 Then
t = 90
End If
If y = -1 Then
t = -90
End If
End If
If alfa >= 0 And alfa <= 360 Then
If teta >= 0 And teta <= 360 ThenIf beta >= -90 And beta <= 90 Then
Text1.Text = alfa
Text2.Text = beta
Text4.Text = p
Text6.Text = t
If teta = 0 Or teta = 360 Then
If p = 90 Or p = -90 Then
Text3.Text = 0
Else
w = 90...