// X方向の時間変化率を与える関数
double dxdt(double x, double y, double z){
	return -S*x + S*y;
}

// Y方向の時間変化率を与える関数
double dydt(double x, double y, double z){
	return -x*z + R*x - y;
}

// Z方向の時間変化率を与える関数
double dzdt(double x, double y, double z){
	return x*y - B*z;
}