for( int i=0; i<LOOPS; i++ ){
	double v_x = dxdt(x, y, z);
	double v_y = dydt(x, y, z);
	double v_z = dzdt(x, y, z);
	x += v_x * DT;
	y += v_y * DT;
	z += v_z * DT;
}