I need a MATLAB program to simulate two classes M/G/1 Queueing System. Exact question is below:

Consider a M/G/1 queueing system model. Messages that arrive to the system belong to two classes: class-1 and class-2. Class 1 messages are shorter, requiring each a service time that is exponentially distributed with mean message service time equal to 1 msec. The class-1 message arrival rate is lambda(1) . Class 2 messages are longer, requiring each a service time that is equal to 2 msec for 60% of the messages and is equal to 10 msec for 40% of the messages. Class-2 message arrival rate is lambda(2). We observe that lambda(1) = 5lambda(2). Assume that the system designer evaluates the performance of the system by increasing the message arrival rate, in inducing a total (including both flow classes) traffic intensity rho that varies over values between 0.1 and 0.9, in increments of 0.1.

The service scheduling policy used is set to be equal to one of the following models:

a. Messages are served on a FIFO basis; no differentiation is made in selecting messages for service based on their class.

b. Class-1 messages are served at higher priority than class-2 messages. A non-preemptive service policy is used.

1. Calculate the average steady-state message waiting time under (a) and the corresponding values for each class under (b).

2. Simulation program for this system and include it in your solution. Then, use the simulation program, and draw the graphs representing the following metrics for each class of message type, by increasing the message arrival rate as noted above, showing the following:

a. The average message waiting time for each class of messages.

b. The standard deviation of the message waiting time for each class of messages.

c. The 95-percentile message waiting time (i.e., the waiting time value that is met by 95% of the served messages) for each class of messages.

d. The ratios of the class-2 to class-1 waiting times for the means, standard-deviations and 95-percentile.