/**
 * X方向の変化率を与える関数
 */
double fX(double x, double y, double z, double s, double b, double r){
	return -s*x + s*y;
}

/**
 * Y方向の変化率を与える関数
 */
double fY(double x, double y, double z, double s, double b, double r){
	return -x*z + r*x - y;
}

/**
 * Z方向の変化率を与える関数
 */
double fZ(double x, double y, double z, double s, double b, double r){
	return x*y - b*z;
}