Метод Эйлера: численное решение начальных задач для дифференциального уравнения - C++ - C++ - Ответ 10475610

#include "stdafx.h"
#include <iostream>
#include <cmath>
#include <fstream>
 
double f_analit (double x) {
    return (exp(x) - exp(-x)) / 2.;
}
 
double f(double x, double y) {
    return -y + exp(x);
}
 
int main()
{
    double r = 0.001;
    double T = 1.0;
    double y = 0., x = y;
 
    std::cout.precision(6);
    std::cout.setf(std::ios::fixed);
    
    std::cout << " Analytical solution in point T : " << std::endl;
    std::cout << f_analit(T) << std::endl;
    std::ofstream  out("out.txt");
    do {
        x += r;
        y += r*f(x, y);
        out << x << " " << y;
        out << std::endl;
    } while (x < T);
    out.close();
    std::cout << " Numerical solution in point T : " << std::endl;
    std::cout << " x = " << x << "\t" << " y = " << y << std::endl;
    system("pause");
}