[b,S]=e721b

function [b,S]=e721b
%
%  The service time distribution:  phase(b,S)
%
b=zeros(1,6); b(1,1)=.7; b(1,5)=.3;
v=6.67778512*ones(1,4); v1=v(1,1:3);
P=diag(-v)+diag(v1,1); P1=zeros(4,2); P=[P P1];
Q=zeros(2,4); v=1.033237*ones(1,2); v1=v(1,1);
Q1=diag(-v)+diag(v1,1); Q=[Q Q1]; S=[P;Q];