Complete Port of integration toolbox to OCTAVE

classic Classic list List threaded Threaded
1 message Options
Reply | Threaded
Open this post in threaded view
|

Complete Port of integration toolbox to OCTAVE

Vinayak Dutt


Hi Octavers:

I have finally finished complete port of the integration toolbox of
MATLAB to OCTAVE. I think most of the functions are working OK.
I get some unpredictable behaviour on the mutlidimensional
quadrature gquadnd(). Once in a while Octave bombs with
error message like "unable to assert" in file strbuf.cc which I
can not fathom. But I have been able to get a lot of my integration
work done with this toolbox. The Gaussian quadrature is also
lot faster than routine quad from quadpack (even though its a
script file!).

So here's final release of the integration toolbox. As usual
uudecode | gunzip | untar to unpack the scripts.

Hope you also find it useful as I did :-)

-Vinayak Dutt

PS: For your info: This toolbox comes from MATHWORKS
ftp site: ftp://ftp.mathworks.com/pub/contrib/integration

------------ cut here --------------------
begin 644 int.tar.gz
M'XL( !8+-3   ^P]:6_;R)+S5?H5O0\P1-F4+%*78S\NQHDS$P-)YD@R\X"L
M9]"2*(D)K_"P*6-__%95-\FF+-ER["B;C A#%KNKJ^^JZCI:CI_8LX@G3N ?
M_O"%GDZGUQGV^_"_TS$&XG]GT!7_Q?-#9] Q^@;D=C'?, >F\4/_ART\:9SP
M"*J\G(2WPKU[<U<G.QW3',I.F3]\(X^CS/_YZ[/G__D2\V_ =/=ZZ^:_WS>[
MY?P;!LY_%U;,#YW=_'_QQWKX4S_W)W;&ID'$E-54?PS,==])CMF_)TYDCY,@
M6OQW_7'0LM>I9T?.F+OLO&PR>QL$[BC(ZJS.7IV^?7GZ-$^ASAFM,YV9^,']
M"?-;9ZNQ8.FW<YL%D3-S?,B#<BP*TL3Q[9A=V1%DC1(.;Q,VC0*/O7[^]N7Y
M4T)*N75V%3E)8OMLM*BSXGD17/%HPOYTW!@K*9XS.^11XME^PH(I>^Y#I3:T
MRI^Q5_9XSGUG'"O0[WSGTHYB)UD@]*G+1]SC2OY3Z.S1L&,.CY3$MVD,G0R"
MF.MY$=;M]XZ&K27 7^>!;S.STV?=WE'+&!A#)?.YQQV7\<DDLN/XF+WX\_SE
MFU]>LQ_9NU/CCU?M=Z?MYV?O\M$#&.I/@M\K@[=R<)Y&B['-?H8!\NU(2?^=
M+]B?;?;"CF">/K*7?!3 +,'4V.J@_)I&D]16QD;)^Q,;\I)/^<*&:G5V_IKU
MAD\Z@^5>'[.N,6SUGO1@3+J&DOL3SXY9F3D\NC$FD#T3+?_1'OOMD%K3MB>I
M'(S89APZCD,QLP$*UA3 A &D*QNN&*9CAN7D\RGEDUG;$R^M%GOAS.:,C\=I
MQ,<+&.70Y6.;%@\N\=_>G9[10L0O1TR#I=NLH#)S9( *MD*E_@ &CW&&6Y7[
ML]3E$7R?.<IBE1B@$X!D/8:\DRM*^TK]_NK2F#Q?A';46M&2QR$?GA./*V2)
M5<@(+%N/)W-F9R$N=4B*'ZGB"8Y!VSL&"L;.:++P*<B0NR@&!'8+]Y46L#1&
MDH!EVLKZ@X8S2Z#2&L:A9AQD?W6;#9WKHR;C81@%F>,A-EU94[<_?#.XPXW
M_I<^#?S8$+Y5/FSRGWM4 D]W4W#-8 <L8ZQ9?ZR>CNH2ZJP-6S[P>,R8>?@D
M;]'$3BCQ1UA7ML_;GI,0@1 KXNBA2^)HS9HXTAH K[6RO\RF6!R:<6@V'[!
M-ED=A_D7L[[)9% +83::&T$KJR/;J(!XH-OUC=I-*Z-^]W1_UFPG$0^OY61S
MF"1ZUQ9ZUL24J0.4&[D$< O.4I#*(GH%8G]ILP7 3[6L*:<=,["X?1TX$R2V
MJ6MC]1F 71)98^/ 1SDEAP;L=FC#!S"+2QXY?.3:;,(3SL( UEC!4!'6=HFI
MQ$@+,P8""0A,0<)&D ,+%E=E' +CF6"9HF%%Q9+O3U-_3 3UDKLIK-YQ$,'2
M#0-_0DT*$"BOM*C0\6'X\^&IC@55 OFB,XD= 5X8ZRF)8J,%E:NS/POYXIB=
MI1SDF1?<CVW7X[Z^)#Z]0BE.9YK9&39;_2.C99J]0?U;$?UOG/^0K7T))<#&
MYW\0'+M]D\Y__7YG=_[_&O,O18PMGO^[0],HYM\<TOF_WS%[N_/_%IZ"QKX_
MO[!H[K6?4E_GALY-G7=UWFO6]^I[0ARH?9:TN5>[EYRY)PO<)D/L;29A;@*F
M2IB;P2])F)M70A+FAN"EA+GW2#T=%2.[/!77-!6-:SD=RJ1>.<F\@AOY+_!H
M9+X SR(44%#(@-EO9 VF(6L%CLI3-VFV-ZE.3P)7K5*6 ?QXDG2!,('H J=W
M*14$;IN=3]D5CT@L\6"Y\1FNOLB6)<7A,;$G.H@W"RG&YDUY8]NT*'4A]>HR
MNP[9RY+87FVM% ;"#W0)^GM2AP99AMWJGM2=*=.8SZ,9"!C_9EV8N-K<=D-&
M>ZI>B^PDC?PZB$_.M )K6:PK)-?,<WR+&R?RA6<6-T_JMAO;2_ ](>I2JA/'
M2:1QH\F:0OZ&QBE($*.I8NS2"R*]I4XQ VX!38TNRFQ4@RC?.Y$])AI"SS-8
M*S,;#\B1D[$@M$D?$^.,2QGNRHGM,@-*G@)%LI* QV/'0>K4/*F?Q@8KDAK[
M#4HRE:1#D=15DOZBI+&G)+4I*0J5)$9IB14'4:*]1Y%:H_HMJ+*I5][-I?=N
M\P**OG=F?A#9^NL+*W:N;2U1TE[)M%/1"V[!,.+WS'I_<5(G#:IE'+^& 92I
M]$_'3XT?:XGF-%M&LWEQ4LP_9<ET]E\6P_[E!Z$JCBC4"V")H9C25?6I<#2'
M0'@M3#PP:%:AM?55K7R%A2@]MA-<G)0+O27Z@]*U(MG#=CM,O/#P4^9!6MLC
MY1/R%M+V0 G<*W;FQ(G6$# -V%CQ(DYL3_M7Y+'6M(K@7\U\R4U][ME6HY(+
M._;4"L8)O[1#V-5:\T1Y>4^PQPW]%*=Q&D9 E:8:H='?-XHF_Y*^#=\EED"I
M-738#GJC^3]^8U4I"=S0:80:)P1T;A&?S3NDXR;2<?,0,3RI*VTZA=<[>DMS
M(1M7K]_W]+%:_CMZ5 'P+OEOD-M_2/X7\I]A[.2_KR/_'3VR 'AT0PZXEU;I
M#G%P$UE0T2?M;2!ZE?JD3: K^J1-"I3ZI+U-VBWT27>!KA#N\E&^EJ-\78SR
M"DFO4@%*?0^3].ZL^_YBWZT2'/ 6!A+<8PEP_UX2WX[ND-_8XPAP._EM)[]]
M5_+;T9<0X+Z@8+:I;"@$N*-'D^"J(MS7U?_-> H<W$\>40*\0_XSNL-"_NMV
MA/ZWW^_O_'^V*_]=<E<?A?K5],+*UX V14D0N+4O9, [8*KR!"^I!+FY<&01
M(RD8@A3AM\AL!) _(RHR0+&KN3.> T\$DL''PB<!S4ONP@\\A[L #.+ !"H!
M,<#<]UM&FST+_+$=HJG)9R$'UO2D"T  V7B%_']"-JB5GCH-]"'98[_.'=<)
MV1F_=&*B7S+A=SYR_ "$H6NV<&P73J-D.1IQ8&_"XB5]=IS9')K*T705MVF@
MH.= =ENL'(*1C=TN!@B1P: )(-6)P75 2L'A\XFL0VX0H>T*^@"#X:4N9YH4
MNW"8S Y*J6C'(E@RE.7F,XD6!XWF# '$4*_N@BH#RM[4V9)>"R67FF^9G9Q(
MIY9F',,\--N'\2?@N)JY+Q.:[;],^ ?,\W(,B^;"@DJTB<-G6JI#^H'\"H  
M @W\>"%X-J6/0DR^FEKF_N58,_2/S0:@.V&9I?&#4?/0/-!&+0[_]T<A$'=K
M:D/7:2EFS7V9=8+#85U-M>-F8Q\_E\[*W\H6W2[]ER;FQZSC+OI/-A])_XUN
M3]A_!IT=_=\F_5]R)]@3$O_-U,=P,Q"'B0=X&NQ]CJO!WH-]#?96.AN4 W5?
M=X,]/+<^V-^@C@.9(94[H?XMQ%?&]N"4]M%F<1K9,#)(Z!?D2!@%5[+W<=W%
M$J[MST!470"%=+/R'85^)/^0!B<3=]%D]1J=R;4&K8GC JL'NPA''+L:@X M
M,8"H+[@$CDZGW=_7V/L.6URP _9^P3H7K,G:^\!>WF?XTL+,#!)/"#Y./8UQ
MH/WF,5:]H][;H?_ET>71F, =]+_30YD_E_^%_7_0Z0]V]'^;]#^$+:<<6^'@
M+\A]F79,[^O@WM@)B*KL9_B'5."79V]/_WC>B-G+7T[/?CU]^Z(@Y 0-9,4/
M@!C-4O(]!DEV!H0."!YJZIB3,*%O*XOK0MFQAWFHERCQ FV&1@#6VLP-1D!B
M\YQZ395=@8(9(+O6L/TYQ$D-"!JBK=7R) MPG4 JDJV3?P#!4?>_[R1?) ;L
M'OY?W4&_)^)_#'/G__45YA].TPF*-H]Y KB=_AM&K].7\5_#86\@[']=Y .[
M^?_BS][MH5",J/4#@J&P^/.)DX#(^P<*LJ@-@N*"/4 >_#TL6&I/#0A2-!C5
M<"DE8Y-X*05\=<"4 J!$3"FI=X1,*9 K8J:4W,V"IHIQW"AL2FU\)6Y*R;@M
M<$IM_(W(*25S;>C4<N^78J>4[!7!4\N#LSYZBM8>CQ<H(J2QK01(D?E220&I
M LU5T G5IBT,H?6]6JW&V-R9H0%4:.-0_0\0<+"*X;"7P RI"C<VY[$L!0=;
M&)B9/6FUL 656J&D"Y/*%D'*KKB/IE99"BVOSG1!2 OUGZHE+.M'BVPBBSED
MA8-C&%0U:0J#J8S_JGU>Z%>!PLR1;![TI9:EB*_:QA%?15$_K_9S8KW4U?S\
M].S5\S8P.$0&B*)@DE+H%@#S"1Q8IXZ43'$7J43H^"8>Q%$A6@H./-BO)3TS
MX5N-Y<\5QPFJKQ#$:=$Y?I@F4ELK9_>FPC9FFA> G.H(?4)E+Y>+"5:C/X&&
M<1BU$&O \$$>P2D=UGRNGRAJIP64ERS5U50KA_X@ !"V,0QT(L[[HI&Y8D84
M5-L8R!K@$']C+ 8;#09G V$K4-9]!14JE02F9WG3XAMM4QNEEHY2O^T&,RS]
M$TX@56EGW M=I'4"]+5]M;PJ'K"[RL)C41A14UM;S^;V:!'/[<M5H9UJL\>W
M=UI!='OW_?6(H--)X+>>!92R"HW<+^;JD7D8T1C?LZQT5I8+D!G"'F&TLC8Z
MH<B7!;ZPV/$_8<DE?C43%:]8DZ:Z)F4+8H!&9=H-8B/1##X'SU(G5JW[1R"K
M,+CW*_IH8UN$_MXY+(7>MV@\5:Z]!7Y'G [6G3"6Y13&\7V[<@6!] JG<HK>
M53%,2DFIJ*O*9O.R$@H6>(G?;3=OD*!\Z6Q(A(@P%&MN!;9BL.Z#,&^T.N W
MZ,9M+;V5<A3HU[9]O$';[UV%[(4D-_XZ]KR!R.!O+#/,!+8BW+RR8OT;*_8V
M%[VUM>9=E;M,VL;E!*]F=N.OW"QE[E8WD+:A< :B!J(P*NP?/EEC0( HAY;<
MA?(WFN!W($WC86N)GXR#%'U2<&)S\1DE$4RE=>6GWDB:R?/^HQN2="R\MJ/@
M[]4H@%,'T>TXL)7YD)!^DHK%4M-):/'PA*L33SSQ.'+"I+HR0<@P)S!'J%C)
MJ6XN8] @K)( 2B&$:L;2L-M 5J%QA2Z@C BM59%CSBN>P(F3]=H&(\/>&Q#W
MQLSHE'C^SKD^0'LMDEMQ-.<V'/5N<OYE[D-5^>,5+:G@WK I^0T5:YI"V:4
M5-3_:;:N?@7A^NKSR<IGNCI=&3%E,57BNY8U+>1V%:AQ.K(E)17?)50ERBLS
M^Q+$[.=8^J(CTRP,KD0F?=4R?4$ %=8YS_+[.\2&FF>(!LUCQ'XK-"M;"'3P
MI4!F'BR6&@ZY1@%G"$ XAL<:.3=FS66<9@%L"F#B^ 7J*C@PY;F$QZ\:@'=O
M )A+$(ME"+<$<!&@HP*XXV@L >@K K2H3;W6C08A1%H"IUK>_@+VGVQ+7-;_
MBA/NX]9QA_W/[&.L[Z!C](9FOS\@_7^OVS%W^M\M/'?<A/6@B[!V]V#M[L&Z
M596[L2;W,Z_!VMV"M>86+/9.**-90%$]I1I\6FJCN51#.^-UJN< 5E@!<X>>
MNLW8.9W;ISRF8[2HZ:,?7 D5-KV2]Q8' 2V92X3KM-Y"H^\D3%S"Y2ZP L)9
M9U+3SF"17Y:[1E=U]%*YZ+JH+I_@9B6UGZK7AZY ]^Q+FH2),YW"ID/:@LU2
MM?X F[I)7*CZV]^4.+&:_[?)N_O3=OA_X?\I^#_=_]H?=+H[_K^%9X^=CN(D
M0F=[[14<CN:P^73CR9..\.KY^?3=FS>_H6_Y[SF-0 KY+/" [F(DD%!5+%M\
MD7-- ]<-KI#_%N0%R!$ZL>-!6XH+ZG8DR^1X";,"T&:G%3L%>@PQ\F5'MZ1U
M#OE(<\APX?AC-YW >0W03((Q.1[)"CF)"?E&)F%$6.JDR[W *X^9>;9P\Q\7
ME(F%47#I3$C9D4M-4'>1"BT&>I/B8">VC"4H6)"0A$@'*;0MHO?E:*CC4*\8
M60LQX)]E\OZ6B.PW1/^EDNQ1Z[C+_]/(_?^-X;#?H?BOWM#8^7]NXY$^DZ_?
MO?K[V2_O7K^MBX ?&;%9)#>:&+=LH!A<I%G%MX/<85P[QK =&99: BYEDTOX
M;NO]_]S_N3)^>_O?Z Y+_[\^[GMT!35W_G_;> H1QK\,7&;EQAA&L70@O]F1
MZWAZ&H;BBX_9(!JF>.T6G-?O4^JDOE?SR3)7V%DF#@@E=$$RB4>Y*0L! 3'"
M5B(:_591H)0=78"&:DF)GOH8<Z-5O5D4-Q89CP-BV1A#()/(YLDJ"QD&L-,%
M%\6Q3X8%48P1R(-JAK =QROQ0(D0-10D_I6%'!]6G3\6"(N@F=)#OI;1.,DP
M&B54:1S ^1,$17%E1C[2JV$IM["1+UNV\\E9799RUY:E^5S3PG)R9ZJ9&NUG
M9!QD-L<(7U5-Z.<,(N\/QB'!P=NU_)-Z9J'Q+M9\'6-)ZW+!Y4;&VQ:<#CF(
M1<_T!IFD&KL8HKOIO[#=/6H==]%_B@T0_O](^DG^&PQWY_^MTO_W,JZ?%H#F
MT?$_I_(W\^"TAU2TM)/#035%6_&R#?OV:'E$5,:)<G<6P&EV[HE@(#S2QH%G
MCX().N6RQAMX>1H$'QLZEJ-8_^Z@#X?-,9_8GC-FO^*)46?&DR&DCN"LKKI'
M<C:FJVNP+&!SZ$=+"NH-O?Q@.0<&-=*S_ .C+7II%WXIJ'Q&QG U_?N#%3J'
MY/;"92QLI9\ . H!: QT2S,_M(SF/L";AYZ,H"6U;TDHU_C$:.KE6DPS/<"3
MZT4+9E/>C2"N3" OR[SERXA_*Q6HU8N["L_E/!HVB'*G'J-6$YX/T&Q,.@?V
MI6F%!QC>\-1D4^T:/B?7#$@SPZ%A;U*/ 2/6;@Y!;EUNE9@_6.3T_1-% GLA
M3_A(.(&04K;44PNEB-*+4K=[SD;P(F\Q2(I9DPL-.Q@(1IMZGG0.BZL]W+!W
MLF?4R7VV00>Q<YCP=,E#5[]I;M'9FY#4,L:3)]UVO?[!>F\<>Q<-8%M74XR[
M0QN]!XP0:@X==L@\R!F%M,Y@@>S+5D &O'@;Q<RNI/_F5N7_3D^5_[M#<?]?
MK[^C_UNF_YD^"A= Z+/%!<KT8B4PC82I3,A4BV6^L%&9$XS.1.!+A+S,&8D
M@ES,6V#>HII')6MY+1>69\?S6>1,-$*%900 (,7""@!60_@D=@7!)/ JY9>*
M8W:E-/;-@H_V/J1\=R&A:_8_&2>WM?^['47_)^]_&IB[^S^VOO_%]M3+S5AN
M:5@0=U&"^Y5^"$W8G>&^V/Z7+HF/6L>=^G]#X?\#J?_O[OC_5O?_M56ZH*H[
M>RF]1@8BV)[7Y*>ZVXO?V?Z7[KO;W?\J_Q_VI/YGQ_^WNO^SN95[8ZO;OYI<
M@]?N;L]_[_O?W.K]?Y#=+?;_(+?_]W;W_W^5_6^N(0!F20%V7/_[WO_N5^3_
M _'[[[W!8+?_M[O_72L/MJIL_TIR#5X[N]W_/>__Q__YQ[OV?\]0[+]F1^[_
M_F[_;W7_T\6UE@C-%GX\F1M<Z1G>>Z2CBW02BZO?]=S"N'$):0\+HDU+Y5Y&
M]RM5MF[9.HV7@Z?%%02JL9-^M"@/3<>"^3T1Y=WUPJ=H9"=7MIW?S\ZP=K+]
M4@OPIP/)KKS:L=TI/-\!<1Y,O^3H+JS22A03NF%(^_,\<)4;A!WT/LH=S9.
MB=YC^3@=Y5!Q89B^$1M$_C<**-4<*T$P>"V <&BGNT>$T36_?%Y>T(]N.)=V
MG'O?YR;K4;C* _]J*GOB Q+N.M?YA0"R]J6@(0PP().Z1(9>3^0OE9 S?\M
M7ZH#HQCTTFT NV%C^/U(F/RK!E_QNT+",SZS1BNO1/A@^>RCY2FSX6ITD?/^
M)&-X:_3$V >,XG+BJZGVH;F/"X0;!Y/]CP>0.0HAK2E^.RJS^)I*#*A$O>?L
M:FY'HL-1Q95 Q]YCVDVW"?EXT":Y_D4T9GE]\RALBIL+"N )9(K+\3T=>H+=
M.32A''[CK8F1^PW<[K0A;GN>.I&X\[GXP4TL+?PLR'TN^#_VKJTY<219GU?S
M*[0/#@0&6:4K-$<;L;L3<V(V8F?/7N;L@\/;@0%CVD;07(PTO_YD9E5)I1L&
MMT?N[A7A"RIE9=U455^FLC*Y&Q/5I$ :DN!,3C7Z8CGYJW1.);W"]OACV,-A
MY4XD[E]TZO0P?IX)]]\P6=>;V?-BM=\^Q:G/JJE165Q<*$JN CV*Z( G29(D
M\> *#]B&QA<*6'Z28V6*]RR:>'18C((GH=MPX:Q%$SXD^ ,.=4A/=^1/=/2T
M7PSE0(9J+)"<CJ_T!Y9:&Q@9*P,<)^X@3+(H<1(6Y]U_/:66&:D92]8W.YHX
MSF 1P1OC;6+C4QPVHY6+,>&AE3D-?8#!'5+S(_69D6;D=^M@,]L^C-<S>-Y[
MZ9./UH(P>9*;L#.(F[ ?D2DA309:OK4^K><P*V@BC6"B6'QNC$1%,I-J!$U$
M+^GZ%*-/=+IDE7'#Q#9TV]&NX)9(Y4\PN^UTNCK:<@B242MBG(-1S@((+"(P
M7F2$SM\I2)PE+K!\;$]_:D$?Y(.5^11>ZOX9'>&G03/@N0528;B?OTFWJ+</
M>.=P7ZCOJ(4;=(#-T7%I/#P;W?OG3F-O^6WB_WKM?QS+,E/\[]C<_L=L[/]K
MQ?_<KEH,?QYAQ_@UIJ]WK'=G]0XLHR4X.;-B+X00::7XG#V51]X(X<W!/N[A
M1^&^M/N,96J<" &G2 "IT>NX>,A5!?YJ<B5*CW"KYYUQ!E3G\1OQ  (,";4!
MAJ6(]6281QRO(H:7*)G$DDH0'_52WL"F',P3+QUQ$6'6-!W3X@[OV^- GTM
M9X+]"(!X&%?#_3^K:#^*]4^]1X[XM0AP/L U_1$CK[\*[E-O%,%]3R$,<4,7
MPQSD+88EEA(N7*,^=X6 <%5E$4L6\8LLXG(648*4^APG65J7GILKD7PEDT$6
M&:_7^&@G4PK'5&6&M=!IPO1CA5F,S"CY*C[*+%:9X:/RYUSMNAGN#@[B>+*Z
M6XP3=X^"K5&'N).UCQ2.'4L,J?B1>2J$LN+Y&;X4*DXS#NA\CN03Z9=BMUBB
M]?/=;#)&KQ?B6'F)4T5\SE9A1A#2,'">D12^G<',G):5#ISO]]!&0)(+\N$A
M&3T))II.<?KP1/SLP.O4R<@Q9QD^RWRX5O/P2QDQJ+W5]&I!J)-X/$3_?$"4
M.-93W8=(N_&<[)' 8SX^P1U%R>6C%-Q90A91AC0H']#2H+\#XHS3)M X8WSJ
MX3OQI<4 L#5+CL].%]NUWO[79@4M^SF9MC^18^@_")EUVT:>F=C$K<N_C-?X
M.$:MSU&0G;:CUEHF725)4? Y B$DNEI'HS1WW/H<!]EY"KEETE62% >?8\@=
M7ZUCR%TF5K3XQHYK*(D&N)60<-"%<C_'C8CP->)_K[[X'Z9E6G:*_UUA_]^<
M_WHW_.]5H^\JU/]2%EB/ _AH_P P^;\$0&%72"]^6.W1]<E/$JSS!?M' 60Q
MISS4I"+]*2JGEN19AKPXKYY@413!]]X*ZTL-ZBNQ/L?Y/$#@_%2LG\7Y'-A"
M5[U&':^"[[Q2O@S2?[E:7D'QR(0GISB^".2WN\U^@KLR.<[!$:43Q[S=J0Y2
MZN$\T>;)&(/JBC-]=+QZ/)FLELO5% ,3C-/#UM04\B4F.BU!)J\$)HF.E3P!
MR3@&BKKVC\8I&MM6TJ)1*X7X^$5<Q^GUB4!2[K5G"N+-%MSL_[G]O\[S/[:=
MQG\S;8O'_W6;^%_OMO_3\./B$<4]V!1I!<'_):K  ]<&6BHN> VGXNDA>LU"
M"L+7LRV>.+JL5!:F#C>G'(PD@1[XGI7ZPA#*IF*XB>4XQBTI<2E>'C@"Y%<U
MW@1Q4GS<:?_8P^XV+K@3H1:1=VS\HJ?!=/Z/^@4/+N-'_M=0C87H1H-]19O&
M!7:(9\C7]NK F9%22I)1W_)8R;VT*-'7. Q)MN0&U2]0JU=Z-CY!(JIBA^LG
M[E?[D.^J@F_IJ;(@C8Y1/%(F5&MB/Z47XGPKE9YC*]R3)#HT@5L43=BVAYX!
MUWB!X<&,X]H:/.GVK@H;48'O5V>#4,]*'0&=I<,Y,*G#.>NQ*M'F#*4VYSE1
MYR!#H<(AW<ZSHMRABX/UQ<H=T8[D932MDZ=I:@0V-%%-@QX1 O:!<X-"I!;G
M65]TN"9'X^)/P)4YZ>+:P?;@(J#<H06W0YF>(TTNTF)]5O G4]9YSHA7"?]>
M:3KV$I2O=8&+..0)RW\KI3&Z6J,OJ@/_>;6>_[ MEMI_VH[/SW^X#?ZK%?^I
M%I9>N84E@;PR74[>'7!6<2/I-=@,=@_:O[@>XD>AA_A)N.75!%G1MK2Z-B\H
MA(C)6>J@W]3.<_ZUV7F^6JM4IY&G<+5\HLXHU1<ML_HB9)%3&8V?</>EKDAK
M<8X>B;<W'>"C[Z<Q>^H2^UNP2N50^B6;T5Y>LBBS2"6@@<@UU;R=IRO+N.)*
MGCF,[T>.D6XT9OC,MAQG:-G^D/GFM&^R$91D&Y[I>\SU;=,9,'M(Z89A)/$4
M',.##/;0=ES?\H8,"6[)=/#F<#^"SKRQ-4MCMYW;4:OTA3'4 1X?J$-_:$ -
MO*'+',NT83<3=>A[AL<L<V@//,?S+,_-U:%O&7"+#=G ,P<V&_JB"G?K )#J
MJ ^#I];A"^P6"U:+TH)/M5D4EH9H:"CM#*69H?ZADS4NO'\.LG:"&0/!A#O4
M6M@'=O7#,W!I=^^?N3/@!LQ]A?BO=O^_?HK_/&[_9[$F_E.M^"]4]&LO^_\M
M> ">-QZ &P_ 7[,'X'GC ?BD];]^_[^697O)^@\(B<O_9G/^L];U/W?(AEMM
M%!./._[]4G>_/XR?%SSCW\<@N(*LM/N57/_^9;9[6$VW0I L#3?8[KW@$ACE
M#\#4RF*"*!\0M@;"VR8 T+P,[A>1KH=7K'.-\';&@K!+EZ/6<ADXW64? /TN
M0.^ONM,-KRQ$R?8'Y\,25BPM# .=]>&'78>=:WW0#;LA8MQH%81A%YVT[N *
M)5^0M3Y@H2@!K!^7+ "VZ\<@6J%2E&3CQ\#Z((.^[1C<,;KK1R1:LV#'^IBI
MKXLOG>O'JQUIG37)C],&2$[I0IE*@=##@!E]9 C%@42'321VO QHQW0=!E-F
M7 ,I99Y:D*!;F,'HPLW^C!&E0F$CA<,)@/I*MX"H0]0JF8-D'B>SD<P39%:&
M;@\5AZMU.-*> [S'+_#60]#?&R"97.F&V]U#5OWY"A.>C>XSWD&VQK5N4\F=
M#FF9M37PNWH JBF6"1D?,".OZ,.U31=8'<."6\##H;R\RNL ,STD]$ENFR<C
M,;#HB-H]]-=8V9$&@QZMKAX2/38\>=##?4BXCZ!W^Y 5NQ&Z'NM@=63M9'U$
MY9RD+-VPJ&11.R@0GEX0J!@(B!C'%EJ-3DKO47=.+ROTY=6RH_T>MNA.ZP*$
MR&4G,).37$3Q.TWG1 &(Z*@>N%@&^(@+FD_XF"[AN6:? IP%_4^?.BB6ZI#0
M"5 NI020D2D!%0]PS:7ARD-Y7&"O/I;7;,_OO_^_E_]GV/]=[O_!959C_UGW
M_E]J8G:6_^?YB?Z?YT=\O<X;_\_O/__?Q_\SSG]+QO_TF_E?]_PO\^ \_R+_
MS_/3_3^?L28T$.&WF_]1_,;2_PG^WVQE_W=,DO]MK]'_UCK_?PU@Z$N\/\M4
M;CWY:T"142(0.*YBE#H L'-7T/RZF9O?_OQG;[T G#7_>?P'QW&:^*_USW]6
MO@"P9 6 N4ZO=K>+7V=ZA$J("UP4V%74-:_B;N,;\CN8_U;]\S_5_P/P%_._
MP?_USW^K?/Y;ZOPO ("+RP8 ?#?S_R%Z<_A_NO]GFO\\_B?S&_Q?\_Q_B&!'
MS\Y]2I+S?K^D:=_(W]_S_)=F%37Z?V#,3?P_V)[+[;\<UL1_J77^HQ7XV58U
M=#QK_S3++!M?Q@E7FQ9G8=+!=S) X%:FDHCGUBDO&9%&41"-R-^-+$#<[,N2
M);$AO>#D"*]*"2-NA/0Y0C>=W$%.^!2HI8MWK72!KU$9=Z\I#U@]T4OY*!+4
M 6>H+_CKXJ2K=MH5OD9==/",DV)1&Q&=>(,LO7.>R/OY_''H,RA2'8K*.B)S
MM6K\=%I*BD>UHN9E[K>W_C]--I/:XW^P5/[SI/__1OZK=_W_-:"AU_/BGTP5
M*+!/XI\%.+ ?-V#P^YS_^_>*_TGSWVOB_[S?_-^7SO^]GM7^--/_NYW_X>3M
M#<!?]/_GI/.?XO[!5]-M[+]JG?]2TN(/@+[L9(P[LK=&13MP[J:"G[O^>7;8
MK<+^GU;D6.<$DW#R.#H.M8$XC9QED+AMHY/A2RUS?BA+*HK1D_,^6*2VW#WP
MD[O) 9WUZBD.5\O%^$F;15"1IYB[0PG1:)9U)#W)6C.T35]V.G3Z\X]G^D61
MSDF6J6R(QSK[;,1NN6>0X(9_3WR,$*FEDIH98@>NKNU\!EO-T&= 0+]*/GN$
MU]=./J>3RVEA>? WS>N/;+BT\*\/#-P\!S?#P;YVD8T[PE^\4"HQ'/GNR#7Q
M![ZPX>TU<PH5\C+L+&@&M<>D%EF95CEL9#%O9/GP8]$_N'*PE9:99^MGV+K7
M/E;51][^"'_Q E.5=KO0:ZX/2;8%)0\'0_Z'+ND&4-Q>#SS'3!R[B/)^)_W]
M2E<O/ZS@N7X,5X?DD#V>X*6H,.C)/1@8FO9WFD/D-F"+C#"U+34 (%^K?>SP
MH8(_#O6,0V.&Z6G]L:[N8# 8]8?68,1,9^B-^HX+M17?,9D(@!*'@OFNM/UN
M]O7_S/V?.Y&:U&G_;=$]OO\S9@K_OZ;;[/]U[O]<=R>&_XCGT-WJJ43;*_.U
M?]R'[4H_P&JXCY/S)25F 0=Z 5DL%T_C#8*'O_WRAQ]@SZ;_ ]JVK;*3PB*R
MP=U^AU[1]EOR^T<GUQ9C !$/L[MX^S![5MV.J4Y+MI.'V1+=M4 S+A;"TU[^
MB')Z*%FT%)HF"7^4/H>/'%!.\D7"11_DXP=Y(WX<-^]*1=.W#ZO]TQ3Y]"/1
MU9*'\!H(//B!WG(>"7V<+S,^H<PX6V:<+S,^7B:,K^P?^#K;X+%HQ9&++@.2
M=:3[#1HP9;3^5APM#OS4X"\R*(,H-G&GJ'Z2Z"[H15'7V35)NZP?@Z3;R2=&
ME"B]2'?0SZ)TL]@OXTY.^%OB,1;.[+2)[$1LY!*=YDC$.(L6VYW>AAYI=Q!.
M</ (EP&;]57XA]T'\.;F-D\@/,;)H]7D*D2X%@XLX>N#CEP5CG%,A'=AJX<2
M/O3BQ]73-#C?J_!D%0(HGL\P9K!\;>&WR%TDE&+]6U_0T<J2@R2R!OB*@A\'
MIQ7C%740D&Q\M]5%0_KPOZ/]-R5!=W7IDF.UI,+B-.7=9C9^5 &8[ MZ_2$]
MU%$!,BN.A4GL./!+?)I/%U.*3"@)^WT\"0M(< -" W3)X^PI;G?JAE\5^W^=
MYS\LWTS>_S*+^<+_;[/_O]_^+]SKO@8"8%:^F[?1 23D$M\P:_M$/' RD[<&
M!U G])V6.@/.4FYF<XH&6@(@E/VH=?DR?OBG&IL,JHOS_T/!@4CI+B7<\^;<
M //^(_>_RL9ET":EY=,CD9[@D0HOP0D>**N = %<"Q@2!6K<QTN4=3^<](>X
MF^N:4U%-#0CF-]G^@^ +-_^B,]Q)Z@RW% 0<7R#*W7F?!P>.UJD"%KRJ5G4#
MA"#XFN!!^?X_KU?^]Q+_7\QF_/PG8\W[_W?;_^>OW/SGKY3_YU^%_/]-2?V_
MM81?AS3_\VJ'3OG'.Q%N?4,N>W%AA$40]L?]8K<0KGPY\,IK(5J7J1;B;K-Z
M1&?[:^XM=SO#,H'_>C&;S+8P>G\(I]ERTO452H&G8JPM5WAKO=ZLUIL%9E8+
M4YZ.M%3R[=^ZU+<\:$,:GUZ!@_?D:!\C""SN%Q-M%Z_IJ5&6=W*\_QL!@S=3
M"LS?72DP_V:5 E_7GG_"_E]K_!]  .G^[SO"_T-C__-^^__KY?_Y6\C_\V]0
M_L\!B-;EEPC^<M,]+N=GPOV\),K+?;^1W!NT\HYHY35:#-A;3Q;R.3HY'62<
MS#@%';2/W]]CH!M>,ZVOB>,0_%D+,BABI*(-S(@P@Q.>B"]$^[]%>-%\OO)/
M&?Z;U&O_;;NI_2?\^L+^N_'_\4[XKV#] 2M9;[<93V:]->NMK=[:[JV=WMKM
MK;W>VN^M![WUL (0%BT[*F%?B1$(!W?GD&?J:<"G!!P>5IO'K?8)!I[6P!P\
M1&A7:1CRTIL=&:7Z;6P4JLP1%,L#LCL01@>"/+,S+^8AH 0>QZC5FC^M[J"D
MR9V5?#U8;[A39QCA0*BL\#I0K3)Y2L)+O9]<)[YH6Q1SCK<1 [%L9SL 3]L=
M6=H"@II(A_W)'C:+9I./G""X:<=!NX>/=5N/VK>C5KA??DQ@WC;@<?WZ;@H6
ML@2X&:OLE(M>N[=N]V 4+;C6%QBB11Z*K,S1&;5OTT8L <-A*Y3B/HTGJ[N
M1WGI1TD$OTLNE>QX*$L>?W)W0/A8\-%?$E(R,SSP#*CZHQNX[L'#(-ZRZ.)D
M!W5\H,--@$M=JM455@=NXAE5I5&*2@@=!0 KHQN++.*,K!AO7F#TD5]B3!G$
M.K&XCNE:%%V%XN8P.H'L<@'F>#S!$*,H!-H-M*_=0[+;!'[QAA-%)S72@($%
M)AH\(3<R"XSHY""_R@Z1UWSH*!_U &06YV0EGV22M36>14NY56260(^S@-I3
MG\M,N:[G++(,HJ0?2\9%X8LPD@^/9*Z,4I&Q[+K,R*5C=R.^C+!PT3 YCC?B
MRRB6]\28?I.JMV(OK)]6.UVTOR?:"D/US<'I,OPWKQ?_.9Z=^'^Q3-=MSO^]
M+_Z;OR'^FY^._^;GX;]Y??A/4:94J/5^P5!"&=P5;G>S\10WY0S[#ZW+"UA,
MQ-F/\62RA_K&(E[@ N,9ZL2'U_$.UF&A!4JJWZ- .KSM,J(A*8&BQ9*"1-_%
M& >1E&7;V68QVV)XY0O</>\!:H1QHAE56L+#3<]7 DP!2PH3J?^$K0H?M:<9
M1JZ\N,C'3:>LAP5"KWB- 3F>8IF9*TEYK/2+BPIDLDMB^FVU)0X(AF=,XET;
MF!/ ]%9HP,;/V#P1H#-3^^V*A]K&/OX?*F\K8UV/Y^,%=-J1*-N('&%?097I
MQ04>T0GQJ1$]IG9]IM444'RQG%'<ROT2>LXX2T>8X/Z7E8*MRU/5@IIVLEI0
MU4:27O"85A!5?@+7I+)#M>B@O:'LH+V%\-"("U\H+N2D!1 64%:8E\@*9XL*
M[RTIG"DH0(Z;5%!0Y(3Y67*"E@@*"3/M-6+".5+":;)!(QK\UJ+!)?7O0D:^
M$:2==Y$9RO!_6+/]GV\E]G\6<]PF_N?[XO^P5 *H /CA.1 _K ;YY[[!#[\R
M\[[O(>:HM"B4!]9>CC]:M%D\-PIIT8KQW%BDS9O\K]+N\-4'$C*6A6'!JD]&
M$E8A'\Q=%?1QI@G<4YDK+_R/L/[FS0-/1P_Y_3^RW/KCOYA>NO]S_Z\.\UFS
M_]>Y_\<!#'W. ;1,(]]?,3EY!YGW A:)<-?8BWP7G\+\G^SO9M-:_7_!;3N=
M_SS^@VM:C?_GNN<_#7UQ"5"2^2)@-U[_OM_YOP40M*EW_IN>H\Q_6_C_:>3_
MNN<_#;V8_UKT;TM.^&[43/C_E/G_ZVRS^D@ K[[XKY:;OO^W7+[_^TZ#_VN=
M__@&ZSY:KPY!^@BH0*#BOGB)\?,O?_GXI[_^\O,_L]J")+G=21U@IIR2VXE^
M($G!%SG-HE/__-_-MKOK-R[#-!W3=]WJ_=\T^?QG/DQ[A\=_]=W_<IOY_S[C
MCW\^OITGZ!?/?XKX?XYON1[W_^ W^*^>S\T=ZQU8XMN9G_SF;_?#U!3XYL[N
M'>PDR:8DIW=PDB2'DMS>P4V27$KR>@<O2?(HR>\=_"3)IZ1![S!(D@:H@[WG
M6EJ]O5V$[9[96R]ZK$>5[?2TR6JSF4UV'\?A-K J2+$1)Y)BXTXDQ4:?2(J=
M<2(I=M*)I-AY)Y)BIQ9(4UK:@Y':DCW;.]*W)=3Y[AU<VT?I\WWL5%8&Q^X<
M]D1?9'^,.C^*MG7M5M8''Y#3^1-U*?]C&=QB!KNR0O@8GE4$92@MXEB&_$/)
MK,&U7UDI?."QC%.+('HO2W^\!)PGY3F.9?#S53I*/<A2N\RZ'E96"&?CZ?R)
M^BS^.(7+<QS+,"QF<(]F8&8VAV4Z@VO&&M#]W7^J\=_;G01X,?X'$_8_C+D4
M]QVH_,;_=RV?^R '' H(X)*3#([1<)+Y>20]-NNS N[2M$OML-@]:#NT0L4
M"6C-,0LGLU;*0UG";F[QAU675Z MPQ6MRU2K,9)-+LE;Q "7I,Y +<7VXR*D
M6:-E5"1)][V.V:",V_RUW#YGN7&# ^*H/8RW9*JSW2GA51+U$!I,;]OHP+*R
M;TKAR)'F5\"7JO95 ;8J^@JP4T5>#G6JZ<LV_2K:B@V\BKQT^S[6DU7;]Y'Y
MTFSXS>?H_K^UIL;3:OX69;RT_]MB_^?Z'_+_X#E^$_^KCL_O?Z]M]J$UI2%O
MM9Y73TP+T ).LPQ/_8!01G?'1V];U7?#X[S#%Y@KJ]71&MC\+C/+[XZ/WG:J
M[X9IZRIO'V%>J'Y%#5Q^UW4,4_U8_*Y7?;?8/T#&<_GYA$$A85S.F!UG/!0)
MKF%G/N)1,:MO%YE6\!!/#+,-6#&&CFL-G8%C>Z(_F%5]N]CCEN%ZGNV;'G-L
M=\CD$VT?+\))<GLF+%"FYP\&OC-@YQ0A1I5YAC]@UM#R?7M@L^%0W/:J;Q>+
M "+799XY8$/FN*X'1./-;"RGEFTPA[E#RW-M=S#TD8=RY >)XF >Q8S,^D36
M<77>XD"5$"$3.?.9X?JF/_1LR_.'SM L*S^Q-) UD#.G+'>Q!D4B6,6FB_$F
M;O#+V^W_X:2F_=^V32?=_RUZ_^^Y_]_>%?6FC03A=WZ%7R)*:SCOKKVV#^5.
MU:F]JY0V5=+3/42M9(@#Z!Q,#"3P[V]F=VW6Q@9#4BKU; B(>&=W9IB9G3'>
M;YVF_C_1_*_]V@?U5!KCP<M\RB#D$4Z9PSR/0K7.R^<# V/E738M6[[-?+R#
MB_F,6O5(;,NA-O5LES.G#@EPYWL8!ACQ*//LFB0^M9CE,Y\X+JG'&-HL1&7&
M&4S2]4DL#TP8XF1MQN#PF$-=<((:)/+3_E9VW<QI0^)5=&SGA#PN/]O'6'X4
M7F<4WK-S*<L1HSC$VSSI@:,0TJ,P&V9/7C?9NZON^J#QCTTV4WHOGWB5J@P:
M4? V8J>O+SL^H<[QXY><.H2>I=Z7'N4FPS2_<GC/1Z0:EQ-&N.>6?F70R-,/
M7QN2ULJR"T-Z/O>Y34!#U"'.@4.6G3J$7D0IB^:,,'Q357A4MM5[\SB!)R20
M1%C$SMZP[>9!9=LLZ6OB__'QOSK_>QB=*/^SJ4VT_(_(_*_!_SQI_B<N%&]9
M,8<:UK*AR#K BI4M^P>3^'YM]]HS6Q[C*]\WI2G[D4B>IM8S.LDNY(CK)<?V
M\I"=AP!7^\>@UE4X7":X -^(PL<P4IC64,T/Q^$M+GG'KGJ&<;VU:+77T#:T
M#>T/IBU4((P1+-P=YAQ:9SVS@H%3;/\T4+\"XH=7(+GBP'UF!<'KI?/Z08Y+
MQW-J>X%T7#X<WW>LP]/QBM[$XQG)_?_] NUW/DKS_V0Y?:G<OT;^[U"F7?^5
MZS^<YOZOTQP2L=XX/S^'US\N/WY^>_7A^O*3<?G>N'IW_??%EVOC_=7E1X&\
M9!I_:F]<DH&K"LSZ8 $)XW0D/IYOW8S>"[MX/\Z9\2%%14)HFLG]+-(@&#&"
M;!)\;CM,Q/G?@ P!5@.H5!10H<2Z$>/"K 8M;@;$,I^(E<+P$:O3-RJ6$(B&
MG<)HV:R"HWU)UH*]5=EHXNZT>+F '/E1H4FF.(N">( 83UMPBL9=@&A&\Y[0
MC^"+YQA3#&ESB>VZU$6=O=5@8A^6DP4BM<:W!DSX4WFK7- RRH\S38A468+)
M3^%*X,K,)[=BNZLLM5="3:9WDRF.-(_B&:A^8:S.+<$Z>ZUT^I LD'=@O$FZ
M?OY$4=JEYSJ0__B<6BUE#*."-91[&.4>H9YG^0Y#@\8-,L(5. 0XR!S!P<"P
M<<Y'8%&@I9D;4^G&Y</(ICE'EK](NY9/9&\.-'&RWASLK:([RS)5ZZW( &F6
MYWIZ8(.Y49J^80HGPW@!86"2BVU2CUIT^R=(IO#?7XV+>(3E--Z(*<<RNK)W
M0^L^G,WE &?&]2*>2;QJ&$2@DHF(@JAPX)2-)?_LWH<& BF\95-<%\\9]WU7
MAO&_XB>87Z)( ($'$-#O9_$<P[:<L3;S)=J_<1L+HN5<3%L*0SMY#*+Y[]+\
M1GGS!J_0G1G9L# ,V!1O<:&4<(KF^:[@R5V">-Q).%]&"F%[OIA$4=4<A0X$
M3H-XZ+,X3GK&YS 9!S.$;E0PZ7$"D]0N<H% /@ZC6:^E<@'!M)8-Y-*!C7QF
MUE0**.7SF6M[+J/4=DSAW@(A3\B#\B$DVSWB[@W"!>A/?A,7(0).&@N,!,&_
MJ.# @"]R%)9=/Y.PCIM,(*=RJ\A4RM6&KU8NB=C0<B QTWB$^(XBP4K*D@"]
M7XL3J/C KJ0H(CQ' 82;.1AK(K:]E-_E? SV)H 3GQ"G?8J:$'&IG84^V<-B
MG,3+T;BM!4(TZ<%:&9_ H2^HHZ=EDC!0LFA*OQ]8_ZD[0D^Q_L=V'9:M_Z=,
MX/]P[C;UWRD.O713L#]OUM^HAEYKW"7Q/09U2#X(SE4K4=9H_UA#"H.+:,Y:
M-X.5.5B;3RKP9KN%/V9[=+='JW7;[&+L%2^* '=\C1-HI[;B++;+SE?U@^/(
M-E,UCTVAT7C5-F^ZI-\E7\T;TL=7VJ=?LY89*/%6PTZKN$1HEW8LU 7=*&?S
M&74C5I[L4$X<L8)<N% #_W9K1[6JT$W:AZ89HZ@:JV\IE>05H^,U%QKF%:-$
MVZ,95M ,2S7CV <8C64R\5=42GFK_6*S/BN523#U'(G4^C681=-=I/$>UQ1$
M.$F!?3/ Z"!*L8VQTXY4"36?4$YX6V>*P1UJ"[H1.V4+P=LCZ"12[^.V4I76
M35%ANTFI9ECZKMQ5!(?K<%WE+X0+H]K+*&TK)SE"2$4KI,PM^^+*H,EV"!3;
M+]I=$*B#_.H?A21 )2(CS4ECOYY-$)4ZM8'A)!E",BYL) EN)Y U4L0V?P7-
MDF^H]&WN"; =#8%.O2\A]-&ZHE<1;\F.G&IXXEU@6D*$SY)X$(7W*N'-$*E#
MXY5 KI87NSH_6&G@,PR4 4:_AO<URWD-XA+T3Z17'*M$K]"+$%[Y1.>E5$4X
M=*^,]E,(G>E5D[[ASTX%09F'UTRWI:25*LIU5J6D*G(Y7L'U-H+D]CK=G8R
M"I1Z\2)Q$00])_]L:#[=#?4]+H.T )/?IT@E("F1#>4^!L&YVFYVTR3/MQH=
I.#V0Z5]H';8S%N6JB6TF,P;U!D468:RFV&F.YFB.W/$?]\RW70"4 0#/
 
end