/*
 * dy/dt = y
 * y(0) = y_init
 *
 * 
 *
 * 
 */


#include 

//
const double y_init  = 1.0;
const double dt      = 0.01; //
const double T_limit = 10.0;

double func(double y){
    return y;
}

int main(){
    //t=0
    double t = 0.0;
    //
    double y = y_init;

    //RKk
    double k1, k2, k3, k4;
    k1 = k2 = k3 = k4 = 0.0;

    //.
    std::cout