The paper presents the stability problem of control systems composed of a fractional-order PI controller and a inertial plant of a fractional order with time delay. Simple and efficient computational method for determining stability regions in the controller and plant parameters space is given. Knowledge of these regions permits tuning of the fractional-order PI controller. The method proposed is based on the classical D-partition method.