Проблема с физической моделью

#include <SFML/Graphics.hpp>
#include <cmath>
#include "imgui.h"
#include "imgui-SFML.h"
#include <vector>
 
using namespace sf;
 
int main()
{
    const int windowWidth = 800;
    const int windowHeight = 600;
    bool simulation = false;
 
    RenderWindow window(VideoMode(windowWidth, windowHeight), "Mathematical Pendulum");
 
    ImGui::SFML::Init(window);
 
    const float g = 9.81;      // ускорение свободного падения
    const float l = 200;       // длина маятника
    float damping = 0.02;      // коэффициент затухания
    float angle = 180;         // начальный угол отклонения
    float angularVelocity = 0; // начальная скорость
 
    CircleShape pendulum(20);
    pendulum.setOrigin(pendulum.getRadius(), 0);
    pendulum.setPosition(windowWidth / 2, windowHeight / 2 + l);
                                                                          //создания маятника и нити
    VertexArray line(Lines, 2);
    line[0].position = Vector2f(windowWidth / 2, windowHeight / 5);
 
    Clock deltaClock, clock;
 
    std::vector<float> angleValues; // Вектор для хранения значений угла
 
    Sprite background;
    Texture back;
    back.loadFromFile("background.png");
    background.setTexture(back);
    while (window.isOpen())
    {
        float deltaTime = clock.getElapsedTime().asMicroseconds();
        clock.restart();
        deltaTime /= 40000;
        Event event;
        while (window.pollEvent(event))
        {
            ImGui::SFML::ProcessEvent(event);
            if (event.type == Event::Closed)
                window.close();
        }
 
        ImGui::SFML::Update(window, deltaClock.restart());
 
        ImGui::Begin("Settings");
 
        if (ImGui::SliderFloat("Angle", &angle, 0, 180))
        {
            angularVelocity = 0;
        }
 
        ImGui::SliderFloat("Damping", &damping, 0, 0.02);
 
        if (ImGui::Button("Start"))
        {
            simulation = true;
        }
 
        ImGui::SameLine();
 
        if (ImGui::Button("Stop"))
        {
            simulation = false;
        }
        ImGui::SameLine();
        if (simulation)
        {
            float angularAcceleration = -g / l * sin(angle * 3.14 / 180) - damping * angularVelocity;    //физ модель маятника
            angularVelocity += angularAcceleration * deltaTime;
            angle += angularVelocity * deltaTime * 180 / 3.14;
            angleValues.push_back(angle);
        }
        if (ImGui::Button("Reset"))
        {
            simulation = false;
            angle = 180;
            angularVelocity = 0;
            angleValues.clear();
        }
 
        ImGui::End();
 
       
        
        float x = l * sin(angle * 3.14 / 180);
        float y = l * cos(angle * 3.14 / 180);
        line[1].position = Vector2f(windowWidth / 2 + x, windowHeight / 5 + y);
        line[1].color = Color::White;
 
        pendulum.setPosition(line[1].position.x, line[1].position.y - 20);
 
        window.clear();
        window.draw(background);
        window.draw(pendulum);
        window.draw(line);
 
        // Рисуем график зависимости угла от времени
        ImGui::Begin("Graph");
        ImGui::PlotLines("Angle and Time", angleValues.data(), angleValues.size(), 0, nullptr, -180, 180, ImVec2(300, 200));
        ImGui::End();
 
        ImGui::SFML::Render(window);
       
        window.display();
    }
 
    ImGui::SFML::Shutdown();
 
    return 0;
}