std::cout << std::fixed << std::setprecision(std::numeric_limits<double>::digits10 + 1);

//t=0からはじめる．

double t = 0.0;

//常微分方程式の初期条件を設定．

double y = y_init;

double C2 = 0.0;

//はじめの値を表示. 「時刻t シミュレーション結果 相対誤差の対数」の順で出力．

C2 = 1.0 - beta / y_init;

std::cout << t << " " << y_init << " " << std::log10(std::abs(y - beta / (1.0 - C2 * std::exp(- alpha * t))) / std::abs(beta / (1.0 - C2 * std::exp(- alpha * t)))) << std::endl;

//漸化式を解くよ．T_limit秒まで解く．

for(int i=1; t<T_limit; i++){

//dt秒後のyの値を求める．

y = y + dt * func(y);

//時間tをdt秒進める．

t = i * dt;

//t秒後のときのyの値を表示する．

if(i%INTV == 0){

std::cout << t << " " << y << " "<< std::log10(std::abs(y - beta / (1.0 - C2 * std::exp(- alpha * t))) / std::abs(beta / (1.0 - C2 * std::exp(- alpha * t)))) << std::endl;

}

}