Public Class Form1
Dim g As Double = -9.81 'Gravity
Function Find_Verical_Velocity(ByVal Vy As Double, ByVal mass As Double, ByVal drag_y As Double) As Double
If Vy > 0 Then
Vy = Vy + ((g - drag_y / mass) * 0.01) 'The equation that gives the velocity is V=Vo-gt
Else
Vy = Vy + ((g + drag_y / mass) * 0.01)
End If
Return (Vy) 'Returns the vertical velocity
End Function
Function Find_Horizontal_Velocity(ByVal Vx As Double, ByVal mass As Double, ByVal drag_x As Double) As Double
Vx = Vx + (-(drag_x / mass) * 0.01)
Return (Vx)
End Function
Function Find_Distance(ByVal Vx As Double) As Double
Static Distance As Double = 0
Distance = Distance + (Vx * 0.01)
Return (Distance)
End Function
Function Find_Height(ByVal Vy As Double) As Double
Static Height As Double = 0
Height = Height + (Vy * 0.01)
Return (Height)
End Function
Private Sub Button1_Click(ByVal sender As System.Object, ByVal e As System.EventArgs) Handles Button1.Click
Dim mass As Double = CInt(TextBox1.Text)
Dim angle As Double = CInt(TextBox2.Text)
Dim initial_velocity As Double = CInt(TextBox3.Text) 'Given initial velocity
Dim density As Double = CInt(TextBox4.Text)
Dim Height As Double = 0.0
Dim time As Double = 0.0
Dim Distance As Double = 0.0
Dim Vx As Double = 0.0 'Horizontal velocity component
Dim Vy As Double = 0.0 'Vertical velocity component
Dim rad As Double
Dim new_velocity As Double
Dim Cd As Double = 0.47 'This is the drag coefficient and we assume that is constant
Dim A As Double = Math.PI * 0.05 ^ 2 'This is the cross-sectional area of the ball
Dim drag As Double
Dim drag_x As Double 'The horizontal drag component
Dim drag_y As Double 'The vertical drag component
Dim ITEMS As ListViewItem
If angle > 180 Then
MsgBox("PLEASE ENTER A SMALLER ANGLE")
ElseIf angle < 0 Then
MsgBox("PLEASE ENTER A POSITIVE ANGLE")
ElseIf mass < 0 Then
MsgBox("PLEASE ENTER A POSITIVE MASS")
Else
rad = (angle / 180) * Math.PI 'Calculate the angle in rads
Vx = initial_velocity * Math.Cos(rad)
Vy = initial_velocity * Math.Sin(rad)
Do
time = time + 0.01
Vy = Find_Verical_Velocity(Vy, mass, drag_y)
Vx = Find_Horizontal_Velocity(Vx, mass, drag_x)
rad = Math.Atan(Vy / Vx) 'Because the angle changes during the flihgt it calculates the angle in different points
new_velocity = Math.Sqrt((Vy ^ 2) + (Vx ^ 2))
drag = 0.5 * Cd * density * new_velocity ^ 2 * A
drag_x = drag * Math.Cos(rad)
drag_y = drag * Math.Sin(rad)
angle = (rad * 180) / Math.PI
Distance = Find_Distance(Vx)
Height = Find_Height(Vy)
ITEMS = New ListViewItem(Math.Round(time, 4))
ITEMS.SubItems.Add(Math.Round(Distance, 4))
ITEMS.SubItems.Add(Math.Round(Height, 4))
ITEMS.SubItems.Add(Math.Round(angle, 4))
ITEMS.SubItems.Add(Math.Round(new_velocity, 4))
ITEMS.SubItems.Add(Math.Round(Vy, 4))
ITEMS.SubItems.Add(Math.Round(Vx, 4))
ListView1.Items.Add(ITEMS)
Loop Until (Height < 0)
End If
End Sub
End Class
Comment