void rk4(double y[], double dydx[], int n, double x, double h, double yout[],void (*derivs)(double x, double yt[], double dyt[])); // nはベクトル（配列）次数

double *vector(long size);
void free_vector(double *v);