{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 58 "assume(u<.1): M:=int(exp(u*t )*exp(1-t/2)/2,t=2..infinity);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\" MG,$*&-%$expG6#,$*&\"\"#\"\"\"%#u|irGF-F-F-,&*&F,F-F.F-F-F-!\"\"F1F1" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "P:=subs(u=5.3*(z-1),M):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 60 "m:=eval(diff(P,z),z=1); \+ sqrt(eval(diff(P,z,z),z=1.)+m-m^2);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6 #>%\"mG$\"++++?@!\")" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+<;ob6!\") " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "1-eval(mtaylor(P,z,21), z=1);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+Y%f44%!#5" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 2 "Q2" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 23 "restart: with(linalg):" }}{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 98 "A:=evalm(matrix([[0,0,0,0,0],[9,-60,23,18,10],[33,11,-64,20,0],[41,0, 27,-77,9],[0,0,0,0,0]])/10.):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "exp(-17/60*7.7);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+]G`G6!#5 " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 27 "exponential(A,17/60.)[4 ,4];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+,t2_9!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "exponential(A,100.)[4,1];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"*Y&f@&)!\"*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 61 "F:=exponential(A,t)[4,1]+exponential(A,t)[4,5]+.2*10^ (-10)*t;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"FG,**&$\"+$o*>LA!#5\" \"\"-%$expG6#,$*&$\"+Bdbcr!\"*F*%\"tGF*!\"\"F*F4*&$\"+P%4v5\"F)F*-F,6# ,$*&$\"+Gm4)R*F2F*F3F*F4F*F4*&$\"+o3HfmF)F*-F,6#,$*&$\"+\\wMXNF2F*F3F* F4F*F4$\"+++++5F2F*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "m:=i nt(diff(F,t)*t,t=0..infinity);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\" mG$\"+!o6#3B!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 43 "sqrt(int (diff(F,t)*(t-m)^2,t=0..infinity));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6 #$\"+6k.GD!#5" }}}{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 44 "lambda:=n->piecewise(n<3,1,0 .82^(n-2))*17.3:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "mu:=n-> 60/(4+26/60.)*piecewise(n<3,n,3):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "p:=n->product(lambda(j-1)/mu(j),j=1..n):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "S:=add(p(n),n=0..30);" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#>%\"SG$\"+y&=*4O!\"*" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 45 "(p(1)/S+2*p(2)/S+3*(1-(p(0)+p(1)+p(2))/S))/3;" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+Yrx9T!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 35 "add(p(n)*(1-0.82^(n-2)),n=3..30)/S;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+Yp(*HM!#6" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "add(p(n)*(n-3),n=4..30)/S;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+`*Rd>'!#6" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "add(p(n),n=6..30)/S;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+icy$ y#!#7" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 2 "Q4" }{MPLTEXT 1 0 0 "" }} }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 59 "restart: with(stats): Zt: =z->1-statevalf[cdf,normald](z):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 7 "c:=5.3:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 " a:=0.5-Zt(8/sqrt(c*3.5));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"aG$\" +Ppw$o%!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "1-2*Zt(4/sqrt (c*3.5))-a;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+0o#fy\"!#5" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 4 "1-a;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+jIB;`!#5" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 2 "Q5 " }{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 63 "A:=matrix([[6,1,-3,2],[10,8, -8,6],[-2,0,8,-2],[-10,-3,13,-4]]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6 #>%\"AG-%'matrixG6#7&7&\"\"'\"\"\"!\"$\"\"#7&\"#5\"\")!\")F*7&!\"#\"\" !F0F37&!#5F,\"#8!\"%" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "w:= eigenvalues(A);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"wG6&\"\"'\"\"%F 'F'" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 136 "solve(\{c[1]+c[2]=1 ,6*c[1]+4*c[2]+c[3]=A,\n36*c[1]+16*c[2]+8*c[3]+2*c[4]=A^2,\n6^3*c[1]+6 4*c[2]+3*16*c[3]+24*c[4]=A^3\},\{seq(c[i],i=1..4)\}):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 10 "assign(%):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 55 "B:=evalf(evalm(c[1]*ln(6)+c[2]*ln(4)+c[3]/4-c[4] /4^2));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"BG-%'matrixG6#7&7&$\"+h VH')=!\"*$\"+++++D!#5$!+++++vF/$\"+++++]F/7&$\"+C`R\"4#F,$\"+C?*>=#F,$ !+qF;w:F,$\"+;-$f=\"F,7&$!*wYge'F,$!*!QBIzF/$\"+#fJ^V#F,$!*%ypScF,7&$! +wYgLDF,$!*QB!owF,$\"+IsL6JF,$!*Bax_&F," }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "exponential(B);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-% 'matrixG6#7&7&$\"+++++g!\"*$\"+,+++5F*$!+-+++IF*$\"+++++?F*7&$\"+++++5 !\")$\"+1+++!)F*$!+2+++!)F*$\"+-+++gF*7&$!+.+++?F*$!*,8'[6!#<$\"+5+++! )F*F<7&$!+++++5F4$!+.+++IF*$\"+,+++8F4$!+.+++SF*" }}}{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, the protected names norm and trace have been redefined an d unprotected\n" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 91 "A:=matri x([[-33,33,0,0,0],[30,-71,41,0,0],[0,41,-50,9,0]\n,[0,0,12,-37,25],[0, 0,0,23,-23]]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"AG-%'matrixG6#7' 7'!#L\"#L\"\"!F,F,7'\"#I!#r\"#TF,F,7'F,F0!#]\"\"*F,7'F,F,\"#7!#P\"#D7' F,F,F,\"#B!#B" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "s:=linsolv e(transpose(A),[0,0,0,0,0]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"sG -%'vectorG6#7',$*&#\"#S\"#L\"\"\"&%#_tG6#F.F.F.,$*&#\"\"%\"\"$F.F/F.F. F2F/,$*&#\"#D\"#BF.F/F.F." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "p:=evalm(s/add(s[i],i=1..5));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#> %\"pG-%'vectorG6#7'#\"$:\"\"$m&#\"$`#\"%K6F,#\"$f(\"%GX#\"$D)F1" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "exponential(A/10.,12.)[3,4]; " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+uF)fn\"!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "3.7*p[4]*24*60;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+*R*)4$*)!\"(" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 2 "Q7" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 65 "r estart: # g(exp(-t)*exp(-exp(z)) = g(t+exp(z)) divided by z" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 69 "z*exp(z)*diff(g(t+exp(z))/z, t)-z*diff(g(t+exp(z))/z,z)-g(t+exp(z))/z:" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 12 "simplify(%);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\" \"!" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 32 "Answer: ln(t+exp(z ))^2/z" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 78 "z*exp(z)*diff(ln(t+exp(z))^2/z,t)-z*diff(ln(t+exp(z))^2/z,z)-ln(t+ exp(z))^2/z:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "simplify(%) ;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"\"!" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 2 "Q8" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 126 "restart: with(stats): \nPr:=proc(m,s,x):statevalf[c df,normald[m,s]](x)end:\nc:=5.3: t1:=7.5: t2:=9.75: x1:=-2.5: x2 :=0.4:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 60 "t:=12: Pr(x2,sqrt (c*(t-t2)),1.6)-Pr(x2,sqrt(c*(t-t2)),-0.3);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+ZI4i@!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 5 "t:=9:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 68 "m:=(x1*(t2-t)+x 2*(t-t1))/(t2-t1): s:=sqrt(c*(t2-t)*(t-t1)/(t2-t1)):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "Pr(m,s,1.6)-Pr(m,s,-0.3);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+=dSLM!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 51 "evalf(2/Pi*arccos(sqrt((10.75-9.75)/(11.5-9.75))));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+;0rVX!#5" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 2 "Q9" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 49 "lambda:=0.37: mu:=1/(1+7/24): a:=2.1: i:=7:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "e:=exp(-(lambda-mu)*t):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 59 "p:=(lambda-mu)*e/(lambda-mu*e): r:=mu*(1-e )/(lambda-mu*e):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 58 "P:=(r+( 1-r)*p*z/(1-(1-p)*z))^i*(p/(1-(1-p)*z))^(a/lambda):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "t:=4: 1-eval(mtaylor(P,z,13),z=1);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#$\"*lD/<$!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 80 "P:=((1-lambda/mu)/(1-lambda/mu*z))^(a/lambda): e val(diff(P,z),z=1); eval(P,z=0);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$ \"+G2`&>&!\"*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+BkC+D!#6" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 5 "a/mu;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+++]7F!\"*" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 3 "Q1 0" }{MPLTEXT 1 0 0 "" }}}{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 76 "for k to 6 do w[k+1]:=(1+mu/lambda)*w[k]-mu/lambda*w[ k-1]-1/k/lambda end do;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\"wG6#\" \"#$\"+&)pJHF!\"*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\"wG6#\"\"$$\" +j?P8M!\"*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\"wG6#\"\"%$\"+h#)zVR !\"*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\"wG6#\"\"&$\"+7V*zP%!\"*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\"wG6#\"\"'$\"+'4sfu%!\"*" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\"wG6#\"\"($\"+yQ[l]!\"*" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "t:=5.: 1-r^i;" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#$\"+!)o.cT!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 50 "add(1./k,k=1..7)/mu; sqrt(add(1./k^2,k=1..7))/mu;" } }{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+Vr5\\L!\"*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+r2<)e\"!\"*" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 3 " Q11" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "re start:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "g:=piecewise(u<20 ,0,exp(1-u/20)/20):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "G:=i nt(g,u):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "lambda:=6.3/60: " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "p:=int(1-G,u=0..t)/t:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 45 "t:=49.: Lx:=lambda*t*p; \+ Ly:=lambda*t*(1-p);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#LxG$\"+&R-u q$!\"*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#LyG$\"+0wfP9!\"*" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 62 "add(Lx^k/k!*exp(-Lx),k=0..2) *(1-add(Ly^k/k!*exp(-Ly),k=0..3));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6# $\"+Ms&)\\;!#6" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 "Linf:=lam bda*int(u*g,u=20..infinity);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%Lin fG$\"+++++U!\"*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "exp(-Lin f)*lambda*60*24;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+:7LnA!\"*" }}} }{MARK "24 0 1" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }