#include <iostream>
#include <fstream>
#include <cmath>
#include <valarray>

using namespace std;

int main() {
    double T, dt;

    cout << "Stabilita: reseni diferencialni rovnice metodou 1. radu" << endl;

    cout << "Zadejte maximalni cas T: "; cin >> T;
    cout << "Zadejte krok dt:         "; cin >> dt;

    int N = T / dt;
    cout << "Pocet kroku:             " << N << endl;

    valarray<double> v(N);
    valarray<double> v_exact(N);

    v[0] = 1;  // Pocatecni podminka v(0) = v0;
    v_exact[0] = 1;
    for(int i = 1 ; i < N ; ++i){
        v[i]        = v[i-1]*(1-dt);
        v_exact[i]  = exp(-i*dt);
    }

    ofstream file;
    file.open("pad_meteoritu1.txt");
    for(int i = 0; i < N; ++i){
        file << i*dt << " " << v[i]
             << " " << v_exact[i] << "\n";
    }
    file.close();

    return 0;
}