Положение равновесия частицы, находящейся в потенциальной яме

Dim XX, YY, Alfa as Double

 

Public Function F(X) ‘ Выражение для силы

F = 6000 * (-0.0005 * Sin(0.0015 * X) ^ 2 + 2 * 0.0015 * Sin(0.0015 * X) * Cos(0.0015 * X)) * (Exp(-0.0005 * X))

End Function

 

Private Sub Command1_Click()

XX = Shape1.Left

YY = Shape1.Top

Timer1.Enabled = True

End Sub

 

Private Sub Command2_Click()

Timer1.Enabled = False

End Sub

 

‘ Рисуется график зависимости потенциальной энергии от

‘ координаты

Private Sub Form_Activate()

Alfa = 1

For i = 0 To 9000

Y = 6000 * Sin(0.0015 * i) ^ 2 * (Exp(-0.0005 * i))

Picture1.Line (i, 5000 - Y)-(i + 3, 5000 - Y + 3)

Next i

End Sub

 

‘ Перемещение частицы мышкой

Private Sub Picture1_MouseMove(Button As Integer, Shift As Integer, X As Single, Y As Single)

If Button = 1 Then Shape1.Left = X

End Sub

 

‘ Движение частицы по направлению силы

Private Sub Timer1_Timer()

XX = XX - Alfa * F(XX)

Shape1.Left = XX

End Sub