cir_model.cpp

#include <iostream>
#include <random>
#include <vector>
using namespace std;

vector<double> generate_random_sequence(int n)
{
    random_device rd;
    mt19937 gen(rd());
    normal_distribution<double> dist(0.0, 1.0);

    vector<double> Z(n, 0.0);
    for (int i = 0; i < n; ++i)
    {
        double z = dist(gen);
        Z[i] = z;
    }
    return Z;
}

vector<double> short_rate(int n, double r0, double T, double volatility, double alpha, double b)
{
    vector<double> Z = generate_random_sequence(n);
    vector<double> r(n + 1);
    double dt = T / n;

    r[0] = r0;

    for (int i = 0; i < n; ++i)
    {
        r[i + 1] = r[i] + alpha * (b - r[i]) * dt + volatility * sqrt(max(r[i], 0.0) * dt) * Z[i];
    }
    return r;
}