{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 35 "L:=11.3: f:=piecewise(t<1, t,2-t):" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 2 "a)" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "int(f,t=0..72/60.);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+++++o!#5" }}}{EXCHG {PARA 0 "" 0 " " {TEXT -1 2 "b)" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "M:=int(f*exp(t*u),t=0..2):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "P:=subs(u=L*(z-1),M):" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 49 "aux:=series(P,z,50): add(coeff(aux,z,i),i=0..12);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+R:LYg!#5" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 2 "c)" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "mu:=limit(diff(P,z),z=1);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#muG$\"++++I6!\")" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 "sqrt(limit(diff(P,z,z),z=1)+mu-mu^2);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+v_.3d!\"*" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 2 "Q2" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 " restart:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 84 "A:=matrix(5,5,[ 0,21,25,10,18,19,0,23,8,21,32,14,0,22,17,40,0,27,0,14,17,8,14,23,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 "" {TEXT -1 2 "a)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 46 "with(linalg): s:=linsolve (transpose(A),[0$5]):" }}{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 29 "S:=evalm(s/add(s[i],i=1..5)) ;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"SG-%'vectorG6#7'#\"(\"4$y\"\" ((R7n#\"'>d'*F+#\"(r2Q\"F+#\"(B!46F+#\"($zt9F+" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 2 "b)" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 48 "B:=evalm(A/10.): for i to 5 do B[3,i]:=0 end do:" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 27 "exponential(B,12/60.)[5,1]; " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#^$$\"+(RUw_\"!#5$\"+Kg%HR\"!#?" }} }{EXCHG {PARA 0 "" 0 "" {TEXT -1 2 "c)" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 57 "w:=linsolve(submatrix(A/10.,[1,3,4, 5],[1,3,4,5]),[-1$4]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"wG-%'vec torG6#7&$\"+5y/gv!#5$\"+10R8#)F+$\"+T<'oD*F+$\"+BKXu*)F+" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 52 "linsolve(submatrix(A/10.,[1,3,4,5], [1,3,4,5]),-2*w);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7&$\" +hQf_7!\"*$\"+u_Gw8F)$\"+Amkp:F)$\"+.S,E:F)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "sqrt(%[4]-w[4]^2);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+.4&))[)!#5" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 2 "Q3" } {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 70 "restart :r:=mu*(1-exp(-(lambda-mu)*t))/(lambda-mu*exp(-(lambda-mu)*t)):" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 67 "p:=(lambda-mu)*exp(-(lambda- mu)*t)/(lambda-mu*exp(-(lambda-mu)*t)):" }}}{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/l ambda):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 35 "lambda:=1./32: m u:=1./27: a:=0.12:" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 2 "a)" } {MPLTEXT 1 0 0 "" }{TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "i:=32: t:=72.:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 52 "a ux:=mtaylor(PGF,z,36): add(coeff(aux,z,i),i=0..27);" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#$\"+_4%f6&!#5" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 2 " b)" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "a* add(1/j,j=1..i)/mu;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+VC&\\J\"!\" )" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 54 "sqrt(a*add(1/j,j=1..i) /mu+a^2*add(1/j^2,j=1..i)/mu^2);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$ \"+
" 0 "" {MPLTEXT 1 0 54 "Digits:=10: w[0]:=
0: w[1]:=ln(mu/(mu-lambda))/lambda:" }}}{EXCHG {PARA 0 "> " 0 ""
{MPLTEXT 1 0 77 "for i to 31 do w[i+1]:=(1+mu/lambda)*w[i]-mu/lambda*w
[i-1]-1/lambda/i end do:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 6 "
w[32];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+d$fI0%!\"(" }}}{EXCHG
{PARA 0 "" 0 "" {TEXT -1 2 "Q4" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "
" 0 "" {TEXT -1 2 "a)" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 ""
{MPLTEXT 1 0 30 "restart: lambda:=n->3*n/(1+n):" }}}{EXCHG {PARA 0 "> \+
" 0 "" {MPLTEXT 1 0 19 "mu:=n->0.6*sqrt(n):" }}}{EXCHG {PARA 0 "> " 0
"" {MPLTEXT 1 0 51 "lomp:=n->product(lambda(k-1)/mu(k),k=2..n)/mu(1): \+
" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 51 "w(0):=0: w(1):=evalf(
add(lomp(n),n=1..300)); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>-%\"
wG6#\"\"\"$\"+%=vK\">!\"&" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0
97 "for k to 27 do w(k+1):=evalf((1+mu(k)/lambda(k))*w(k)-mu(k)/lambda
(k)*w(k-1)-1/lambda(k)) end do:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT
1 0 6 "w(28);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+!3bqL$!\"&" }}}
{EXCHG {PARA 0 "" 0 "" {TEXT -1 2 "b)" }{MPLTEXT 1 0 0 "" }}}{EXCHG
{PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "evalf(exp(-32/60.*(lambda(28)+mu(28
))));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+lY!Q#R!#6" }}}{EXCHG
{PARA 0 "" 0 "" {TEXT -1 2 "c)" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "
> " 0 "" {MPLTEXT 1 0 88 "evalf(mu(28)/(lambda(28)+mu(28))*mu(27)/(lam
bda(27)+mu(27))*mu(26)/(lambda(26)+mu(26)));" }}{PARA 11 "" 1 ""
{XPPMATH 20 "6#$\"+K\")3&R\"!#5" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 2
"Q5" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "r
estart: L:=8.7:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "g:=u/2*
exp(-u^2/4): G:=1+int(g,u):" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 2 "a
)" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 10 "t:=
47/60.:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 52 "pt:=int(1-G,u=0.
.t)/t: Lx:=L*t*pt; Ly:=L*t*(1-pt);" }}{PARA 11 "" 1 "" {XPPMATH 20 "
6#>%#LxG$\"+]'))>['!\"*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#LyG$\"+,
N6IL!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "Lx^7/7!*exp(-Lx)
*Ly*exp(-Ly);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+p^`&[$!#6" }}}
{EXCHG {PARA 0 "" 0 "" {TEXT -1 2 "b)" }{MPLTEXT 1 0 0 "" }}}{EXCHG
{PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "Linf:=evalf(L*int(u*g,u=0..infinity
));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%LinfG$\"+][.U:!\")" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "1-add(Linf^i/i!,i=0..20)*exp
(-Linf);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+tN'y,\"!#5" }}}{EXCHG
{PARA 0 "" 0 "" {TEXT -1 2 "c)" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "
> " 0 "" {MPLTEXT 1 0 17 "int(g,u=0..2.)^3;" }}{PARA 11 "" 1 ""
{XPPMATH 20 "6#$\"+yX!e_#!#5" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 2 "Q6
" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "rest
art: with(linalg):" }}{PARA 7 "" 1 "" {TEXT -1 80 "Warning, the protec
ted names norm and trace have been redefined and unprotected\n" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 101 "A:=matrix(5,5,[3/2,-49/2,-6
,-13/2,-21/2,-9/2,-21/2,-2,-13/2,-11/2,1,9,4,3,4,6,38,8,16,16,6,6,0,6,
5]):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "w:=eigenvalues(A);
" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"wG6'\"\"$\"\"%\"\"&\"\"#F)" }}
}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 87 "f:=(i,j)->w[j]^(i-1): M:
=matrix(5,5,f): for i to 5 do M[i,5]:=(i-1)*2^(i-2) end do:" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "C:=linsolve(M,[seq(a^k,k=0..
4)]):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 61 "R:=evalf(evalm(sub
s(a=A,add(log(w[i])*C[i],i=1..4)+C[5]/2)));" }}{PARA 11 "" 1 ""
{XPPMATH 20 "6#>%\"RG-%'matrixG6#7'7'$\"*Xog$R!\"*$!+-lI<')F,$!+T:Wz?F
,$!+x=)*yBF,$!+