#include <iostream>
#include <cmath>
#include <fstream>
#include <vector>

int main()
{
    // Constantes
    const double m = 2.6;  // Massa em g
    const double L = 1.0;  // Comprimento em m
    const double R = 0.03; // Raio em m

    // Valores iniciais
    double theta0 = 0.05;  // Ângulo inicial em rad
    double w0 = 0.0;       // Velocidade angular inicial em rad/s
    double t0 = 0.0;       // Tempo inicial em s
    double h = 0.1;        // Passo de tempo em s

    // Listas para armazenar os resultados
    std::vector<double> t;
    std::vector<double> theta;
    std::vector<double> w;

    t.push_back(t0);
    theta.push_back(theta0);
    w.push_back(w0);

    // Loop de tempo
    while (t.back() < 100.0)
    {
        // Cálculo dos valores intermediários
        double k1x = h * w.back();
        double k1v = h * (-1.0 if w.back() < 0.0 else 1.0);

        // Atualização dos valores de θ e w
        theta.push_back(theta.back() + k1x);
        w.push_back(w.back() + k1v);

        // Atualização do tempo
        t.push_back(t.back() + h);
    }

return 0;
}