Computation of the stationary distribution of the queue size in an M/G/1 queueing system with variable service rate