1function [ percentileRTs ] = returnRT1( arrival, service, pers )
3% The following function needs to use the QMAM toolbox, which can be
8[ ql, res_alpha, wait_alpha, Smat ] = Q_CT_MAP_MAP_1( arrival.lambda0, arrival.lambda1,service.ST, service.St*service.tau_st );
10percentileRTs = returnPer( res_alpha, Smat, pers );