{VERSION 5 0 "IBM INTEL NT" "5.0" }
{USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 
1 0 0 0 0 1 }{CSTYLE "2D Math" -1 2 "Times" 0 1 0 0 0 0 0 0 2 0 0 0 0 
0 0 1 }{CSTYLE "2D Output" 2 20 "" 0 1 0 0 255 1 0 0 0 0 0 0 0 0 0 1 }
{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 
2 2 1 1 1 1 }1 1 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Warning" -1 7 1 
{CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 2 2 2 2 2 1 1 1 3 1 }1 1 0 
0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Maple Output" -1 11 1 {CSTYLE "" -1 
-1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }3 3 0 0 0 0 1 0 1 0 2 2 
0 1 }}
{SECT 0 {EXCHG {PARA 0 "" 0 "" {TEXT -1 2 "Q1" }{MPLTEXT 1 0 0 "" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart:" }}}{EXCHG {PARA 0 "
> " 0 "" {MPLTEXT 1 0 45 "int(x^3*exp(-x/36)/3!/36.^4,x=120..infinity)
;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+?*f)Hd!#5" }}}{EXCHG {PARA 0 
"> " 0 "" {MPLTEXT 1 0 29 "P:=(1-36*13.7/60*(z-1))^(-4):" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "1-eval(mtaylor(P,z,51),z=1);" }}
{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+g'R0\\\"!#5" }}}{EXCHG {PARA 0 "> \+
" 0 "" {MPLTEXT 1 0 24 "mu:=eval(diff(P,z),z=1);" }}{PARA 11 "" 1 "" 
{XPPMATH 20 "6#>%#muG$\"++++)G$!\")" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 36 "sqrt(eval(diff(P,z,z),z=1)+mu-mu^2);" }}{PARA 11 "" 
1 "" {XPPMATH 20 "6#$\"+n18T<!\")" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 
2 "Q2" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 
"restart: with(linalg):" }}{PARA 7 "" 1 "" {TEXT -1 80 "Warning, the p
rotected names norm and trace have been redefined and unprotected\n" }
}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 89 "A:=matrix([ [0,21,35,10,8]
,[9,0,23,18,20],[33,11,0,20,17],[41,0,27,0,4],[7,12,13,21,0]]):" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 49 "for i to 5 do A[i,i]:=-add(A
[i,j],j=1..5) end do:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "s:
=linsolve(transpose(A),[0,0,0,0,0]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "
6#>%\"sG-%'vectorG6#7',$*&#\"(C(*)>\"(TH:\"\"\"\"&%#_tG6#F.F.F.F/,$*&#
\"('[L>F-F.F/F.F.,$*&#\"(R2`\"F-F.F/F.F.,$*&#\"(56Z\"F-F.F/F.F." }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "s:=evalm(s/add(s[i],i=1..5))
; " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"sG-%'vectorG6#7'#\"'Ju\\\"(+
&>?#\"(TH:\"\"(+!y!)#\"'Vn'*\"(+!RS#\"'x'=#\"(+S:\"#\"'6r9\"'+y!)" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "exponential(A*.1,17/60.)[5,1
];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#^$$\"+BFo`=!#5$!*9Dp;#!#?" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "for i to 5 do A[3,i]:=0 end \+
do:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "exponential(A*.1,17/
60.)[5,1];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#^$$\"+_5\\#3\"!#5$\"+XJ2
hQ!#@" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 2 "Q3" }{MPLTEXT 1 0 0 "" }}
}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart:" }}}{EXCHG {PARA 0 
"> " 0 "" {MPLTEXT 1 0 25 "emi:=exp(-(lambda-mu)*t):" }}}{EXCHG {PARA 
0 "> " 0 "" {MPLTEXT 1 0 65 "r:=piecewise(mu=lambda,mu*t/(1+mu*t),mu*(
1-emi)/(lambda-mu*emi)):" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 67 "p:=piec
ewise(mu=lambda,1/(1+mu*t),(lambda-mu)*emi/(lambda-mu*emi)):" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 60 "PGF:=(r+(1-r)*p*z/(1-(1-p)*z
))^i*(p/(1-(1-p)*z))^(a/lambda):" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 47 "lambda:=1/3.: mu:=1/(2+35/60.): a:=2.3: i:=75: " }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "t:=9: eval(mtaylor(PGF,z,50)
,z=1);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+!>1uH#!#5" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "add(1./j,j=1..i)/mu;" }}{PARA 11 "
" 1 "" {XPPMATH 20 "6#$\"+QN=m7!\")" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 28 "sqrt(add(1./j^2,j=1..i))/mu;" }}{PARA 11 "" 1 "" 
{XPPMATH 20 "6#$\"+$y())*H$!\"*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 
1 0 42 "w[0]:=0:  w[1]:=ln(mu/(mu-lambda))/lambda:" }}}{EXCHG {PARA 0 
"> " 0 "" {MPLTEXT 1 0 84 "for j to 74 do w[j+1]:=(1+mu/lambda)*w[j]-m
u/lambda*w[j-1]-1/j/lambda end do: w[75];" }}{PARA 11 "" 1 "" 
{XPPMATH 20 "6#$\"+vGC'e&!\")" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 2 "Q
4" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "rest
art:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "lambda:=n->.6*sqrt(
n):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "mu:=n->3*n/(1+n):   \+
" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "dbn:=n->product(mu(j)/l
ambda(j),j=1..n):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "S:=sum
(dbn(n),n=0..120);  " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"SG$\"+sF$*
ya!\"&" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 46 "a:=n->evalf(1-sum
(dbn(k),k=0..n-1)/S):  a(30);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+P
-2K>!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 102 "lambda(30.)/(la
mbda(30.)+mu(30.))*lambda(31.)/(lambda(31.)+mu(31.))*lambda(32.)/(lamb
da(32.)+mu(32.));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+*)o\"*G:!#5" 
}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "exp(-(lambda(30.)+mu(30.)
)*22/60);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+TDjL5!#5" }}}{EXCHG 
{PARA 0 "" 0 "" {TEXT -1 2 "Q5" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "
> " 0 "" {MPLTEXT 1 0 80 "restart: with(stats): Ztail:=z->1-statevalf[
cdf,normald](z): c:=12: t:=2+24/60.:" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 42 "statevalf[cdf,normald[3.2,sqrt(c*t)]](-2);" }}{PARA 
11 "" 1 "" {XPPMATH 20 "6#$\"+u(=Gm\"!#5" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 29 "2*Ztail((5.3-3.2)/sqrt(c*t));" }}{PARA 11 "" 1 "" 
{XPPMATH 20 "6#$\"+oqmbp!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 
56 "Ztail((2+3.2)/sqrt(c*t))-Ztail((2-3.2+2*5.3)/sqrt(c*t));" }}{PARA 
11 "" 1 "" {XPPMATH 20 "6#$\"+oTfj7!#5" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 38 "evalf(2/Pi*arccos(sqrt(t/(3+37/60))));" }}{PARA 11 "
" 1 "" {XPPMATH 20 "6#$\"+Qk)*QR!#5" }}}{EXCHG {PARA 0 "" 0 "" {TEXT 
-1 2 "Q6" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 
22 "restart: with(linalg):" }}{PARA 7 "" 1 "" {TEXT -1 80 "Warning, th
e protected names norm and trace have been redefined and unprotected\n
" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 68 "A:=matrix([[11,19,15,4]
,[7,5,9,-1],[-12,-15,-15,-6],[-2,-7,-3,-7]]):" }}}{EXCHG {PARA 0 "> " 
0 "" {MPLTEXT 1 0 15 "eigenvalues(A);" }}{PARA 11 "" 1 "" {XPPMATH 20 
"6&\"\"$!\"'F#F$" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 10 "w:=[3,-
6]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 62 "eq:=seq(add(c[i]*w[i
]^k+d[i]*k*w[i]^(k-1),i=1..2)=A^k,k=0..3);" }}{PARA 11 "" 1 "" 
{XPPMATH 20 "6#>%#eqG6&/,&&%\"cG6#\"\"\"F+&F)6#\"\"#F+F+/,**&\"\"$F+F(
F+F+&%\"dGF*F+*&\"\"'F+F,F+!\"\"&F4F-F+%\"AG/,**&\"\"*F+F(F+F+*&F6F+F3
F+F+*&\"#OF+F,F+F+*&\"#7F+F8F+F7*$)F9F.F+/,**&\"#FF+F(F+F+*&FHF+F3F+F+
*&\"$;#F+F,F+F7*&\"$3\"F+F8F+F+*$)F9F2F+" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 34 "solve(\{eq\},\{c[1],c[2],d[1],d[2]\});" }}{PARA 11 "
" 1 "" {XPPMATH 20 "6#<&/&%\"cG6#\"\"#,*#\"\"(\"#F\"\"\"*&#\"\"%F,F-%
\"AGF-!\"\"*&#F-\"#\")F-*$)F1F(F-F-F-*&#F(\"$H(F-*$)F1\"\"$F-F-F-/&%\"
dGF',(#F(F=F-*&#F-F=F-F1F-F2*&F4F-F;F-F-/&F@6#F-,(#F0F=F2*&#F-\"\"*F-F
6F-F-*&F4F-F;F-F-/&F&FH,*#\"#?F,F-*&#F0F,F-F1F-F-*&#F-F5F-F6F-F2*&#F(F
:F-F;F-F2" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "assign(%); w:=
w*1.:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 71 "evalm(c[1]*tan(w[1
])+d[1]/cos(w[1])^2+c[2]*tan(w[2])+d[2]/cos(w[2])^2);" }}{PARA 11 "" 
1 "" {XPPMATH 20 "6#-%'matrixG6#7&7&$\"+#H.eo%!\"*$\"+([K4j\"!\")$\"+X
'eq>\"F-$!+W,U1aF*7&$\"*v#>pmF*$\"+1P$z;&F*$\"+4$evq#F*$!*IW&f&*F*7&$!
+e0JMCF*$!+-))y\\:F-$!+.Ux)p)F*$\"*-;9I$F-7&$\"+osc$H\"F*$!+k1\\lSF*$
\"*^xCt#F*$!+j\\V!e)!#5" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 2 "Q7" }
{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 60 "restart
:  c:=7:  lambda:=15.2:  mu:=60./13:  rho:=lambda/mu:" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 45 "G:=add(rho^i/i!,i=0..c-1)+rho^c/c!/
(1-rho/c):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 6 "rho/c;" }}
{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+1>w/Z!#5" }}}{EXCHG {PARA 0 "> " 
0 "" {MPLTEXT 1 0 7 "lambda;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"$_
\"!\"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "rho^4/4!/G*(lamb
da+4*mu);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+Vn%z2'!\"*" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 56 "(rho+rho^c/c!/G*rho/c/(1-rho
/c)^2)/lambda-1/mu;  %*3600;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\")hE
!R$!#5" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+'z&\\?7!\")" }}}{EXCHG 
{PARA 0 "" 0 "" {TEXT -1 2 "Q8" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "
> " 0 "" {MPLTEXT 1 0 8 "restart:" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 59 "pdsolve(z^2*diff(P(z,t),t)=z*diff(P(z,t),z)+P(z,t),P(
z,t));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"PG6$%\"zG%\"tG*&-%$_F1G
6#,&F(\"\"\"*&\"\"#!\"\"F'F0F.F.F'F1" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 22 "P(z,t):=(2*t+z^2)^2/z:" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 53 "simplify(z^2*diff(P(z,t),t)-z*diff(P(z,t),z)-P(z,t));
" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"\"!" }}}{EXCHG {PARA 0 "" 0 "" 
{TEXT -1 2 "Q9" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 49 "restart:  lambda:=7:  mu:=5.6:  N:=14:  t:=4/60.:" }}
}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 96 "p1:=(lambda+mu*exp(-(lambda
+mu)*t))/(lambda+mu): p2:=lambda*(1-exp(-(lambda+mu)*t))/(lambda+mu):
" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 62 "PGF:=(1-p1+p1*z)^12*(1-
p2+p2*z)^2: add(coeff(PGF,z,i),i=0..9);" }}{PARA 11 "" 1 "" {XPPMATH 
20 "6#$\"+1Dx/Y!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 119 "with
(linalg): \nM:=matrix([[0,0,0,0,0],[61.6,-82.6,21,0,0],[0,67.2,-81.2,1
4,0],\n[0,0,72.8,-79.8,7],[0,0,0,78.4,-78.4]]):" }}{PARA 7 "" 1 "" 
{TEXT -1 80 "Warning, the protected names norm and trace have been red
efined and unprotected\n" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 
"exponential(M,4/60.)[3,1];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+?(
\\fH)!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 64 "Pinf:=((mu+lamb
da*z)/(lambda+mu))^14:  coeff(Pinf,z,10)*(28+56);" }}{PARA 11 "" 1 "" 
{XPPMATH 20 "6#$\"+#=)y)=*!\"*" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 3 "
Q10" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "r
estart: with(linalg):" }}{PARA 7 "" 1 "" {TEXT -1 80 "Warning, the pro
tected names norm and trace have been redefined and unprotected\n" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 111 "A:=matrix([[0,0,0,0,0,0],[9
,0,23,18,20,0],[0,11,0,20,17,14],[41,0,27,0,4,28],[21,32,0,9,0,31],\n[
0,0,0,0,0,0]]):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 49 "for i to
 6 do A[i,i]:=-add(A[i,j],j=1..6) end do:" }}}{EXCHG {PARA 0 "> " 0 "
" {MPLTEXT 1 0 60 "exponential(A*.1,23/60.)[3,1]+exponential(A*.1,23/6
0.)[3,6];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+i(=+N'!#5" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 68 "F:=exponential(A*.1,t)[3,1]+exponen
tial(A*.1,t)[3,6]-.326*10^(-9)*t:" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 35 "mu:=int(t*diff(F,t),t=0..infinity);" }}{PARA 11 "" 1 
"" {XPPMATH 20 "6#>%#muG$\"+wP9;P!#5" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 44 "sqrt(int((t-mu)^2*diff(F,t),t=0..infinity));" }}
{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+_/ZOL!#5" }}}{EXCHG {PARA 0 "> " 
0 "" {MPLTEXT 1 0 66 "linsolve(submatrix(A,2..5,2..5),submatrix(-A,2..
5,[1,6]))[2,2]*1.;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+K4N*3'!#5" }
}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 3 "Q11" }{MPLTEXT 1 0 0 "" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 47 "restart:  with(stats):  d:=-
2.3/60:  c:=12/60.:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 90 "t1:=
23: t2:=161: t:=90: mu:=(3.2*(t2-t)-(t-t1))/(t2-t1): s:=sqrt(c*(t-t1)*
(t2-t)/(t2-t1)):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "1-state
valf[cdf,normald[mu,s]](0);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+\"o
#*zq'!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 92 "t1:=41: t2:=77:
 t:=60: mu:=(-1*(t2-t)-0.8*(t-t1))/(t2-t1): s:=sqrt(c*(t-t1)*(t2-t)/(t
2-t1)):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "1-statevalf[cdf,
normald[mu,s]](0);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+i7f@D!#5" }}
}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "mu:=-0.8+d*43: s:=sqrt(c*43
):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "1-statevalf[cdf,norma
ld[mu,s]](0);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+,7%*=?!#5" }}}}
{MARK "87" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 
1 1 }