Imports System

Public Module Module1

    Public Sub Main()

        Dim py = 1.0

        Dim px = 0.0

        Dim t = 0.0

        Dim rate = 0.06851

        Dim fuelmass = 3.0

        Dim tf = ((-fuelmass)/(-rate)) 'Final Time

        Dim th = 167.969 'Thrust

        Dim im = 5.865 'Initial mass

        Console.WriteLine("Interval:")

        Dim interval = Console.ReadLine()

        Dim lastang = 0.0

        Dim ang = 0.0

        Dim term1 = 0.0

        Dim term2 = 0.0

        Dim dy = 0.0

        Dim dx = 0.0

        Dim vy = 0.0

        Dim vx = 0.0

        Dim vxy = 0.0

        Console.WriteLine("Setting (1/2):")

        Dim Setting = Console.ReadLine()

        Do Until Math.Round(t) = Math.Round(tf)

            'vy=((th/(im-rate*t))*py-9.81)*t

            'vx=((th/(im-rate*t))*px)*t

            'vxy=Math.Sqrt((vx^2)+(vy^2))

                'ang = (t/tf)*90

            'ang = vxy/10

            ang = (dx/48000)*90

                py=(90-ang)/90

                px=ang/90

            dy=0.5*((th/(im-rate*t))*py-9.81)*(t^2)

            dx=0.5*((th/(im-rate*t))*px)*(t^2)

            lastang = ang

                t += interval

            If Setting = 1 Then

                Console.WriteLine("t="+Math.Round(t, 4, MidPointRounding.AwayFromZero).ToString()+", Angle="+Math.Round(ang, 2, MidpointRounding.AwayFromZero).ToString()+", py="+Math.Round(py, 2, MidpointRounding.AwayFromZero).ToString()+", px="+Math.Round(px, 2, MidpointRounding.AwayFromZero).ToString()+"")

            Else If Setting = 2 Then

                    Console.WriteLine("t="+Math.Round(t, 4, MidPointRounding.AwayFromZero).ToString()+", dy="+Math.Round(dy, 3, MidpointRounding.AwayFromZero).ToString()+", dx="+Math.Round(dx, 3, MidpointRounding.AwayFromZero).ToString()+"")

                    End If

                Loop

            If Setting = 1 Then

            t=tf

            ang = (t/tf)*90

            py=(90-ang)/90

            px=ang/90

            Console.WriteLine("tf="+t.ToString()+", Angle="+Math.Round(ang, 2, MidpointRounding.AwayFromZero).ToString()+", py="+Math.Round(py, 2, MidpointRounding.AwayFromZero).ToString()+", px="+Math.Round(px, 2, MidpointRounding.AwayFromZero).ToString()+"")

        Else If Setting = 2 Then

            t=tf

            ang = (t/tf)*90

            py=(90-ang)/90

            px=ang/90

                    dy=0.5*((th/(im-rate*t))*py-9.81)*(t^2)

            dx=0.5*((th/(im-rate*t))*px)*(t^2)

                    Console.WriteLine("t="+Math.Round(t, 4, MidPointRounding.AwayFromZero).ToString()+", dy="+Math.Round(dy, 3, MidpointRounding.AwayFromZero).ToString()+", dx="+Math.Round(dx, 3, MidpointRounding.AwayFromZero).ToString()+"")

                    End If

    End Sub

End Module