I added the corrected code $"$ p$\_$state$($i$)=($c$-$i$)/($c $-($c $-$ i$))$;$"$ but I am still gettting the same error message: $\%$ Part $1$: Define parameters

c $=2$; $\%$ Number of servers

lambda $=2$; $\%$ Average Arrival Rate per minute

mu $=0\mathrm{.}1$; $\%$ Service rate

Nwait $=15$; $\%$ Number of waiting customers

$\%$ Part $2$: Call the MMCQ function and retrieve results

$[$Ws$,$ Wq$,$ c$\_$util, p$\_$drop, p$\_$state$]=$ MMCQ$(1,0\mathrm{.}5,4,5)$;

$\%$ Part $3$: Define a relative error function

rel$\_$error $=$ @$($x$,$ y$)$ abs$($x $-$ y$)/$ y;

$\%$ Part $4$: Validate and print results

$\%$ p$\_$state$($i$)=($c$-$i$)/($c $-($c $-$ i$))$

fprintf$("$Ws: $\%6\mathrm{.}3$f

$",$ Ws$)$;

fprintf$("$Wq: $\%6\mathrm{.}3$f

$",$ Wq$)$;

fprintf$("$c$\_$util: $\%6\mathrm{.}3$f

$",$ c$\_$util$)$;

fprintf$("$p$\_$drop: $\%8\mathrm{.}5$f

$",$ p$\_$drop$)$;

fprintf$("$State Probabilities: $")$;

for i $=1$:Nwait $+1$

i$=0$:$1$:$10$;

p$\_$state$($i$)=($c$-$i$)/($c $-($c $-$ i$))$;

fprintf$("$p$\_$state$(\%$d$)$: $\%$f$,",$ i$,$ p$\_$state$($i$))$;

end

fprintf$("$

$")$;

$\%$ Part $5$: Validate and print "PASS" or "FAIL" for each metric

if rel$\_$error$($Ws$,2\mathrm{.}1557)<0\mathrm{.}005$

fprintf$("$Ws $\%6\mathrm{.}3$f PASS

$",$ Ws$)$;

else

fprintf$("$Ws $\%6\mathrm{.}3$f FAIL

$",$ Ws$)$;

end

if rel$\_$error$($Wq$,0\mathrm{.}1557)<0\mathrm{.}005$

fprintf$("$Wq$=\%6\mathrm{.}3$f PASS

$",$ Wq$)$;

else

fprintf$("$Wq$=\%6\mathrm{.}3$f FAIL

$",$ Wq$)$;

end

if rel$\_$error$($c$\_$util, $0\mathrm{.}4986)<0\mathrm{.}005$

fprintf$("$c$\_$util $\%6\mathrm{.}3$f PASS

$",$ cutil$)$;

else

fprintf$("$c$\_$util $\%6\mathrm{.}3$f FAIL

$",$ c$\_$util$)$;

end

if rel$\_$error$($p$\_$drop, $0\mathrm{.}00272)<0\mathrm{.}005$

fprintf$("$p$\_$drop $=\%8\mathrm{.}5$f PASS

$",$ p$\_$drop$)$;

else

fprintf$("$p$\_$drop $=\%8\mathrm{.}5$f FAIL

$",$ p$\_$drop$)$;

end

$\%$ Part $6$: Corrected p$\_$state initialization

p$\_$state $=[-0\mathrm{.}010609,-0\mathrm{.}021217,-0\mathrm{.}021217,-0\mathrm{.}014145,-0\mathrm{.}0070725,-0\mathrm{.}05658,0\mathrm{.}0,0\mathrm{.}0,0\mathrm{.}0,0\mathrm{.}01]$;

$\%$ Part $7$: Check if the sum of probabilities is close to $1$

if abs$($sum$($p$\_$state$)-1)<1$e$-5$

disp$("$The probabilities sum to $1.$ Test PASSED."$)$;

else

disp$("$The probabilities do not sum to $1.$ Test FAILED."$)$;

end

$\%$ Part $8$: Assuming $10$ states

num$\_$states $=10$;

fprintf$("$The probability of being in state $\%$d is $\%$f$.$

$",$ Nwait, p$\_$state$($Nwait $+1))$;

disp$("$Extra Tests PASSED."$)$;function $[$Ws$,$ Wq$,$ c$\_$util, p$\_$drop, p$\_$state$]=$ MMCQ$($lambda$,$ mu$,$ c$,$ Nwait$)$

$\%$ INPUTS:

$\%$ lambda $=$ arrival rate

$\%$ mu $=$ service rate

$\%$ Nwait $=$ Number of queue elements, not including servers

$\%$ c $=$ Number of servers

$\%$ OUTPUTS:

$\%$ Ws $=$ Average wait in the system

$\%$ Wq $=$ Average wait in the queue

$\%$ c$\_$util $=$ Average percentage of servers busy

$\%$ p$\_$drop $=$ probability the queue is full

$\%$ p$\_$state $=$ probability of each state

$\%$ INSERT YOUR CODE HERE

$\%$ Calculate traffic intensity

rho $=$ lambda $/($c $*$ mu$)$;

$\%\%$ Compute the Erlang B formula to calculate the probability of blocking

$\%$ p$\_$block $=$ erlangb$($c$,$ rho$)$;

$\%$ Calculate the utilization factor

c$\_$util $=$ rho $/$ c;

$\%$ Display the result

$\%$ Calculate p$0($probability of no customers in the system$)$

p$0\_$denom $=$ sum$(($c$*$rho$).^(0$:$($c$-1))./$ factorial$(0$:$($c$-1)))$;

p$0\_$denom $=$ p$0\_$denom $+(($c$*$rho$)^$c $/($factorial$($c$)*(1-$rho$)))*(1-($rho$^$Nwait $/(1-$ rho$)))$;

p$0=1/$ p$0\_$denom;

fprintf$(\text{'}$p$0$: $\%\mathrm{.}4$f

$\text{'},$ p$0)$;

$\%$ Calculate pN $($probability of having N customers in the system$)$

p$\_$state $=$ zeros$(1,$ Nwait$+1)$;

for i $=1$:Nwait

if i $<=$ c

p$\_$state$($i$+1)=($rho$^$i $/$ factorial$($i$))*$ p$0$;

else

p$\_$state$($i$+1)=($rho$^$i $/($factorial$($c$)*($c$^($i$-$c$))))*$ p$0$;

end

pN $=$ p$\_$state$($Nwait$+1)$;

$\%$ Calculate p$\_$drop $($probability of a customer being dropped$)$

p$\_$drop $=$ pN;

$\%$ Calculate c$\_$util $($utilization of the servers$)$

c$\_$util $=($rho$^$c $*$ p$0)/($factorial$($c$)*(1-$ rho$))*(1-$ rho$^($Nwait $-$ c $+1))$;

$\%$ Calculate Wq $($average waiting time in the queue$)$

Wq $=(($rho$^($c$+1))*$ p$0)/($factorial$($c$)*(1-$ rho$/$Nwait$)^2)$; $\%$ Average waiting time in the queue

$\%$ Print the result

fprintf$(\text{'}$Average waiting time in the queue $($Wq$)$: $\%\mathrm{.}2$f$\text{'},$ Wq$)$;

$\%$ Calculate Ws $($average time spent in the system$)$

Ws $=$ Wq $+(1/$ mu$)$;

$\%$ Print the results

fprintf$(\text{'}$p$0$: $\%\mathrm{.}4$f $\text{'},$ p$0)$; fprintf$(\text{'}$pN: $\%\mathrm{.}2$f$\text{'},$ p$0)$;

end

end