{VERSION 2 3 "DEC ALPHA UNIX" "2.3" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 }{CSTYLE "2D Math" -1 2 "Times" 0 1 0 0 0 0 0 0 2 0 0 0 0 0 0 }{CSTYLE "2D Comment" 2 18 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 } {CSTYLE "2D Output" 2 20 "" 0 1 0 0 255 1 0 0 0 0 0 0 0 0 0 }{PSTYLE " Normal" -1 0 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple Output" 0 11 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 3 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 11 12 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }1 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple Plot" 0 13 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }} {SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 568 "invGAMMA:=proc(x::a lgebraic)\n local z,f,n;\n if type(x,`!`) then op(1,x)+1\n elif typ e(x,posint) then\n z:=`evalf/invGAMMA`(x);\n f:=1;\n for n from 2 to z do\n f:=n*f;\n if f=x then RETURN(n+1) fi\n od;\n \+ 'invGAMMA(x)'\n elif type(x/Pi^(1/2),rational) and x/Pi^(1/2)>1 then \n z:=`evalf/invGAMMA`(x);\n f:=1/2;\n for n from 2 to z do\n f:=(n-1/2)*f;\n if f=x/Pi^(1/2) then RETURN(n+1/2) fi\n od; \n 'invGAMMA(x)'\n elif type(x,float) then\n `evalf/invGAMMA`(x) \n else\n 'invGAMMA'(x)\n fi\nend:\ninvGAMMA(1):=2:\ninvGAMMA(3/4 *Pi^(1/2)):=5/2:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 251 "`evalf /invGAMMA`:=proc(x::algebraic)\n local xx;\n xx:=evalf(x);\n if xx> =1 then \n fsolve(GAMMA(t)=xx,t=2..infinity)\n else\n ERROR(`ev alf/invGAMMA expects its argument to be not smal ler than 1, but received `,x)\n fi\nend:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "seq(invGAMMA(n),n=1..6);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6(\"\"#\"\"$-%)invGAMMAG6#F$-F&6#\"\"%-F&6#\"\"&F*" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "invGAMMA(N!);" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#,&%\"NG\"\"\"F%F%" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 35 "seq(invGAMMA(GAMMA(n+1/2)),n=2..5);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6&#\"\"&\"\"##\"\"(F%#\"\"*F%#\"#6F%" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "evalf(invGAMMA(10));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+^w2!R%!\"*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "invGAMMA(10.);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+^w2!R%!\"* " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "GAMMA(invGAMMA(x));" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#-%&GAMMAG6#-%)invGAMMAG6#%\"xG" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 41 "subsop(3=NULL,readlib(`simpl ify/GAMMA`));" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#:6#%\"fG6#%\"gG6\"F(C )>8$--%(readlibG6#.%-normal/GAMMAG6#9$>F+-%(collectG6%F+%&GAMMAG%'norm alG>F+-%7simplify/GAMMA/combineG6#F+>F+F5>F+-%8simplify/GAMMA/products GF=>F+-F-F=F5F(F(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 385 "`mysi mplify/GAMMA`:=proc(f)\n local g,h;\n if type(f,\{table, procedure, \+ name, numeric\}) then f\n elif not has(f,\{GAMMA,invGAMMA\}) then f \+ \n else \n g:=readlib(`simplify/GAMMA`)(f);\n if type(g,functio n) and op(0,g)=GAMMA then\n h:=op(1,g);\n if type(h,function ) and op(0,h)=invGAMMA then\n op(1,h)\n else g\n fi\n else\n map(procname,g)\n fi\n fi\nend:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 39 "`mysimplify/GAMMA`(GAMMA(invGAMMA(x)));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#%\"xG" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "invGAMMA(GAMMA(x));" }}{PARA 11 "" 1 "" {XPPMATH 20 " 6#-%)invGAMMAG6#-%&GAMMAG6#%\"xG" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 39 "`mysimplify/GAMMA`(invGAMMA(GAMMA(x)));" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#-%)invGAMMAG6#-%&GAMMAG6#%\"xG" }}}{EXCHG {PARA 0 "" 0 "" {XPPEDIT 18 0 "sqrt(x)^2 = x;" "/*$-%%sqrtG6#%\"xG\"\"#F'" } }}{EXCHG {PARA 0 "" 0 "" {XPPEDIT 18 0 "sqrt(x^2) = x;" "/-%%sqrtG6#*$ %\"xG\"\"#F'" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 421 "`mysimplif y/GAMMA/symbolic`:=proc(f)\n local g,h;\n if type(f,\{table, procedu re, name, numeric\}) then f\n elif not has(f,\{GAMMA,invGAMMA\}) then f \n else \n g:=readlib(`simplify/GAMMA`)(f);\n if type(g,func tion) then\n h:=op(1,g);\n if type(h,function) then\n \+ if \{op(0,g),op(0,h)\}=\{GAMMA,invGAMMA\} then\n op(1,h)\n \+ else g\n fi\n fi\n else\n map(procname,g)\n \+ fi\n fi\nend:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 48 "`mysim plify/GAMMA/symbolic`(GAMMA(invGAMMA(x)));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#%\"xG" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 48 "`my simplify/GAMMA/symbolic`(invGAMMA(GAMMA(x)));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#%\"xG" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 50 "`my simplify/GAMMA/symbolic`(invGAMMA(x*GAMMA(x)));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,&%\"xG\"\"\"F%F%" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 64 "`diff/invGAMMA`:=proc(a,x)\n 1/a/Psi(invGAMMA(a))*di ff(a,x)\nend:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 27 "diff(invGA MMA(t*exp(t)),t);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#**%\"tG!\"\"-%$ex pG6#F$F%-%$PsiG6#-%)invGAMMAG6#*&F$\"\"\"F&F0F%,&F&F0F/F0F0" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "asympt(ln(GAMMA(z)),z,5);" } }{PARA 11 "" 1 "" {XPPMATH 20 "6#,,*&,&-%#lnG6#%\"zG\"\"\"!\"\"F*F*F)F *F*-F'6#*&\"\"##F*F/%#PiGF0F*F&#F+F/*$F)F+#F*\"#7-%\"OG6#*$F)!\"$F*" } }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 46 "lb:=convert(asympt(ln(GAMM A(z)),z,3),polynom);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#lbG,(*&,&-% #lnG6#%\"zG\"\"\"!\"\"F,F,F+F,F,-F)6#*&\"\"##F,F1%#PiGF2F,F(#F-F1" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 46 "ub:=convert(asympt(ln(GAMMA( z)),z,5),polynom);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#ubG,**&,&-%#l nG6#%\"zG\"\"\"!\"\"F,F,F+F,F,-F)6#*&\"\"##F,F1%#PiGF2F,F(#F-F1*$F+F-# F,\"#7" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "triple:=[lb,ln(GA MMA(z)),ub];" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%'tripleG7%,(*&,&-%#l nG6#%\"zG\"\"\"!\"\"F-F-F,F-F--F*6#*&\"\"##F-F2%#PiGF3F-F)#F.F2-F*6#-% &GAMMAGF+,*F'F-F/F-F)F5*$F,F.#F-\"#7" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 61 "plot(triple,z=2..10,color=[blue,red,blue],thickness=[ 1,3,1]);" }}{PARA 13 "" 1 "" {INLPLOT "6'-%'CURVESG6%7S7$$\"\"#\"\"!$! 3,H4a&fpS8%!#>7$$\"3LLLLLsPu@!#<$\"36D*pwJ\"47XF-7$$\"3!pmm;$>5EBF1$\" 3Y#)*f&*o1UM\"!#=7$$\"3LLLLBAt'\\#F1$\"3(ffR#H)==\\#F97$$\"3aLLLjN\\oE F1$\"3V8CG#))>&)y$F97$$\"3[mmm^&Q%RGF1$\"3:xqa])G&4_F97$$\"3GLLLQk#z*H F1$\"3k&zDUWw`j'F97$$\"35+++l9.iJF1$\"3#oJY?&zw9#)F97$$\"3?LLL=\"\\#>&**F97$$\"3\")*****\\[A4]$F1$\"3:))y3\\9Dy6F17$$\"3xmmm '3Q\\n$F1$\"3!yr(4qvHw8F17$$\"3NLLLB6@GQF1$\"3(HP@b65&e:F17$$\"3&)**** **f-w+SF1$\"3,Q,0!zF>x\"F17$$\"3&*********y,uTF1$\"3Bo!zG,[Y*>F17$$\"3 U+++SP)4M%F1$\"3yI8*o#o!p@#F17$$\"3ILLL=Zg#\\%F1$\"369FA)\\`\\U#F17$$ \"36mmmEn*Gn%F1$\"33,@,JqszEF17$$\"3Tmmm1xiD[F1$\"3\">^-/&Gf,HF17$$\"3 O*****\\9!H.]F1$\"3?c'H#*fwj;$F17$$\"3%emmm]^0;&F1$\"3h5k'GdalS$F17$$ \"3<+++X@4L`F1$\"3))HHDp,6wOF17$$\"3h*****\\j\"R(\\&F1$\"32\"pz!)=q%QR F17$$\"3Kmmm;4#)ocF1$\"3)*>#p$[W#z@%F17$$\"3jmmm6lCEeF1$\"3=oo80\\_zWF 17$$\"3ELLL$G^g*fF1$\"34!)*y6$4(ow%F17$$\"3oKLL=2VshF1$\"3%>`a#ew&32&F 17$$\"3f*****\\`pfK'F1$\"3]?,'*Q7#*R`F17$$\"3\"HLLLm&z\"\\'F1$\"3PJ\") y!QJ]j&F17$$\"3s******z-6jmF1$\"3Xy9r'H*pWfF17$$\"3<******4#32$oF1$\"3 \\b)*>7EB_iF17$$\"3P*****\\#y'G*pF1$\"3MQCyl#*)Rb'F17$$\"3H******H%=H< (F1$\"3ghLw*RbP*oF17$$\"35mmm1>qMtF1$\"3?CO'*p#yJ?(F17$$\"3%)*******HS u](F1$\"3%HLN^2fx`(F17$$\"3(HLL$ep'Rm(F1$\"3>qAzlKgWyF17$$\"3()******R >4NyF1$\"3YKBLe$fR=)F17$$\"3$emm;@2h*zF1$\"3DQ4yv&**o])F17$$\"3^***** \\c9W;)F1$\"3[I$fT)o9[))F17$$\"3Lmmmmd'*G$)F1$\"36&3G&3\"R`=*F17$$\"3j *****\\iN7])F1$\"3k,Epr&R?a*F17$$\"3lKLLt>:n')F1$\"3L&eO4$*z!*))*F17$$ \"3AKLL.a#o$))F1$\"367'e))GXZ-\"!#;7$$\"3ammm^Q40!*F1$\"3k()yhIeig5Fgx 7$$\"3z******z]rf\"*F1$\"3FSFP%4*)Q4\"Fgx7$$\"3rlmmc%GpL*F1$\"3bf/i*z] B8\"Fgx7$$\"3:KLL8-V&\\*F1$\"3,y^'p%=0n6Fgx7$$\"3>+++XhUk'*F1$\"3xF$[* yqN/7Fgx7$$\"3=+++:o$)F-7$F5$\"3s5![+Dc.q\"F97$F;$\"3:\\w$4JxQ#GF97$F@$ \"3v+[bHvR*4%F97$FE$\"3/xr$*o-%=]&F97$FJ$\"3,j)egzXB\"pF97$FO$\"3&>'GV OkXx%)F97$FT$\"3$Rn=!*yI,-\"F17$FY$\"35^-N)Q\"*>?\"F17$Fhn$\"3;M[OW)=* )R\"F17$F]o$\"3@!Gt8uH-e\"F17$Fbo$\"3Ge\\EFXr#z\"F17$Fgo$\"3J`e-$>vX,# F17$F\\p$\"3$f&\\!HCqgB#F17$Fap$\"3\\v#3ZGsMW#F17$Ffp$\"3F`5k([Lvp#F17 $F[q$\"39V!\\dLP)=HF17$F`q$\"3\"G7\"\\t.,$=$F17$Feq$\"3xS$znr#oAMF17$F jq$\"3Uxr\"p$F17$F_r$\"3'o8voK7O&RF17$Fdr$\"3i\"eWli4EB%F17$Fir$\" 3<:iUmS\"Q\\%F17$F^s$\"3k7FZ\"=c2y%F17$Fcs$\"3j6h3FoM%3&F17$Fhs$\"3[B+ TaN3``F17$F]t$\"3EM3j5!eyk&F17$Fbt$\"3\")p]\"zj'>dfF17$Fgt$\"3#y9gxvBW E'F17$F\\u$\"3$HV!3$4)*ec'F17$Fau$\"3R9s^'ol`!pF17$Ffu$\"3vGUC$yKX@(F1 7$F[v$\"3>!fXDk_)[vF17$F`v$\"3flF$)=0ZbyF17$Fev$\"3^YIz>&*e%>)F17$Fjv$ \"3U?eP/fJ<&)F17$F_w$\"3^n@-)p[$e))F17$Fdw$\"3.3>Rk&R`>*F17$Fiw$\"3\"* p?\"[cP=b*F17$F^x$\"3M+FvE0p)*)*F17$Fcx$\"3l!3+7\"zoD5Fgx7$Fix$\"3)=H; P&3bh5Fgx7$F^y$\"3+he)\\^)z%4\"Fgx7$Fcy$\"3K]M#G(HCL6Fgx7$Fhy$\"383CnQ \"Hz;\"Fgx7$F]z$\"3omJ(4/>_?\"Fgx7$Fbz$\"3B())z3#>?T7Fgx7$Fgz$\"34Z\"3 ![F=!G\"Fgx-F\\[l6&F^[lF_[lF*F*-Fc[l6#\"\"$-F$6%7S7$F($\"3dUPd7rqfK!#@ 7$F/$\"3b.#)H/jgW$)F-7$F5$\"3/y3([$)fCq\"F97$F;$\"3:SmQO%)eDGF97$F@$\" 3E;')RRf!35%F97$FE$\"3.4cf+T,.bF97$FJ$\"3Yr9D^jM8pF97$FO$\"3x0^D3;Jy%) F97$FT$\"3OEz`zS??5F17$FY$\"3F/G9*pa??\"F17$Fhn$\"3J*\\(3*ot*)R\"F17$F ]o$\"3(*oLFL$y-e\"F17$Fbo$\"3PE%[W;dFz\"F17$Fgo$\"39LU5*y7Y,#F17$F\\p$ \"3eW#[Wq.hB#F17$Fap$\"3JS`4+D]VCF17$Ffp$\"3ylv8m.c(p#F17$F[q$\"3:Ns(4 wh)=HF17$F`q$\"3\"*z#efIKI=$F17$Feq$\"31h`tA*)GO&RF17$Fdr$\"3H/DBUZiKUF17$Fir$\"3!Qv9g*z#Q\\%F17$ F^s$\"3?U_jm*o2y%F17$Fcs$\"3qXFL_&eV3&F17$Fhs$\"3:-cD]W4``F17$F]t$\"3Q GK='4oyk&F17$Fbt$\"3#\\f&\\of?dfF17$Fgt$\"3BRJ)3UKWE'F17$F\\u$\"3o0xgp h!fc'F17$Fau$\"3>!R.@)F17$Fjv$\"3Ha+l88K<&)F17$F_w $\"3l(zA0y`$e))F17$Fdw$\"3B=kN_VM&>*F17$Fiw$\"3'GO&Ho?%=b*F17$F^x$\"3W l,:xZp)*)*F17$Fcx$\"38B3G7$)oD5Fgx7$Fix$\"3#4W\"yK7bh5Fgx7$F^y$\"3*=g= _())z%4\"Fgx7$Fcy$\"39/X(HJVK8\"Fgx7$Fhy$\"3x**p5i%Hz;\"Fgx7$F]z$\"3%p ()oxM>_?\"Fgx7$Fbz$\"3@T " 0 "" {MPLTEXT 1 0 39 "2*int(t/(exp(2*Pi*t)-1),t=0..infini ty);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6##\"\"\"\"#7" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "asympt(LambertW(x),x);" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#-%)LambertWG6#%\"xG" }}}{EXCHG {PARA 0 "" 0 "" {XPPEDIT 18 0 "z*exp(z) = x;" "/*&%\"zG\"\"\"-%$expG6#F$F%%\"xG" }}} {EXCHG {PARA 0 "" 0 "" {XPPEDIT 18 0 "z = ln(x)-ln(z);" "/%\"zG,&-%#ln G6#%\"xG\"\"\"-F&6#F#!\"\"" }}}{EXCHG {PARA 0 "" 0 "" {XPPEDIT 18 0 "z = ln(x)-ln(ln(x)-ln(z));" "/%\"zG,&-%#lnG6#%\"xG\"\"\"-F&6#,&-F&6#F(F )-F&6#F#!\"\"F1" }}}{EXCHG {PARA 0 "" 0 "" {XPPEDIT 18 0 "z = ln(x)-ln (ln(x))-ln(1-ln(z)/ln(x));" "/%\"zG,(-%#lnG6#%\"xG\"\"\"-F&6#-F&6#F(! \"\"-F&6#,&\"\"\"F)*&-F&6#F#F)-F&6#F(F.F.F." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 50 "subs(L[1]=ln(x),asympt(LambertW(exp(L[1])),L[1])); " }}{PARA 12 "" 1 "" {XPPMATH 20 "6#,2-%#lnG6#%\"xG\"\"\"-F%6#F$!\"\"* &F)F(F$F+F(*(F)F(,&F)F(!\"#F(F(F$F/#F(\"\"#*(F)F(,(F)!\"*\"\"'F(*$F)F1 F1F(F$!\"$#F(F5*(F)F(,**$F)\"\"$F " 0 "" {MPLTEXT 1 0 36 "equ:=L[2]-ln(lambda)-ln(1-1/lambda);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%$equG,(&%\"LG6#\"\"#\"\"\"-%#lnG6#%' lambdaG!\"\"-F,6#,&F*F**$F.F/F/F/" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 111 "da:=L[2]:\nto 6 do\n da:=asympt(map(collect,subs(ln (L[2])=L[3],asympt(subs(lambda=da,equ),L[2])),L[3]),L[2])\nod;" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%#daG,2&%\"LG6#\"\"#\"\"\"&F'6#\"\"$! \"\"*$F&F.F**$F&!\"##F*F)*$F&!\"$#F*F-*$F&!\"%#F*\"\"%*$F&!\"&#F*\"\"& -%\"OG6#*$F&!\"'F*" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%#daG,2&%\"LG6# \"\"#\"\"\"&F'6#\"\"$!\"\"*&,&F+F*F*F*F*F&F.F**&,(#F.F)F**$F+F)#F*F)F+ F*F*F&!\"#F**&,(*$F+F-#F*F-F4F*#!\"(\"\"'F*F*F&!\"$F**&,,*$F+\"\"%#F*F BF9F*F4F5F+#F>F)#!#8\"#7F*F*F&!\"%F**&,.*$F+\"\"&#F*FLFAF*F9F*F4FDF+#F %#daG,2&%\"LG6#\"\"#\"\"\"&F'6#\"\"$!\"\"*&,&F+F*F*F*F *F&F.F**&,&#F.F)F**$F+F)#F*F)F*F&!\"#F**&,**$F+F-#F*F-F4F3F+F6#F.\"\"' F*F*F&!\"$F**&,,*$F+\"\"%#F*FAF9#!\"&FF<#\"#\\\"#?F*F*F&F DF*-%\"OG6#*$F&!\"'F*" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%#daG,2&%\"L G6#\"\"#\"\"\"&F'6#\"\"$!\"\"*&,&F+F*F*F*F*F&F.F**&,&#F.F)F**$F+F)#F*F )F*F&!\"#F**&,**$F+F-#F*F-F4F3F+F.#F.\"\"'F*F*F&!\"$F**&,,*$F+\"\"%#F* FAF9#!\"&F%#daG,2&%\"LG6#\"\"#\"\"\"&F'6#\"\"$!\"\"*&,&F+ F*F*F*F*F&F.F**&,&#F.F)F**$F+F)#F*F)F*F&!\"#F**&,**$F+F-#F*F-F4F3F+F.# F.\"\"'F*F*F&!\"$F**&,,*$F+\"\"%#F*FAF9#!\"&F%#daG,2&%\"LG6#\"\"#\"\" \"&F'6#\"\"$!\"\"*&,&F+F*F*F*F*F&F.F**&,&#F.F)F**$F+F)#F*F)F*F&!\"#F** &,**$F+F-#F*F-F4F3F+F.#F.\"\"'F*F*F&!\"$F**&,,*$F+\"\"%#F*FAF9#!\"&F " 0 "" {MPLTEXT 1 0 25 "DA:=asympt(exp(da),L[2]);" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%#DAG*&,0-%$expG6#,$&%\"LG6#\"\"$!\"\"\"\"\"*(F'F0,&F+ F0F0F0F0&F,6#\"\"#F/F0*(F'F0,&*$F+F5F0F+F0F0F3!\"#F0*(F'F0,**$F+F.F0F8 #F0F5F+F/#F/F5F0F0F3!\"$F0*(F'F0,,*$F+\"\"%F0F<#F/F.F8#!\"&F5F+F/#F0\" \"'F0F0F3!\"%F0*(F'F0,,F+#\"#6FHF<#!#BFHFB#!#<\"#7*$F+\"\"&F0#FTFRF0F0 F3FFF0-%\"OG6#*$F3!\"'F0F0-F(6#F3F0" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 46 "DA:=subs(exp(L[2])=L[1],exp(-L[3])=1/L[2],DA);" }} {PARA 12 "" 1 "" {XPPMATH 20 "6#>%#DAG*&,0*$&%\"LG6#\"\"#!\"\"\"\"\"*& F(!\"#,&&F)6#\"\"$F-F-F-F-F-*&F(!\"$,&*$F1F+F-F1F-F-F-*&F(!\"%,**$F1F3 F-F7#F-F+F1F,#F,F+F-F-F-*&F(!\"&,,*$F1\"\"%F-F;#F,F3F7#F?F+F1F,#F-\"\" 'F-F-F-*&F(!\"',,F1#\"#6FFF;#!#BFFFA#!#<\"#7*$F1\"\"&F-#FRFPF-F-F--%\" OG6#*$F(FHF-F-&F)6#F-F-" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 " approx:=convert(DA,polynom);" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%'app roxG*&,.*$&%\"LG6#\"\"#!\"\"\"\"\"*&F(!\"#,&&F)6#\"\"$F-F-F-F-F-*&F(! \"$,&*$F1F+F-F1F-F-F-*&F(!\"%,**$F1F3F-F7#F-F+F1F,#F,F+F-F-F-*&F(!\"&, ,*$F1\"\"%F-F;#F,F3F7#F?F+F1F,#F-\"\"'F-F-F-*&F(!\"',,F1#\"#6FFF;#!#BF FFA#!#<\"#7*$F1\"\"&F-#FRFPF-F-F-F-&F)6#F-F-" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 113 "plot([[GAMMA(x),x,x=2..5],[x,subs(z=approx,L[3] =ln(L[2]),L[2]=ln(L[1]),L[1]=ln(x),z),x=3..30]],color=[red,blue]);" }} {PARA 13 "" 1 "" {INLPLOT "6$-%'CURVESG6$7S7$$\"\"\"\"\"!$\"\"#F*7$$\" 3Cyv_,>E!e5F0$\"34+]PC#)GA@F07$$\"3MF KnY6m$4\"F0$\"3%****\\Peui=#F07$$\"3$)eU0XxWL6F0$\"3')***\\i3&o]AF07$$ \"3R(e$Ha8=x6F0$\"3#)**\\(oX*y9BF07$$\"3@#3hmH[;A\"F0$\"3:+]P9CAuBF07$ $\"3q[cQJ.$>F\"F0$\"3\")**\\P*zhdV#F07$$\"3\"[jr9J&yG8F0$\"3,+]P>fS*\\ #F07$$\"3ChVbcYt!R\"F0$\"34+](=$f%Gc#F07$$\"3ovPKLrQg9F0$\"3$*****\\#y ,\"GEF07$$\"3/r)4Ys4r_\"F0$\"37++Dr\"zbo#F07$$\"3q%=\"p#*\\t3;F0$\"3%* ****\\(4&G]FF07$$\"39d,#o*H@)p\"F0$\"34++]7nD:GF07$$\"3rvE2OF>#z\"F0$ \"3%*****\\-*oy(GF07$$\"3;H['fj)o%)=F0$\"35+]PpnsMHF07$$\"3PQp3i&=V+#F 0$\"3,++]siL-IF07$$\"3h!o$p6!pX6#F0$\"3,+++!R5'fIF07$$\"3mP/kOI0aAF0$ \"3)***\\P/QBEJF07$$\"3W0Dym\"3&)Q#F0$\"3\"******\\\"o?&=$F07$$\"3TSW` %\\A!\\DF0$\"3%)**\\Pa&4*\\KF07$$\"3I9-1jto:FF0$\"3&)**\\7j=_6LF07$$\" 3)*fcLbyM0HF0$\"34++vVy!eP$F07$$\"3XyGKx-/&4$F0$\"3K+](=WU[V$F07$$\"3! *pR/Wy#zJ$F0$\"3(****\\7B>&)\\$F07$$\"3l?571(>:d$F0$\"3)***\\P>:mkNF07 $$\"3y6It_bR7QF0$\"3'***\\iv&QAi$F07$$\"3`tUg1um&4%F0$\"31++vtLU%o$F07 $$\"3Fu$o#y'>gT%F0$\"37+++bjm[PF07$$\"3?rp_A+YfZF0$\"3\"****\\(yb^6QF0 7$$\"3XYo#)y4$H7&F0$\"3****\\PMaKsQF07$$\"3:X3JsY?mbF0$\"3=++D6W%)RRF0 7$$\"3.C&)*>UpQ+'F0$\"3z*****\\@80+%F07$$\"3\"ys9g]In^'F0$\"3h****\\7, HlSF07$$\"3c%oa&))=HEqF0$\"3')**\\P4w)R7%F07$$\"3z*HB'zgKPwF0$\"3<++]x %f\")=%F07$$\"3l%*4iO:3p#)F0$\"3!)**\\P/-a[UF07$$\"3EkZwO#pY**)F0$\"3[ +](=Yb;J%F07$$\"3\"=Es@K&Rv(*F0$\"3V****\\i@OtVF07$$\"3$fm<[\\iw1\"!#; $\"3')**\\PfL'zV%F07$$\"3.(Q'*e)oYj6Fiw$\"3u*******)>=+XF07$$\"3O$>ZJ5 .;F\"Fiw$\"3.++DE&4Qc%F07$$\"3L3))RUy7!R\"Fiw$\"3<+]P%>5pi%F07$$\"3SJ \\)Q.M+^\"Fiw$\"39+++bJ*[o%F07$$\"3+@JMZZ'=m\"Fiw$\"34++Dr\"[8v%F07$$ \"3S**RH]197=Fiw$\"3++++Ijy5[F07$$\"3\"o)39Z+:*)>Fiw$\"32+]P/)fT([F07$ $\"3%H1q!=igw@Fiw$\"3^+]i0j\"[$\\F07$$\"3O+++++++CFiw$\"\"&F*-%'COLOUR G6&%$RGBG$\"*++++\"!\")F*F*-F$6$7in7$$\"\"$F*$!3$e))\\c(4o(4\"!\"*7$$ \"3)[7G)[8R=IF0$!38c(G<,'>-o!#57$$\"3v\\il(p#yOIF0$!3Z)*QJ&zS.L%F^\\l7 $$\"3juV[YS))!#67$$\"3!***\\i!zIr9$F0$!30))y=\"eC VS%Fc]l7$$\"3%)*\\Pf='p?KF0$!3$f]qA=\"*QF\"Fc]l7$$\"3z***\\7ehUH$F0$!3 zk79:\\P$G%!#77$$\"3u*\\il(p#yO$F0$!3&z#*G)4H*)4;Fc^l7$$\"3p**\\(=P#RT MF0$!3\"y'p_6!*Q\"e'!#87$$\"33+v=nx&\\^$F0$!3M)>!>%QB\"oGF^_l7$$\"3.++ ]iJ_)e$F0$!3mgHw\\[]68F^_l7$$\"3'**\\P4feX%QF0$!3!)f#f(>`)G/\"!#97$$\" 3M+]P>Sf+TF0$!35+E'=F8_7'Fiw7$$\"3&****\\PDrkn%F0$\"3cw(**3&)Q&))fF07$ $\"3[++Dwd;c_F0$!3CrsP\"\\;Wp\"F07$$\"3p**\\(=60J$eF0$!3@K*3.\"p'y@\"! #=7$$\"3)***\\PH<+ojF0$\"30_K#3=-F\">F07$$\"3m**\\P%>c=#pF0$\"3.m/\\v3 ::MF07$$\"3^+]PuKl%\\(F0$\"3NHSeJW())Q%F07$$\"3S+](oQ8c1)F0$\"3yz:*\\o gC(\\F07$$\"3N****\\Ug\"Hl)F0$\"3'f46%*>/nK&F07$$\"3&))**\\7a7-<*F0$\" 3'R%f0&*GU7bF07$$\"3P++]xec_(*F0$\"3UR!\\lFz[j&F07$$\"3!****\\7/JP.\"F iw$\"3\"*=.F)3#z+dF07$$\"33++D7?3!4\"Fiw$\"3S#R*pT^GJdF07$$\"3'**\\PC4 a79\"Fiw$\"3_;iB)HE:u&F07$$\"3#****\\_k-@?\"Fiw$\"3yUI.InVSdF07$$\"3(* *****4N\\OD\"Fiw$\"3>)3NR;]Et&F07$$\"3%**\\PRU5OJ\"Fiw$\"3bn\"R<\"Q6>d F07$$\"3,++]Lhom8Fiw$\"3KO'4yJ**\\q&F07$$\"3/+v$*)f=\\U\"Fiw$\"3f0CMn] d)o&F07$$\"30+D\"onp.[\"Fiw$\"3ly+>)=TGn&F07$$\"3/+]PfqAQ:Fiw$\"3*R8MN =vol&F07$$\"36+vo(>e8f\"Fiw$\"3)))ypyp3Hk&F07$$\"3-+]73tm[;Fiw$\"3i4R( Rcq(GcF07$$\"32+vVn`>3(*zt=Fiw$\"38; T\"yn5Oe&F07$$\"3,+](3-k.$>Fiw$\"3E[JmFiw$ \"3i6dP\">wrc&F07$$\"3!)**\\7q*fe/#Fiw$\"3&e8yY/m'fbF07$$\"3j****\\$*= Y+@Fiw$\"37fcghSs`bF07$$\"3=++D,6we@Fiw$\"3.mafY@;[bF07$$\"3))*\\P%[)) e6AFiw$\"3PL76e1xVbF07$$\"3C++vHNMpAFiw$\"3d4nr[FiRbF07$$\"3\"**\\PR=' oBBFiw$\"3k5d^UoHObF07$$\"3!**\\(o:**[!Q#Fiw$\"3>^A#eUsL`&F07$$\"3-++D Yf-OCFiw$\"3vW3#y9<5`&F07$$\"3'**\\PM-nT\\#Fiw$\"3#p-sr=V!HbF07$$\"3I+ ++\"zj,b#Fiw$\"3B//)**=!eFbF07$$\"36+]it&Gug#Fiw$\"3+$**4lF$\\EbF07$$ \"3;+v$\\<>Um#Fiw$\"3-F))[q;zDbF07$$\"3A++]RQS;FFiw$\"3**y\"H>s^a_&F07 $$\"3<+]7aL@wFFiw$\"3@*QvJ$GRDbF07$$\"3\"******pp2(HGFiw$\"32#=:^M;c_& F07$$\"3:+v$R#Qu')GFiw$\"3fULlQ\">h_&F07$$\"3$**\\i]nM8%HFiw$\"3p@G6#Q Oo_&F07$$\"#IF*$\"3%fZ_a7Wy_&F0-Fjz6&F\\[lF*F*F][l" 2 258 258 258 2 0 1 0 2 9 0 4 2 1.000000 45.000000 45.000000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 30 0 0 0 0 0 0 }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 113 "plot([[GAMMA(x),x,x=2..5],[x,subs(z=approx,L[3]=ln(L [2]),L[2]=ln(L[1]),L[1]=ln(x),z),x=5..30]],color=[red,blue]);" }} {PARA 13 "" 1 "" {INLPLOT "6$-%'CURVESG6$7S7$$\"\"\"\"\"!$\"\"#F*7$$\" 3Cyv_,>E!e5F0$\"34+]PC#)GA@F07$$\"3MF KnY6m$4\"F0$\"3%****\\Peui=#F07$$\"3$)eU0XxWL6F0$\"3')***\\i3&o]AF07$$ \"3R(e$Ha8=x6F0$\"3#)**\\(oX*y9BF07$$\"3@#3hmH[;A\"F0$\"3:+]P9CAuBF07$ $\"3q[cQJ.$>F\"F0$\"3\")**\\P*zhdV#F07$$\"3\"[jr9J&yG8F0$\"3,+]P>fS*\\ #F07$$\"3ChVbcYt!R\"F0$\"34+](=$f%Gc#F07$$\"3ovPKLrQg9F0$\"3$*****\\#y ,\"GEF07$$\"3/r)4Ys4r_\"F0$\"37++Dr\"zbo#F07$$\"3q%=\"p#*\\t3;F0$\"3%* ****\\(4&G]FF07$$\"39d,#o*H@)p\"F0$\"34++]7nD:GF07$$\"3rvE2OF>#z\"F0$ \"3%*****\\-*oy(GF07$$\"3;H['fj)o%)=F0$\"35+]PpnsMHF07$$\"3PQp3i&=V+#F 0$\"3,++]siL-IF07$$\"3h!o$p6!pX6#F0$\"3,+++!R5'fIF07$$\"3mP/kOI0aAF0$ \"3)***\\P/QBEJF07$$\"3W0Dym\"3&)Q#F0$\"3\"******\\\"o?&=$F07$$\"3TSW` %\\A!\\DF0$\"3%)**\\Pa&4*\\KF07$$\"3I9-1jto:FF0$\"3&)**\\7j=_6LF07$$\" 3)*fcLbyM0HF0$\"34++vVy!eP$F07$$\"3XyGKx-/&4$F0$\"3K+](=WU[V$F07$$\"3! *pR/Wy#zJ$F0$\"3(****\\7B>&)\\$F07$$\"3l?571(>:d$F0$\"3)***\\P>:mkNF07 $$\"3y6It_bR7QF0$\"3'***\\iv&QAi$F07$$\"3`tUg1um&4%F0$\"31++vtLU%o$F07 $$\"3Fu$o#y'>gT%F0$\"37+++bjm[PF07$$\"3?rp_A+YfZF0$\"3\"****\\(yb^6QF0 7$$\"3XYo#)y4$H7&F0$\"3****\\PMaKsQF07$$\"3:X3JsY?mbF0$\"3=++D6W%)RRF0 7$$\"3.C&)*>UpQ+'F0$\"3z*****\\@80+%F07$$\"3\"ys9g]In^'F0$\"3h****\\7, HlSF07$$\"3c%oa&))=HEqF0$\"3')**\\P4w)R7%F07$$\"3z*HB'zgKPwF0$\"3<++]x %f\")=%F07$$\"3l%*4iO:3p#)F0$\"3!)**\\P/-a[UF07$$\"3EkZwO#pY**)F0$\"3[ +](=Yb;J%F07$$\"3\"=Es@K&Rv(*F0$\"3V****\\i@OtVF07$$\"3$fm<[\\iw1\"!#; $\"3')**\\PfL'zV%F07$$\"3.(Q'*e)oYj6Fiw$\"3u*******)>=+XF07$$\"3O$>ZJ5 .;F\"Fiw$\"3.++DE&4Qc%F07$$\"3L3))RUy7!R\"Fiw$\"3<+]P%>5pi%F07$$\"3SJ \\)Q.M+^\"Fiw$\"39+++bJ*[o%F07$$\"3+@JMZZ'=m\"Fiw$\"34++Dr\"[8v%F07$$ \"3S**RH]197=Fiw$\"3++++Ijy5[F07$$\"3\"o)39Z+:*)>Fiw$\"32+]P/)fT([F07$ $\"3%H1q!=igw@Fiw$\"3^+]i0j\"[$\\F07$$\"3O+++++++CFiw$\"\"&F*-%'COLOUR G6&%$RGBG$\"*++++\"!\")F*F*-F$6$7]o7$Fgz$!3]'3S!HJ%Q2&!#=7$$\"3!Q$3x1h 6o]F0$!3.Ea[Ms(H1\"F07$$\"3qm;a8ABO^F0$!3**zuxB5#QU\"F07$$\"3h*\\7.K[V ?&F0$!3\"3Lq@bgki\"F07$$\"3SLL3FWYs_F0$!3dq^lAzl-y&F0$!3#)oa]-wB)H$Ff[l7$$\"3p;a)ehL0! fF0$\"3pB\"=0Ua&Q:Ff[l7$$\"3iL$ek`o!>gF0$\"3hnTp6P_DjFf[l7$$\"3i[1^cy?F07$$\"3±z>)G_lF0$\"3T:\"G1B')\\[#F07$$\"3Jn;zW nZ'o'F0$\"3%F07$$\"3&GLeRZXKi(F0$\"3ZlN+]xkZXF07$$\"3q** \\(oz#)3(yF0$\"3AX'Hg6&>1[F07$$\"3am;z>,_=\")F0$\"3=(fK9KnF,&F07$$\"33 +]7G$[8j)F0$\"3+/B9>Q\"oJ&F07$$\"3(om\"z%*frh\"*F0$\"3,2tJFp05bF07$$\" 3X**\\ilFQ!p*F0$\"3$)[u:u\"4]i&F07$$\"3^LL3_\"=M-\"Fiw$\"3.#yFPy)=#p&F 07$$\"3vm;/wfJr5Fiw$\"3;QXC?\"HQs&F07$$\"3A++D\"eP_7\"Fiw$\"3?Yf1\\TnR dF07$$\"3-++v$f!Qz6Fiw$\"3P\"*\\5hp6UdF07$$\"36++v=ubJ7Fiw$\"3S`tM$fal t&F07$$\"3lm\"zW(*Q*y7Fiw$\"3sav$RKxts&F07$$\"3\\LL3F-GN8Fiw$\"3)eum3_ 7Nr&F07$$\"3LLLLe'3IQ\"Fiw$\"3^\\\"\\1C]/q&F07$$\"34+DJq\"G&Q9Fiw$\"3D 0T*zm(p%o&F07$$\"3TLL$eMsw[\"Fiw$\"3k:q3o?zqcF07$$\"3;+DJ&H\"fT:Fiw$\" 3M(4%\\V)ofl&F07$$\"31+v$f)[$Hf\"Fiw$\"3,(>#[&G1Dk&F07$$\"3^L$ek`1lk\" Fiw$\"3LIr+d[GHcF07$$\"3ULe*[.-dp\"Fiw$\"3;T,hZQ&zh&F07$$\"3ym;/Egw[Lgef&F07$$\"38+voza' =&=Fiw$\"3_hXF2uG(e&F07$$\"3tm;zWho.>Fiw$\"3C-?&R1R)ybF07$$\"3D++]i>Ad >Fiw$\"3A-)[kqb4d&F07$$\"32+]i:jf4?Fiw$\"3c-VeA)HSc&F07$$\"3B+DJ&>r-1# Fiw$\"37:*R8%f-ebF07$$\"3\"***\\P4q`;@Fiw$\"3baONTA6_bF07$$\"3MLL$eM%4 n@Fiw$\"3g)e9$)oHua&F07$$\"3@++v$4v5A#Fiw$\"3*H/xx(Q/VbF07$$\"3bm\"zWn *)*pAFiw$\"3itI&p+!eRbF07$$\"3H++DJiYBBFiw$\"3V\"QC\\B4j`&F07$$\"3>Lek .NytBFiw$\"36KCtX(*oLbF07$$\"3B+Dc^&zjU#Fiw$\"3KRNjJARJbF07$$\"3FLL3-= !yZ#Fiw$\"3UW(HuG]&HbF07$$\"3:+D\"G8O;`#Fiw$\"3Lr1(=$)=!GbF07$$\"3*pmm ;*\\[$e#Fiw$\"3LV2^a&**o_&F07$$\"3sm;aQz]OEFiw$\"3[\"))Hy9*3EbF07$$\"3 rLekG=4*o#Fiw$\"3T]!ojm%fDbF07$$\"3F++]i4TPFFiw$\"3?n7$4@#RDbF07$$\"3V L$3F9!z#z#Fiw$\"3tH4IN_VDbF07$$\"3(omm;%>KUGFiw$\"3Z^QZ&y/d_&F07$$\"3L +DJqJ8&*GFiw$\"3e\"oh[#[@EbF07$$\"3P+voa-oXHFiw$\"3([m12&G!p_&F07$$\"# IF*$\"3%fZ_a7Wy_&F0-Fjz6&F\\[lF*F*F][l" 2 258 258 258 2 0 1 0 2 9 0 4 2 1.000000 45.000000 45.000000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 117 "pl ot([[GAMMA(x),x,x=5..6.5],[x,subs(z=approx,L[3]=ln(L[2]),L[2]=ln(L[1]) ,L[1]=ln(x),z),x=20..300]],color=[red,blue]);" }}{PARA 13 "" 1 "" {INLPLOT "6$-%'CURVESG6$7S7$$\"3O+++++++C!#;$\"\"&\"\"!7$$\"3\"=QAY^S9 _#F*$\"3e***\\7t&pK]!#<7$$\"3\\!H;&pqfKEF*$\"3#)*\\(=7T9h]F37$$\"33z: \\$[_Sw#F*$\"3>+](=HPJ4&F37$$\"3fn\"4FOgO!HF*$\"3$***\\7VDMD^F37$$\"39 I:x\"z\"G]IF*$\"3p*\\P%GZRd^F37$$\"30WY&*e;X$>$F*$\"3&)*\\(=276(=&F37$ $\"3w5&)HCMW\\LF*$\"3p*\\(o**3)y@&F37$$\"3spgot*H&>NF*$\"3B+vofHq\\_F3 7$$\"3g8*pQ(zW)p$F*$\"3#)*\\Pf'HU\"G&F37$$\"3kgVf5G$G*QF*$\"3'****\\7* 309`F37$$\"3Y%z'>'pNK2%F*$\"3%)**\\i&e*yU`F37$$\"3J@ZL(*3A(G%F*$\"3U++ v[D9v`F37$$\"3Rm/(*G#oV^%F*$\"30++Dc$GwS&F37$$\"3?Q\"p2-4cu%F*$\"3v*** \\7XM*QaF37$$\"3883fl]pm\\F*$\"3F+vo%QjtY&F37$$\"35-sp`&QTC&F*$\"3c*** \\i8o6]&F37$$\"3XUr')4+C#\\&F*$\"3y*****\\>0)HbF37$$\"3E[2CL;$oz&F*$\" 3x*\\(=-p6jbF37$$\"3Qw4cT+i\"3'F*$\"3S++]2Mg#f&F37$$\"3w8c;4YX6kF*$\"3 #**\\(=xZ&\\i&F37$$\"3(p))3%[UUVnF*$\"3$**\\i:$4wbcF37$$\"3?k;i=u_4rF* $\"3#)**\\(=#R!zo&F37$$\"3g%3)Qc6^kuF*$\"3;+v$4A@ur&F37$$\"3R@J'=7t(oy F*$\"3@+]i:'f#\\dF37$$\"3Od^DM,e8$)F*$\"3V+vof2L#y&F37$$\"3d;DopQnA()F *$\"3w*\\7yG>6\"eF37$$\"3>>[$*)y_()=*F*$\"3D+](oo6A%eF37$$\"3gnV?2ZE)p *F*$\"3h****\\xJLueF37$$\"3-jZPkXfA5!#:$\"3S+]P*ydd!fF37$$\"3,HH)fNol2 \"Fht$\"3***\\(=\\?a)o7Fht$\"3!)***\\i0XE.'F37 $$\"3s*p%fW3MM8Fht$\"3$**\\(o/Q*>1'F37$$\"3V3%HYld+T\"Fht$\"3`++vQ(zS4 'F37$$\"3KTBa:#pa[\"Fht$\"3M+v=-,FChF37$$\"3hhC@1<))o:Fht$\"3C+v$4tFe: 'F37$$\"3sZQk1+Cb;Fht$\"3r***\\73\"o'='F37$$\"3e*fF#G*f5v\"Fht$\"3P+vo z;)*=iF37$$\"3$3G[#)4/*[=Fht$\"3()*****\\*44]iF37$$\"3(eK`\"4+'\\&>Fht $\"3Y+]7jZ!>G'F37$$\"3+%zz8J,l1#Fht$\"34+v=(4bMJ'F37$$\"3+.RVmf#\\<#Fh t$\"3_++]xlWUjF37$$\"3vk]sbQh1BFht$\"3\\+]i&3ucP'F37$$\"3#[oROg([k_&F37$$\"3GL$3xUzIu#F*$ \"3X4F>rMQDbF37$$\"3))***\\(Q&Qe(GF*$\"30k0T,H+EbF37$$\"3[m;z\\wf3IF*$ \"3EMO#GD6!GbF37$$\"33LL$3wc89$F*$\"3Le:g_];JbF37$$\"3o***\\7xf*RMF*$ \"3H5tY!\\h9a&F37$$\"3)pmm;yi&QPF*$\"31>`.,&))[b&F37$$\"3@nmmrusRVF*$ \"3(R[`]OGse&F37$$\"3%HLL3$\\.Q\\F*$\"3d3!f\"\\#QDi&F37$$\"3Dmm;MDu#\\ &F*$\"3^kd<#Q#*el&F37$$\"3\")****\\F,6ngF*$\"3\"Hoa4c0+p&F37$$\"3hmm;9 >7hmF*$\"3*oWv4O\"HCdF37$$\"3_++](pGKD(F*$\"3OZ'=[@1sv&F37$$\"3hLLL.LG iyF*$\"3Q`eILBl*y&F37$$\"3ZmmmJ*Q()R)F*$\"30YyoV,/VS+\"\\(eF37$$\"3&**** ***3VM>5Fht$\"3v(QH=`M5!fF37$$\"3cmmT^;Ts5Fht$\"3s?U,3G\"Q#fF37$$\"3BL LLaQ^N6Fht$\"3IvSTF8!)\\fF37$$\"3OLLL(pp*)=\"Fht$\"3kIWz`J%4(fF37$$\"3 !****\\2b^6D\"Fht$\"318Fht$\"3u#z)*3dRZ,'F37$ $\"3=++v]Aem8Fht$\"3\">0%HaH/OgF37$$\"3,++Dsq3C9Fht$\"3oNFB)*4gbgF37$$ \"3[LL$3K(3%[\"Fht$\"3!ySi&R!*HvgF37$$\"3\\LL3zi=R:Fht$\"3=)p%)3<*y#4' F37$$\"3omm;\\zh)f\"Fht$\"36(4zrc]56'F37$$\"3`mmT^2Ng;Fht$\"33SN(pk%RH hF37$$\"3(****\\sL*39'F37$$\"3!)***\\(QP]Z>Fht$\"3xa6%4&Q>xiF37$$\"37LL3Cvj)H#Fht$\"3\"p_gsJ`&*G'F37$$\"3-++v(4Xv N#Fht$\"3ONo>X9<-jF37$$\"3/LLL=!Q^T#Fht$\"3NGUA`NA9jF37$$\"3m***\\(oCV vCFht$\"3#)\\z`zHbEjF37$$\"3&omm1>.N`#Fht$\"3LoP&R:h\"QjF37$$\"3nmm;\" *)))Gf#Fht$\"3qj!4jCt(\\jF37$$\"3JLL3[Gy^EFht$\"31TA_68/hjF37$$\"3I+++ y-!fq#Fht$\"3%o;#RNl=rjF37$$\"36LL$)f\\#zw#Fht$\"3dryVUxd#Q'F37$$\"3ym mmu0SBGFht$\"3k%>WRYgDR'F37$$\"3&)***\\2:\\D)GFht$\"3gl^\\E%)*HS'F37$$ \"3&****\\_)=;RHFht$\"3hGe&fO(z7kF37$$\"$+$F-$\"3C@\\,Am7BkF3-F[[l6&F] [lF-F-F^[l" 2 258 258 258 2 0 1 0 2 9 0 4 2 1.000000 45.000000 45.000000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 -16416 -24486 0 0 0 0 0 0 }}}{EXCHG {PARA 0 "" 0 "" {XPPEDIT 18 0 "f(z) = z-(GAMMA(z)- x)/(Psi(z)*GAMMA(z));" "/-%\"fG6#%\"zG,&F&\"\"\"*&,&-%&GAMMAG6#F&F(%\" xG!\"\"F(*&-%$PsiG6#F&F(-F,6#F&F(F/F/" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 "f:= z-(GAMMA(z)-x)/(Psi(z)*GAMMA(z));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"fG,&%\"zG\"\"\"*(,&-%&GAMMAG6#F&F'%\"xG!\"\"F' -%$PsiGF,F.F*F.F." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 114 "plot( [[GAMMA(x),x,x=2..5],[x,subs(z=4.5,f),x=1..20]], \+ color=[red,blue]);" }}{PARA 13 "" 1 "" {INLPLOT " 6$-%'CURVESG6$7S7$$\"\"\"\"\"!$\"\"#F*7$$\"3Cyv_,>E!e5F0$\"34+]PC#)GA@F07$$\"3MFKnY6m$4\"F0$\"3%****\\Peui= #F07$$\"3$)eU0XxWL6F0$\"3')***\\i3&o]AF07$$\"3R(e$Ha8=x6F0$\"3#)**\\(o X*y9BF07$$\"3@#3hmH[;A\"F0$\"3:+]P9CAuBF07$$\"3q[cQJ.$>F\"F0$\"3\")** \\P*zhdV#F07$$\"3\"[jr9J&yG8F0$\"3,+]P>fS*\\#F07$$\"3ChVbcYt!R\"F0$\"3 4+](=$f%Gc#F07$$\"3ovPKLrQg9F0$\"3$*****\\#y,\"GEF07$$\"3/r)4Ys4r_\"F0 $\"37++Dr\"zbo#F07$$\"3q%=\"p#*\\t3;F0$\"3%*****\\(4&G]FF07$$\"39d,#o* H@)p\"F0$\"34++]7nD:GF07$$\"3rvE2OF>#z\"F0$\"3%*****\\-*oy(GF07$$\"3;H ['fj)o%)=F0$\"35+]PpnsMHF07$$\"3PQp3i&=V+#F0$\"3,++]siL-IF07$$\"3h!o$p 6!pX6#F0$\"3,+++!R5'fIF07$$\"3mP/kOI0aAF0$\"3)***\\P/QBEJF07$$\"3W0Dym \"3&)Q#F0$\"3\"******\\\"o?&=$F07$$\"3TSW`%\\A!\\DF0$\"3%)**\\Pa&4*\\K F07$$\"3I9-1jto:FF0$\"3&)**\\7j=_6LF07$$\"3)*fcLbyM0HF0$\"34++vVy!eP$F 07$$\"3XyGKx-/&4$F0$\"3K+](=WU[V$F07$$\"3!*pR/Wy#zJ$F0$\"3(****\\7B>&) \\$F07$$\"3l?571(>:d$F0$\"3)***\\P>:mkNF07$$\"3y6It_bR7QF0$\"3'***\\iv &QAi$F07$$\"3`tUg1um&4%F0$\"31++vtLU%o$F07$$\"3Fu$o#y'>gT%F0$\"37+++bj m[PF07$$\"3?rp_A+YfZF0$\"3\"****\\(yb^6QF07$$\"3XYo#)y4$H7&F0$\"3**** \\PMaKsQF07$$\"3:X3JsY?mbF0$\"3=++D6W%)RRF07$$\"3.C&)*>UpQ+'F0$\"3z*** **\\@80+%F07$$\"3\"ys9g]In^'F0$\"3h****\\7,HlSF07$$\"3c%oa&))=HEqF0$\" 3')**\\P4w)R7%F07$$\"3z*HB'zgKPwF0$\"3<++]x%f\")=%F07$$\"3l%*4iO:3p#)F 0$\"3!)**\\P/-a[UF07$$\"3EkZwO#pY**)F0$\"3[+](=Yb;J%F07$$\"3\"=Es@K&Rv (*F0$\"3V****\\i@OtVF07$$\"3$fm<[\\iw1\"!#;$\"3')**\\PfL'zV%F07$$\"3.( Q'*e)oYj6Fiw$\"3u*******)>=+XF07$$\"3O$>ZJ5.;F\"Fiw$\"3.++DE&4Qc%F07$$ \"3L3))RUy7!R\"Fiw$\"3<+]P%>5pi%F07$$\"3SJ\\)Q.M+^\"Fiw$\"39+++bJ*[o%F 07$$\"3+@JMZZ'=m\"Fiw$\"34++Dr\"[8v%F07$$\"3S**RH]197=Fiw$\"3++++Ijy5[ F07$$\"3\"o)39Z+:*)>Fiw$\"32+]P/)fT([F07$$\"3%H1q!=igw@Fiw$\"3^+]i0j\" [$\\F07$$\"3O+++++++CFiw$\"\"&F*-%'COLOURG6&%$RGBG$\"*++++\"!\")F*F*-F $6$7S7$F($\"+Y6*=%Q!\"*7$$\"+Hf999Ff[l$\"+op_nQFf[l7$$\"+)3#\\upFf[l$\"+2zL3UFf[l7$$\"+D(H\"[tFf[ l$\"+\"GV[B%Ff[l7$$\"+.e'3r(Ff[l$\"+6oHdUFf[l7$$\"+%49G8)Ff[l$\"+.cT$G %Ff[l7$$\"+G)4j])Ff[l$\"+g^`1VFf[l7$$\"+WQ4;*)Ff[l$\"+y4!>L%Ff[l7$$\"+ K^I1$*Ff[l$\"+N_0cVFf[l7$$\"+w'\\Mr*Ff[l$\"+pwD\"Q%Ff[l7$$\"+YNt35F_[l $\"+q8S/WFf[l7$$\"+!=i!\\5F_[l$\"+z\\OHWFf[l7$$\"+'H_44\"F_[l$\"+`^HbW Ff[l7$$\"+kxTF6F_[l$\"+Wu'yZ%Ff[l7$$\"+q9!o;\"F_[l$\"+RhC-XFf[l7$$\"+ \"p)[27F_[l$\"+:X9F_[l$\"+X7duYFf[l7$$\"+OV $e[\"F_[l$\"+V*G(*p%Ff[l7$$\"+ja2C:F_[l$\"+w.SBZFf[l7$$\"+f%[Sc\"F_[l$ \"+VQ9[ZFf[l7$$\"+p$HJg\"F_[l$\"+0^LsZFf[l7$$\"+hM/W;F_[l$\"+:6m(z%Ff[ l7$$\"+%f[Mo\"F_[l$\"+tI0A[Ff[l7$$\"+LguB0#\\Ff[l7$$\"+wY;!)=F_[l$\"+G7#Q%\\Ff[l7$$\"+47I?>F_[l$\"+Efmo\\Ff[l7 $$\"+$*pre>F_[l$\"+saW#*\\Ff[l7$$\"#?F*$\"+U)**z,&Ff[l-Fjz6&F\\[lF*F*F ][l" 2 258 258 258 2 0 1 0 2 9 0 4 2 1.000000 45.000000 45.000000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 -16423 -8411 0 0 0 0 0 0 }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "ff:=subs(z=f,f);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#ffG,(%\"zG\"\"\"*(,&-%&GAMMAG6#F&F'%\"xG! \"\"F'-%$PsiGF,F.F*F.F.*(,&-F+6#,&F&F'F(F.F'F-F.F'-F0F4F.F3F.F." }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 115 "plot([[GAMMA(x),x,x=2..5],[ x,subs(z=4.5,ff),x=1..20]], \+ color=[red,blue]);" }}{PARA 13 "" 1 "" {INLPLOT "6$-%'CURVESG6$7S7$$\" \"\"\"\"!$\"\"#F*7$$\"3Cyv_,>E!e5F0$ \"34+]PC#)GA@F07$$\"3MFKnY6m$4\"F0$\"3%****\\Peui=#F07$$\"3$)eU0XxWL6F 0$\"3')***\\i3&o]AF07$$\"3R(e$Ha8=x6F0$\"3#)**\\(oX*y9BF07$$\"3@#3hmH[ ;A\"F0$\"3:+]P9CAuBF07$$\"3q[cQJ.$>F\"F0$\"3\")**\\P*zhdV#F07$$\"3\"[j r9J&yG8F0$\"3,+]P>fS*\\#F07$$\"3ChVbcYt!R\"F0$\"34+](=$f%Gc#F07$$\"3ov PKLrQg9F0$\"3$*****\\#y,\"GEF07$$\"3/r)4Ys4r_\"F0$\"37++Dr\"zbo#F07$$ \"3q%=\"p#*\\t3;F0$\"3%*****\\(4&G]FF07$$\"39d,#o*H@)p\"F0$\"34++]7nD: GF07$$\"3rvE2OF>#z\"F0$\"3%*****\\-*oy(GF07$$\"3;H['fj)o%)=F0$\"35+]Pp nsMHF07$$\"3PQp3i&=V+#F0$\"3,++]siL-IF07$$\"3h!o$p6!pX6#F0$\"3,+++!R5' fIF07$$\"3mP/kOI0aAF0$\"3)***\\P/QBEJF07$$\"3W0Dym\"3&)Q#F0$\"3\"***** *\\\"o?&=$F07$$\"3TSW`%\\A!\\DF0$\"3%)**\\Pa&4*\\KF07$$\"3I9-1jto:FF0$ \"3&)**\\7j=_6LF07$$\"3)*fcLbyM0HF0$\"34++vVy!eP$F07$$\"3XyGKx-/&4$F0$ \"3K+](=WU[V$F07$$\"3!*pR/Wy#zJ$F0$\"3(****\\7B>&)\\$F07$$\"3l?571(>:d $F0$\"3)***\\P>:mkNF07$$\"3y6It_bR7QF0$\"3'***\\iv&QAi$F07$$\"3`tUg1um &4%F0$\"31++vtLU%o$F07$$\"3Fu$o#y'>gT%F0$\"37+++bjm[PF07$$\"3?rp_A+YfZ F0$\"3\"****\\(yb^6QF07$$\"3XYo#)y4$H7&F0$\"3****\\PMaKsQF07$$\"3:X3Js Y?mbF0$\"3=++D6W%)RRF07$$\"3.C&)*>UpQ+'F0$\"3z*****\\@80+%F07$$\"3\"ys 9g]In^'F0$\"3h****\\7,HlSF07$$\"3c%oa&))=HEqF0$\"3')**\\P4w)R7%F07$$\" 3z*HB'zgKPwF0$\"3<++]x%f\")=%F07$$\"3l%*4iO:3p#)F0$\"3!)**\\P/-a[UF07$ $\"3EkZwO#pY**)F0$\"3[+](=Yb;J%F07$$\"3\"=Es@K&Rv(*F0$\"3V****\\i@OtVF 07$$\"3$fm<[\\iw1\"!#;$\"3')**\\PfL'zV%F07$$\"3.(Q'*e)oYj6Fiw$\"3u**** ***)>=+XF07$$\"3O$>ZJ5.;F\"Fiw$\"3.++DE&4Qc%F07$$\"3L3))RUy7!R\"Fiw$\" 3<+]P%>5pi%F07$$\"3SJ\\)Q.M+^\"Fiw$\"39+++bJ*[o%F07$$\"3+@JMZZ'=m\"Fiw $\"34++Dr\"[8v%F07$$\"3S**RH]197=Fiw$\"3++++Ijy5[F07$$\"3\"o)39Z+:*)>F iw$\"32+]P/)fT([F07$$\"3%H1q!=igw@Fiw$\"3^+]i0j\"[$\\F07$$\"3O+++++++C Fiw$\"\"&F*-%'COLOURG6&%$RGBG$\"*++++\"!\")F*F*-F$6$7S7$F($\"+\"peH=$! \"*7$$\"+Hf999Ff[l$\"+\"fvVF$Ff[l7$$\"+)3#\\upFf[l$\"+SJ#)=TFf[l7$$\"+D(H\"[tFf[l$\"+16\"G;%Ff[l7$$ \"+.e'3r(Ff[l$\"+![*H)>%Ff[l7$$\"+%49G8)Ff[l$\"+#RbwB%Ff[l7$$\"+G)4j]) Ff[l$\"+D^'3F%Ff[l7$$\"+WQ4;*)Ff[l$\"+M@l0VFf[l7$$\"+K^I1$*Ff[l$\"+o9F PVFf[l7$$\"+w'\\Mr*Ff[l$\"+UmzoVFf[l7$$\"+YNt35F_[l$\"+R4^'R%Ff[l7$$\" +!=i!\\5F_[l$\"+LX;DWFf[l7$$\"+'H_44\"F_[l$\"+tgl`WFf[l7$$\"+kxTF6F_[l $\"+UZZxWFf[l7$$\"+q9!o;\"F_[l$\"+$>UA]%Ff[l7$$\"+\"p)[27F_[l$\"+qy&o_ %Ff[l7$$\"++KHZ7F_[l$\"+gR0]XFf[l7$$\"+3h!eG\"F_[l$\"++esrXFf[l7$$\"+F \"o&G8F_[l$\"+L(o\\f%Ff[l7$$\"+.<*pO\"F_[l$\"+!zv^h%Ff[l7$$\"+rq,39F_[ l$\"+Y=5OYFf[l7$$\"+_@>X9F_[l$\"+kB`aYFf[l7$$\"+OV$e[\"F_[l$\"+(R`Tn%F f[l7$$\"+ja2C:F_[l$\"+X`:#p%Ff[l7$$\"+f%[Sc\"F_[l$\"+=-a5ZFf[l7$$\"+p$ HJg\"F_[l$\"+G%G\"GZFf[l7$$\"+hM/W;F_[l$\"+9UF_[l$\"+@>2h[Ff[l7$$\"+$*pre>F_[l$\"++ JVw[Ff[l7$$\"#?F*$\"+mF&G*[Ff[l-Fjz6&F\\[lF*F*F][l" 2 258 258 258 2 0 1 0 2 9 0 4 2 1.000000 45.000000 45.000000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "fff:=subs(z=f,ff):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 116 "plot([[GAMMA(x),x,x=2..5],[x,subs(z=4.5,fff),x=1..20 ]], color=[red,blue]);" }} {PARA 13 "" 1 "" {INLPLOT "6$-%'CURVESG6$7S7$$\"\"\"\"\"!$\"\"#F*7$$\" 3Cyv_,>E!e5F0$\"34+]PC#)GA@F07$$\"3MF KnY6m$4\"F0$\"3%****\\Peui=#F07$$\"3$)eU0XxWL6F0$\"3')***\\i3&o]AF07$$ \"3R(e$Ha8=x6F0$\"3#)**\\(oX*y9BF07$$\"3@#3hmH[;A\"F0$\"3:+]P9CAuBF07$ $\"3q[cQJ.$>F\"F0$\"3\")**\\P*zhdV#F07$$\"3\"[jr9J&yG8F0$\"3,+]P>fS*\\ #F07$$\"3ChVbcYt!R\"F0$\"34+](=$f%Gc#F07$$\"3ovPKLrQg9F0$\"3$*****\\#y ,\"GEF07$$\"3/r)4Ys4r_\"F0$\"37++Dr\"zbo#F07$$\"3q%=\"p#*\\t3;F0$\"3%* ****\\(4&G]FF07$$\"39d,#o*H@)p\"F0$\"34++]7nD:GF07$$\"3rvE2OF>#z\"F0$ \"3%*****\\-*oy(GF07$$\"3;H['fj)o%)=F0$\"35+]PpnsMHF07$$\"3PQp3i&=V+#F 0$\"3,++]siL-IF07$$\"3h!o$p6!pX6#F0$\"3,+++!R5'fIF07$$\"3mP/kOI0aAF0$ \"3)***\\P/QBEJF07$$\"3W0Dym\"3&)Q#F0$\"3\"******\\\"o?&=$F07$$\"3TSW` %\\A!\\DF0$\"3%)**\\Pa&4*\\KF07$$\"3I9-1jto:FF0$\"3&)**\\7j=_6LF07$$\" 3)*fcLbyM0HF0$\"34++vVy!eP$F07$$\"3XyGKx-/&4$F0$\"3K+](=WU[V$F07$$\"3! *pR/Wy#zJ$F0$\"3(****\\7B>&)\\$F07$$\"3l?571(>:d$F0$\"3)***\\P>:mkNF07 $$\"3y6It_bR7QF0$\"3'***\\iv&QAi$F07$$\"3`tUg1um&4%F0$\"31++vtLU%o$F07 $$\"3Fu$o#y'>gT%F0$\"37+++bjm[PF07$$\"3?rp_A+YfZF0$\"3\"****\\(yb^6QF0 7$$\"3XYo#)y4$H7&F0$\"3****\\PMaKsQF07$$\"3:X3JsY?mbF0$\"3=++D6W%)RRF0 7$$\"3.C&)*>UpQ+'F0$\"3z*****\\@80+%F07$$\"3\"ys9g]In^'F0$\"3h****\\7, HlSF07$$\"3c%oa&))=HEqF0$\"3')**\\P4w)R7%F07$$\"3z*HB'zgKPwF0$\"3<++]x %f\")=%F07$$\"3l%*4iO:3p#)F0$\"3!)**\\P/-a[UF07$$\"3EkZwO#pY**)F0$\"3[ +](=Yb;J%F07$$\"3\"=Es@K&Rv(*F0$\"3V****\\i@OtVF07$$\"3$fm<[\\iw1\"!#; $\"3')**\\PfL'zV%F07$$\"3.(Q'*e)oYj6Fiw$\"3u*******)>=+XF07$$\"3O$>ZJ5 .;F\"Fiw$\"3.++DE&4Qc%F07$$\"3L3))RUy7!R\"Fiw$\"3<+]P%>5pi%F07$$\"3SJ \\)Q.M+^\"Fiw$\"39+++bJ*[o%F07$$\"3+@JMZZ'=m\"Fiw$\"34++Dr\"[8v%F07$$ \"3S**RH]197=Fiw$\"3++++Ijy5[F07$$\"3\"o)39Z+:*)>Fiw$\"32+]P/)fT([F07$ $\"3%H1q!=igw@Fiw$\"3^+]i0j\"[$\\F07$$\"3O+++++++CFiw$\"\"&F*-%'COLOUR G6&%$RGBG$\"*++++\"!\")F*F*-F$6$7U7$F($\"+E_A)f#!\"*7$$\"+lH227Ff[l$\" +4HpVf\"Ff[l$\"+KPm9HFf[l7 $$\"+)3#\\upFf[l$\"+\"z)=7TFf[l7$$ \"+D(H\"[tFf[l$\"+9SdeTFf[l7$$\"+.e'3r(Ff[l$\"+?b[&>%Ff[l7$$\"+%49G8)F f[l$\"+[\"zfB%Ff[l7$$\"+G)4j])Ff[l$\"+LR&)pUFf[l7$$\"+WQ4;*)Ff[l$\"+E! 3^I%Ff[l7$$\"+K^I1$*Ff[l$\"+@T*pL%Ff[l7$$\"+w'\\Mr*Ff[l$\"+,]noVFf[l7$ $\"+YNt35F_[l$\"+hAY'R%Ff[l7$$\"+!=i!\\5F_[l$\"+]2:DWFf[l7$$\"+'H_44\" F_[l$\"+tRl`WFf[l7$$\"+kxTF6F_[l$\"+AYZxWFf[l7$$\"+q9!o;\"F_[l$\"+#>UA ]%Ff[l7$$\"+\"p)[27F_[l$\"+6w&o_%Ff[l7$$\"++KHZ7F_[l$\"+Y20]XFf[l7$$\" +3h!eG\"F_[l$\"+(*=rrXFf[l7$$\"+F\"o&G8F_[l$\"+_[#\\f%Ff[l7$$\"+.<*pO \"F_[l$\"+A(y]h%Ff[l7$$\"+rq,39F_[l$\"+f$3fj%Ff[l7$$\"+_@>X9F_[l$\"+U[ ?aYFf[l7$$\"+OV$e[\"F_[l$\"+Z`htYFf[l7$$\"+ja2C:F_[l$\"+$[X8p%Ff[l7$$ \"+f%[Sc\"F_[l$\"+-iN4ZFf[l7$$\"+p$HJg\"F_[l$\"+&Rwks%Ff[l7$$\"+hM/W;F _[l$\"+'z4Ru%Ff[l7$$\"+%f[Mo\"F_[l$\"+J8DgZFf[l7$$\"+LguBF_[l$\" +>\\P][Ff[l7$$\"+$*pre>F_[l$\"+d*\\Q'[Ff[l7$$\"#?F*$\"+KA-y[Ff[l-Fjz6& F\\[lF*F*F][l" 2 258 258 258 2 0 1 0 2 9 0 4 2 1.000000 45.000000 45.000000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 959 "`evalf/invGAMMA`:=proc(x ::algebraic)\n local xx,oldDigits,L,oldz,newz,GAMMAoldz;\n xx:=evalf (x);\n if not type(xx,float) then RETURN(`invGAMMA`(x)) fi;\n oldDig its:=Digits;\n Digits:=5;\n if xx>=20 then\n L[1]:=ln(xx):\n L [2]:=ln(L[1]);\n L[3]:=ln(L[2]);\nnewz:=L[1]*(1/L[2]+(L[3]+1)/(L[2] ^2)+(L[3]^2+L[3])/(L[2]^3)+(L[3]^3+1/2*L[3]^2-L[3]-1/2)/(L[2]^4)+(L[3] ^4-1/3*L[3]^3-5/2*L[3]^2-L[3]+1/6)/(L[2]^5)+(11/6*L[3]-23/6*L[3]^3-17/ 12*L[3]^4+L[3]^5+5/12)/(L[2]^6));\n elif xx>=1. then\n newz:=4.5; \n else\n ERROR(`evalf/invGAMMA expects its argument to be\n \+ not smaller than 1, but received `,x)\n fi;\n while evalf(newz, Digits-1)<>evalf(oldz,Digits-1) do\n oldz:=newz;\n newz:=oldz-(G AMMA(oldz)-xx)/Psi(oldz)/GAMMA(oldz);\n od;\n Digits:=3;\n while ol dDigits>Digits do\n Digits:=min(2*Digits,oldDigits+5);\n oldz:=n ewz;\n GAMMAoldz:=GAMMA(oldz);\n newz:=oldz-(GAMMAoldz-xx)/Psi(o ldz)/GAMMAoldz;\n od; \n evalf(newz,oldDigits)\nend:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "Digits:=50:" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 10 "xxx:=19.5;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6# >%$xxxG$\"$&>!\"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 56 "start :=time():\nzzz:=`evalf/invGAMMA`(xxx);\ntime()-start;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$zzzG$\"SN&)GM3P#pL97S]MXW:_t(yr`pg[!#\\" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#$\"$&**!\"$" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "GAMMA(zzz);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\" S+++++++++++++++++++++++]>!#[" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 10 "xxx:=20.5;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$xxxG$\"$0#!\"\" " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 56 "start:=time():\nzzz:=`e valf/invGAMMA`(xxx);\ntime()-start;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6 #>%$zzzG$\"Sprf6))Hw8no911#4=_BT8R%p^%*[!#\\" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"$I'!\"$" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "GAMMA(zzz);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"S*************** *******************************\\?!#[" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "Digits:=10:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 68 "plot([[GAMMA(z),z,z=2..4],[x,invGAMMA(x),x=1..6]],color=[red,blu e]);" }}{PARA 13 "" 1 "" {INLPLOT "6$-%'CURVESG6$7S7$$\"\"\"\"\"!$\"\" #F*7$$\"3K57202A>5!#<$\"3ALLL3VfV?F07$$\"3!RA'>RCDP5F0$\"3smm\"H[D:3#F 07$$\"3I$==UtF!f5F0$\"3WLL$e0$=C@F07$$\"3LeGp]))f#3\"F0$\"3RLL$3RBr;#F 07$$\"3LV:6S?w26F0$\"3%om;zjf)4AF07$$\"3+&3!y_pmK6F0$\"3VLLe4;[\\AF07$ $\"3dt?3.S6g6F0$\"3-++Dmy]!H#F07$$\"3SH;yF![.>\"F0$\"3>LLezs$HL#F07$$ \"3b^;7tlVA7F0$\"31++D@1BvBF07$$\"3xb[BL*pvD\"F0$\"3\"pmm;_M(=CF07$$\" 3ji*Q!3PQ!H\"F0$\"37LL$3y_qX#F07$$\"3HoM9s!=&H8F0$\"3'******\\1!>+DF07 $$\"3sPlzmcEr8F0$\"3w*****\\Z/Na#F07$$\"3m\"zRt[URT\"F0$\"36+++NfC&e#F 07$$\"3&fg/MX\")[X\"F0$\"35LLez6:BEF07$$\"3h(QXF0$\"3)om;zihl&HF07$$\"3pGO@2$z\")*>F0$\"3JLL$3#G,**HF0 7$$\"3v1S!*>8#>3#F0$\"3i.E,VBF0$\"3$***** **pvxlJF07$$\"3z4sZu0lUCF0$\"3z****\\_qn2KF07$$\"3*)GWIV*RYa#F0$\"3%)* **\\i&p@[KF07$$\"3G%fS!=TskEF0$\"3#)****\\2'HKH$F07$$\"3wR::zM:zFF0$\" 3`mmmwanLLF07$$\"3Z%H==j%f3HF0$\"3u*****\\2goP$F07$$\"3zIp4:&oF.$F0$\" 3CLLeR<*fT$F07$$\"3v2(>%zm[wJF0$\"3(******\\)HxeMF07$$\"3\"GbY>\\&y>LF 0$\"3Bmm\"H!o-*\\$F07$$\"3n*=\"[\">z%yMF0$\"3))***\\7k.6a$F07$$\"3$*\\ S.k[*Hk$F0$\"3emmmT9C#e$F07$$\"33>EDQ5!e#QF0$\"3\"****\\i!*3`i$F07$$\" 3)y#>B\\Ww7SF0$\"3;LLL$*zymOF07$$\"3P#oy&\\Lw:UF0$\"30LL$3N1#4PF07$$\" 3_6(f/=b'HWF0$\"3jmm\"HYt7v$F07$$\"3j?cSMI'zj%F0$\"3&*******p(G**y$F07 $$\"3p$y$3j;\\\"*[F0$\"3Vmm;9@BMQF07$$\"3ZWbu>7^K^F0$\"3/LLL`v&Q(QF07$ $\"3\"*z&fKeL`S&F0$\"30++DOl5;RF07$$\"3c\"4H(=Pz#o&F0$\"3/++v.UacRF07$ $\"\"'F*$\"\"%F*-%'COLOURG6&%$RGBG$\"*++++\"!\")F*F*-F$6$7U7$F($\"&++# !\"%7$$\"+')G\\a5!\"*$\"+yUb:@Fi[l7$$\"+rd)*36Fi[l$\"+O_'=@#Fi[l7$$\"+ R(*Rc6Fi[l$\"+?Q5&G#Fi[l7$$\"+2P\"Q?\"Fi[l$\"+7h+^BFi[l7$$\"+SwX58Fi[l $\"+Yq]zCFi[l7$$\"+x%3yT\"Fi[l$\"+f5\"*)e#Fi[l7$$\"+&4\\Y_\"Fi[l$\"+Iv `$o#Fi[l7$$\"+CSqB;Fi[l$\"+U&f:w#Fi[l7$$\"+m'pis\"Fi[l$\"+Q(3X$GFi[l7$ $\"+*>VB$=Fi[l$\"+0n3.HFi[l7$$\"+`l2Q>Fi[l$\"+uOmlHFi[l7$$\"+0j$o/#Fi[ l$\"+:B&\\-$Fi[l7$$\"+_>jU@Fi[l$\"+Zn]tIFi[l7$$\"+j^Z]AFi[l$\"+I!)fCJF i[l7$$\"+)=h(eBFi[l$\"+>:csJFi[l7$$\"+Q[6jCFi[l$\"+DP*f@$Fi[l7$$\"+\\z (yb#Fi[l$\"+2xJ`KFi[l7$$\"+b/cqEFi[l$\"+\\zM&H$Fi[l7$$\"+*fa$Fi[l7$$\"++#pxg$Fi[l$\"+s>jtNFi[l7$$\"+g4t.PFi[l$\"+PR%of$Fi [l7$$\"+!Hst!QFi[l$\"+7Y3@OFi[l7$$\"+ERW9RFi[l$\"+'p!GXOFi[l7$$\"+KE>> SFi[l$\"+>9#[Z&Fi[l $\"+$*=YERFi[l7$$\"+(G!e&e&Fi[l$\"+B$[E%RFi[l7$$\"+&)Qk%o&Fi[l$\"+If!o &RFi[l7$$\"+UjE!z&Fi[l$\"+4AerRFi[l7$$\"+60O\"*eFi[l$\"+m\"Ha)RFi[l7$F dz$\"&++%Fe[l-Fiz6&F[[lF*F*F\\[l" 2 258 258 258 2 0 1 0 2 9 0 4 2 1.000000 45.000000 45.000000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 16328 14682 0 0 0 0 0 0 }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{MARK "69 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 }