a0a2a1a4a3a6a5a8a7a10a9a12a11a14a13a16a15a18a17a19a5a21a20a23a22a25a24a27a26a29a28a31a30a32a24a33a7a35a34a37a36a39a38a40a28 a41a43a42a45a44a47a46a49a48a51a50a53a52a55a54a57a56a8a58a12a59a45a52a60a44a49a61a63a62a60a42a25a64a65a46a12a44a67a66a45a62a69a68a60a70a72a71a74a73a75a46a72a76a12a44a77a59a78a50a79a46a16a80 a81a16a59a78a82a83a46a43a50a84a61a63a76a72a52a60a54a86a85a87a59a45a52a32a68a60a50a53a52a55a44a77a59a78a50a79a46 a88a27a89a57a90a92a91a32a93a27a94a96a95a18a97a75a98a23a99a75a100a92a101a57a102a103a98a105a104a106a91a107a101a108a102a74a109a105a110a23a111a10a112a35a113a74a114a77a90a115a91 a116a118a117 a100a115a114a49a119a8a120a57a121a23a122a123a120a57a124a57a124a126a125 a0 a1a3a2a5a4a7a6a9a8a9a10a12a11 a13a15a14a17a16a19a18a21a20a23a22a25a24a27a26a7a28a30a29a32a31a34a33a30a22a36a35a38a37a40a39a5a24a41a37a43a42a45a44a46a33a47a20a48a29a32a28a30a49a50a29a51a28a52a20a23a22a25a49a53a39a54a14a17a55a51a16 a56a58a57a60a59a23a61a19a59a19a62a63a57a65a64a66a59a67a61a9a68a70a69a53a64a66a71a73a72a74a69a75a71a76a72a74a77a53a77a17a78a79a68a75a72a41a61a81a80 a82 a14a17a33a30a16a23a16a48a29a32a14a53a28a84a83a85a33a7a14a53a86a47a39a54a14a87a20a48a33a47a39a54a22 a88 a68a41a69a25a89a70a59a19a78a73a90a41a59a19a69a40a61a9a59a91a77a53a78a79a68a70a77a40a59a23a78a66a71a73a57a92a59a48a64 a93a94a64a73a64a66a59a23a69a75a71a79a57a65a72a74a95a34a78a73a68a41a95a60a59a15a68a41a96a7a68a70a78a66a71a73a80a17a68a41a90a70a68a41a69a53a72a41a95a40a77a43a68a41a95a60a97a25a69a17a68a41a98a99a57a65a72a74a95a65a64 a100a102a101a17a95a92a71a79a57a65a103a87a57a60a98a99a59a19a69a53a64a79a57a92a68a70a69a53a72a74a95a34a104a105a69a75a71a73a59a19a90a41a78a81a72a74a95a65a64 a106 a22a75a24a27a26a30a28a7a29a51a31a34a33a7a22a25a16a107a35a38a37a40a39a5a16a23a29a51a28a7a49a43a33a30a55a32a14a87a39a107a108a43a22a70a39a54a28a30a22a25a55a32a16 a109 a103a17a72a41a77a87a71a73a72a74a71a79a57a60a68a41a69a110a72a74a69a53a103a111a89a54a72a74a78a73a57a60a72a41a112a17a95a92a59a15a71a73a78a73a72a41a69a53a64a113a96a114a68a70a78a79a98a107a72a74a71a79a57a60a68a41a69 a115 a57a60a69a17a90a41a101a53a95a60a72a41a78a117a116a75a101a53a72a70a103a87a78a81a72a54a71a79a101a53a78a79a59a70a118 a119 a120a122a121 a123a125a124a67a126a110a6a9a124a122a127a128a11a122a129a113a130a132a131a134a133a99a2a111a124a135a4a30a8a9a136a122a10 a137a117a138a66a139 a140a125a141a143a142a25a144a66a142a146a145a76a147a149a148a91a150a43a151a17a144a113a152a153a148a155a154a157a156a149a156a91a158a17a152a94a141a46a150a40a159 a160a162a161a32a163a43a161a32a163 a164 a22a25a28a52a20a19a39a54a37a43a29a32a86 a164 a37a40a55a51a55a32a37a162a24a70a14a87a20a23a29a51a37a40a28 a165a43a166a75a167a169a168a47a170a21a171 a56a58a101a87a71a91a61a19a68a41a95a60a95a92a68a87a61a23a72a54a71a79a57a60a68a41a69a12a77a43a68a41a57a60a69a70a71a81a64a63a72a54a71a76a77a53a72a41a69a17a59a19a95a162a61a19a59a19a69a75a71a79a78a73a68a41a57a65a103a17a64 SMA-HPC ?1999 MIT Laplace’s Equation in 3-D Basis Function Approach () 1 , 1 i i n cj j c ij panel j x dSxx A α = ′Ψ= ′ ?"!"###$###% () () 11,1 1, 1 ,1 , n cn nnnnc xAA AA x α α # $Ψ#$#$ % &%&%& % &%&= % &%& % &%&%& % &%& Ψ’(’(’ ( && ’( ’ ’’ ’(’ ’’ && Put collocation points at panel centroids ic x Collocationpoint a172a174a173a34a175a177a176a179a178 a104a105a69a180a71a79a80a17a59a91a95a60a72a70a64a113a71a58a95a60a59a23a61a9a71a79a101a17a78a73a59a63a62a153a59a63a62a128a59a19a78a73a59a63a57a60a69a70a71a73a78a79a68a87a103a87a101a40a61a9a59a23a103a181a71a73a68a67a57a92a69a75a71a79a59a23a90a41a78a81a72a74a95a87a59a48a116a70a101a40a72a54a71a79a57a60a68a41a69a40a64a143a72a74a69a53a103a99a64a79a59a19a89a41a59a23a78a73a72a41a95a17a103a87a57a92a182a43a59a23a78a79a59a23a69a70a71 a71a79a59a48a61a81a80a17a69a17a57a65a116a75a101a17a59a23a64a122a96a114a68a41a78a181a103a17a57a60a64a73a61a9a78a73a59a9a71a73a57a92a183a23a57a92a69a17a90a110a71a73a80a17a59a19a98a184a62a128a59a19a78a73a59a85a103a87a59a23a64a73a61a9a78a73a57a60a112a40a59a48a103a34a118a102a104a38a71a181a62a117a72a41a64a67a77a43a68a41a57a60a69a75a71a79a59a48a103a45a68a41a101a87a71a185a71a79a80a53a72a74a71a185a68a41a69a17a59 a68a74a96a94a71a79a80a17a59a99a98a99a68a70a64a66a71a15a77a43a68a41a77a53a101a17a95a60a72a41a78a135a98a180a59a48a72a74a69a53a64a15a68a41a96a153a68a41a112a17a71a73a72a74a57a60a69a17a57a60a69a17a90a110a72a111a103a87a57a65a64a73a61a9a78a73a59a9a71a79a59a99a64a79a59a9a71a135a68a74a96a153a59a48a116a75a101a53a72a54a71a73a57a92a68a70a69a53a64a91a57a65a64a15a71a79a68a12a101a53a64a79a59a180a72 a77a17a57a60a59a23a61a19a59a19a62a63a57a65a64a66a59a185a61a19a68a41a69a53a64a66a71a73a72a41a69a75a71a149a61a9a59a23a69a70a71a73a78a79a68a70a57a60a103a102a61a19a68a41a95a60a95a92a68a87a61a23a72a54a71a79a57a60a68a41a69a36a64a73a61a81a80a17a59a23a98a180a59a70a118a76a186a174a59a181a61a9a68a70a69a53a64a79a57a60a103a87a59a23a78a149a72a111a64a66a57a60a98a99a77a17a95a92a59a181a77a17a78a73a68a41a112a17a95a60a59a19a98 a68a74a96a58a64a79a68a41a95a60a89a25a57a92a69a17a90a85a187a7a72a74a77a17a95a65a72a41a61a19a59a41a188 a64a63a59a23a116a75a101a53a72a74a71a79a57a60a68a41a69a110a57a60a69a179a189a41a190a181a118a40a191a117a80a17a59a185a77a40a68a41a71a79a59a19a69a75a71a73a57a60a72a41a95a51a192a34a193a194a57a65a64a76a77a53a78a79a59a48a64a79a61a19a78a79a57a60a112a40a59a48a103a110a68a41a69a12a71a79a80a53a59a181a64a79a101a17a78a66a195 a96a32a72a41a61a19a59a99a68a74a96a128a71a79a80a17a59a5a61a19a101a17a112a43a59a85a72a74a69a53a103a174a62a128a59a107a69a17a59a19a59a48a103a179a71a73a68a36a61a9a68a70a98a180a77a53a101a87a71a79a59a99a71a73a80a17a59a85a61a81a80a53a72a74a78a73a90a41a59a107a103a87a57a65a64a66a71a79a78a73a57a92a112a17a101a17a71a79a57a60a68a41a69a45a196a46a118a85a104a105a69a125a68a41a78a81a103a87a59a19a78 a71a79a68a45a103a87a68a50a71a79a80a40a72a54a71a85a62a128a59a111a112a17a78a73a59a23a72a41a197a125a71a79a80a17a59a36a64a79a101a17a78a79a96a32a72a41a61a19a59a111a68a41a96a91a71a79a80a53a59a110a61a19a101a17a112a43a59a102a101a17a77a21a57a60a69a75a71a79a68a125a198a199a77a40a72a74a69a17a59a23a95a60a64a99a72a74a69a40a103a21a72a70a64a79a64a79a101a17a98a99a59a12a72 a61a9a68a70a69a53a64a66a71a73a72a74a69a75a71a122a61a81a80a53a72a41a78a79a90a70a59a99a103a17a57a60a64a66a71a79a78a73a57a92a112a53a101a87a71a79a57a60a68a41a69a125a68a41a69a125a59a48a72a41a61a81a80a125a77a53a72a74a69a53a59a19a95a200a118a111a100a36a72a74a71a79a80a17a59a23a98a99a72a74a71a79a57a65a61a19a72a41a95a92a95a60a97a41a192a47a71a79a80a17a57a65a64a185a61a19a68a41a78a73a78a79a59a48a64a66a77a43a68a41a69a53a103a53a64 a71a79a68a174a72a41a64a73a64a66a101a53a98a180a57a60a69a17a90a36a77a53a57a92a59a48a61a9a59a19a62a63a57a65a64a79a59a5a61a9a68a70a69a53a64a113a71a81a72a74a69a75a71a181a112a40a72a41a64a79a57a60a64a122a96a114a101a17a69a53a61a177a71a73a57a92a68a70a69a53a64a23a192a143a59a23a72a41a61a81a80a45a112a40a72a41a64a79a57a60a64a122a96a114a101a17a69a53a61a177a71a73a57a92a68a70a69a162a192a94a201a46a202a149a112a40a59a23a57a92a69a17a90 a61a9a68a70a98a99a77a53a72a41a61a9a71a79a95a60a97a181a64a79a101a17a77a17a77a43a68a41a78a79a71a79a59a23a103a85a68a41a69a99a71a79a80a53a59a91a203a38a204a32a205a67a77a53a72a41a69a17a59a19a95a200a118a143a191a117a80a17a59a76a78a73a59a23a64a79a101a17a95a92a71a79a57a60a69a17a90a185a64a79a59a19a98a99a57a206a195a105a103a87a57a65a64a79a61a19a78a79a59a19a71a79a59a63a59a48a116a70a101a40a72a54a71a79a57a60a68a41a69a107a57a65a64a94a72 a96a114a101a17a69a53a61a9a71a79a57a60a68a41a69a85a68a41a96a162a71a79a80a53a59a149a64a79a77a53a72a54a71a73a57a60a72a41a95a40a89a54a72a74a78a73a57a65a72a74a112a17a95a60a59a76a207a47a118a143a104a105a69a5a68a41a78a81a103a87a59a19a78a94a71a79a68a185a68a70a112a87a71a73a72a41a57a92a69a111a72a185a103a87a57a65a64a79a61a19a78a79a59a19a71a79a59a91a64a79a59a9a71a117a68a74a96a162a59a48a116a75a101a53a72a54a71a73a57a92a68a70a69a53a64a19a192 a62a128a59a185a72a41a64a73a64a66a101a53a98a180a59a122a71a79a80a53a72a74a71a76a71a79a80a17a57a65a64a15a64a66a59a23a98a180a57a92a195a105a103a87a57a60a64a73a61a9a78a73a59a9a71a73a59a122a59a23a116a75a101a53a72a74a71a79a57a60a68a41a69a102a57a65a64a91a64a73a72a54a71a73a57a60a64a66a208a53a59a23a103a110a59a19a209a17a72a41a61a177a71a73a95a92a97a110a72a54a71a91a71a79a80a17a59a180a61a9a59a19a69a75a71a73a78a79a68a70a57a60a103a17a64 a68a74a96a143a71a73a80a17a59a180a77a53a72a74a69a17a59a23a95a60a64a23a118a135a191a117a80a17a57a60a64a149a90a41a57a60a89a41a59a48a64a91a78a73a57a65a64a66a59a67a71a73a68a5a71a73a80a17a59a99a64a66a59a19a71a149a68a74a96a153a198a45a59a23a116a75a101a53a72a54a71a73a57a92a68a70a69a53a64a15a61a9a68a41a78a73a78a73a59a23a64a79a77a40a68a70a69a53a103a87a57a60a69a17a90a99a71a79a68a5a71a79a80a17a59a180a198 a77a53a72a41a69a17a59a19a95a65a64a19a118a85a100a179a72a54a71a79a80a53a59a19a98a107a72a54a71a73a57a60a61a23a72a74a95a60a95a92a97a70a192a34a71a79a80a17a57a65a64a67a77a17a78a73a68a87a61a9a59a48a64a79a64a135a68a41a96a117a61a19a68a41a95a60a95a92a68a87a61a19a72a74a71a79a57a60a68a41a69a174a61a9a68a41a78a73a78a73a59a23a64a79a77a40a68a70a69a53a103a17a64a149a71a79a68a36a64a79a59a9a71a79a71a79a57a60a69a17a90a110a71a79a80a17a59 a78a73a59a23a64a79a57a60a103a87a101a40a72a74a95a52a68a41a78a79a71a79a80a17a68a70a90a41a68a70a69a53a72a74a95a43a71a79a68a107a72a99a64a79a59a9a71a76a68a74a96a46a103a87a59a19a95a92a71a73a72a180a96a114a101a53a69a53a61a177a71a73a57a92a68a70a69a53a64a63a95a60a68a25a61a23a72a54a71a73a59a23a103a12a72a54a71a117a71a79a80a17a59a67a77a40a72a74a69a17a59a23a95a47a61a9a59a23a69a70a71a73a78a79a68a70a57a60a103a53a64a19a118 a210 a191a117a80a17a57a65a64a91a77a53a78a79a68a87a61a9a59a48a64a79a64a149a95a60a59a23a72a27a89a41a59a48a64a76a101a53a64a149a62a63a57a92a71a79a80a50a72a111a64a66a59a19a71a15a68a41a96a58a198a146a95a60a57a92a69a17a59a48a72a74a78a135a72a74a95a60a90a41a59a23a112a17a78a73a72a41a57a60a61a135a59a48a116a75a101a53a72a54a71a73a57a92a68a70a69a53a64a76a71a73a68a111a64a79a68a41a95a60a89a41a59a70a118a15a191a117a80a17a59 a116a75a101a53a72a74a69a75a71a73a57a206a71a73a57a92a59a48a64a46a68a74a96a52a57a92a69a75a71a79a59a23a78a79a59a48a64a113a71a143a112a43a59a19a57a60a69a17a90a135a71a79a80a53a59a91a61a9a68a70a95a92a95a60a68a87a61a19a72a54a71a73a57a92a68a70a69a185a77a43a68a41a57a60a69a75a71a58a62a153a59a23a57a92a90a70a80a70a71a81a64a1a0a3a2a41a188 a64a23a118a7a104a105a69a180a71a79a68a87a103a17a72a27a97a52a188 a64a46a95a92a59a48a61a177a71a73a101a17a78a79a59 a62a128a59a15a62a63a57a60a95a60a95a7a61a9a68a70a69a53a61a9a59a23a69a75a71a79a78a81a72a54a71a79a59a15a68a41a69a110a68a70a112a87a71a73a72a41a57a92a69a17a57a60a69a17a90a180a71a79a80a53a59a122a59a19a69a75a71a79a78a73a57a60a59a23a64a117a68a74a96a30a71a79a80a53a59a67a198a52a195a200a112a25a97a75a195a200a198a50a98a107a72a54a71a79a78a73a57a92a209a111a64a66a80a53a68a54a62a63a69a111a57a60a69a110a71a79a80a17a59 a64a79a95a92a57a65a103a87a59a111a72a41a112a40a68a54a89a70a59a41a118a50a104a105a69a146a71a79a80a53a59a111a61a9a59a23a69a75a71a79a78a73a68a41a57a65a103a146a61a9a68a70a95a92a95a60a68a87a61a19a72a54a71a73a57a92a68a70a69a125a71a79a59a48a61a81a80a17a69a17a57a65a116a70a101a53a59a41a192a46a71a79a80a53a59a111a98a107a72a54a71a73a78a79a57a92a209a45a59a19a69a75a71a79a78a73a57a60a59a23a64a181a57a92a69a25a89a41a68a70a95a92a89a70a59 a57a60a69a70a71a73a59a19a90a70a78a73a72a41a95a60a64a135a68a41a96a128a71a73a80a17a59a5a4a15a78a73a59a19a59a23a69a162a188 a64a67a96a114a101a17a69a53a61a177a71a73a57a92a68a70a69a146a68a54a89a41a59a23a78a135a71a79a80a17a59a5a77a53a72a74a69a17a59a23a95a60a64a23a118a12a56a58a80a25a97a25a64a79a57a65a61a19a72a74a95a60a95a60a97a41a192a7a71a79a80a17a59a5a71a79a59a19a78a73a98a7a6 a202a8a2 a68a41a96 a71a79a80a53a57a60a64a63a98a107a72a74a71a79a78a73a57a206a209a111a61a19a68a41a78a73a78a79a59a48a64a66a77a43a68a41a69a40a103a17a64a153a71a79a68a99a71a79a80a53a59a67a77a40a68a41a71a79a59a19a69a75a71a73a57a60a72a41a95a162a72a54a71a63a71a79a80a53a59a122a61a9a59a19a69a75a71a73a78a79a68a70a57a60a103a111a68a74a96a30a71a79a80a17a59a67a203a38a204a32a205a99a77a40a72a74a69a17a59a23a95a162a103a87a101a17a59a135a71a79a68 a101a17a69a17a57a92a71a91a61a81a80a53a72a41a78a79a90a70a59a15a103a87a59a23a69a53a64a79a57a206a71a113a97a5a103a17a57a60a64a66a71a79a78a73a57a92a112a53a101a87a71a79a57a60a68a41a69a110a68a41a69a111a71a79a80a53a59a10a9a41a204a32a205a107a77a53a72a74a69a17a59a23a95a51a118 a160a162a161a32a163a43a161a65a160 a164 a14a17a55a51a24a41a33a30a55a32a14a87a20a48a29a32a28a30a49a12a11a194a14a87a20a19a39a54a29a14a13a21a13a15a55a32a22a25a42a146a22a25a28a52a20a23a16 a165a43a166a75a167a169a168a47a170a16a15 3-D Laplace’s Equation Basis Function Approach Panel j ic x Collocationpoint , 1 i ij pa cnel j xx A dS′′?= ! , ijccentr j id i o Panel Area xx A ? ≈One point quadrature Approximation x y z t 4 , 1 in 0.25* ijco ij j p Ara xx A e = ? ≈"Four point quadrature Approximation a172a174a173a34a175a177a176a18a17 a19 a69a17a59a99a89a41a59a23a78a79a97a50a64a66a57a60a98a99a77a17a95a92a59a5a62a128a72a27a97a102a68a41a96a117a61a19a68a41a98a99a77a17a101a87a71a73a57a92a69a17a90a110a71a73a80a17a59a85a57a60a69a75a71a79a59a19a90a70a78a73a72a41a95a46a57a92a69a20a6a149a202a8a2a70a192a47a96a114a68a41a78a185a72a12a77a53a72a41a69a17a59a19a95a21a9a102a62a63a80a17a57a65a61a81a80a125a57a60a64 a96a32a72a74a78a107a78a73a59a19a98a99a68a54a89a70a59a23a103a45a96a114a78a73a68a41a98 a77a53a72a74a69a17a59a23a95a63a203a81a192a153a57a65a64a180a71a79a68a45a64a79a57a92a98a99a77a17a95a60a97 a78a73a59a19a77a17a95a65a72a41a61a19a59a5a71a73a80a17a59a102a57a60a69a70a71a73a59a19a90a70a78a73a72a41a95a128a112a25a97a45a71a73a80a17a59a102a57a60a69a70a71a73a59a19a90a70a78a73a72a41a69a53a103 a59a19a89a54a72a41a95a92a101a53a72a74a71a79a59a48a103a45a72a54a71a185a71a73a80a17a59a111a61a19a59a19a69a75a71a79a78a73a68a41a57a65a103a174a68a74a96a63a71a73a80a17a59a5a77a53a72a41a69a17a59a19a95a1a9a40a118 a19 a96a91a61a9a68a70a101a17a78a81a64a66a59a70a192a47a71a73a80a17a57a60a64a181a57a65a64a67a71a79a68a25a68a174a64a66a57a60a98a99a77a17a95a92a57a65a64a66a71a79a57a65a61a107a96a114a68a41a78 a77a53a72a41a69a17a59a19a95a65a64a15a62a63a80a17a57a60a61a81a80a50a57a92a69a75a71a73a59a19a78a81a72a41a61a177a71a23a22a9a64a66a71a79a78a73a68a41a69a53a90a41a95a60a97a3a22a67a62a63a57a206a71a73a80a179a71a73a80a17a59a99a77a53a72a74a69a17a59a23a95a143a203a81a118a25a24a17a68a70a78a149a71a73a80a17a59a23a64a79a59a99a77a53a72a74a69a53a59a19a95a65a64a149a62a63a80a53a57a60a61a81a80a174a72a74a78a73a59 a61a9a95a60a68a70a64a79a59a19a78a135a112a25a97a41a192a30a62a128a59a180a98a107a72a27a97a50a101a53a64a66a59a107a72a12a96a114a68a70a101a17a78a66a195a38a77a43a68a41a57a60a69a70a71a67a57a60a69a75a71a79a59a19a90a70a78a73a72a74a71a79a57a60a68a41a69a125a64a79a61a81a80a53a59a19a98a99a59a41a118a107a104a105a69a45a68a70a101a17a78a67a98a99a57a92a69a53a103a53a64a67a62a153a59a85a64a79a77a17a95a60a57a206a71 a71a79a80a53a59a5a77a53a72a41a69a17a59a19a95a128a101a17a77a146a57a60a69a75a71a79a68a36a96a114a68a70a101a17a78a185a59a48a116a75a101a53a72a74a95a128a64a66a101a17a112a53a77a53a72a74a69a17a59a23a95a60a64a180a72a74a69a40a103a45a62a63a78a79a57a92a71a79a59a85a71a79a80a17a59a12a57a92a69a75a71a73a59a19a90a41a78a81a72a74a95a94a68a54a89a41a59a19a78a122a71a79a80a17a59a12a59a19a69a75a71a79a57a60a78a73a59 a77a53a72a41a69a17a59a19a95a26a9 a72a41a64a99a71a73a80a17a59a36a64a66a101a53a98 a68a74a96a15a96a114a68a41a101a53a78a107a57a92a69a75a71a79a59a23a90a41a78a81a72a74a95a65a64a99a68a41a69a21a71a79a80a17a59a48a64a66a59a110a96a114a68a70a101a17a78a5a64a79a101a17a112a17a77a53a72a41a69a17a59a19a95a65a64a23a118a134a191a117a80a17a59a19a69a157a101a53a64a79a59a110a71a79a80a17a59 a64a73a72a74a98a99a59a117a71a79a78a73a57a60a61a81a197a99a72a41a64a58a112a40a59a19a96a114a68a41a78a73a59a41a192a41a78a73a59a19a77a53a95a60a72a70a61a9a57a60a69a17a90a67a59a23a72a70a61a81a80a180a68a41a96a52a71a79a80a17a59a48a64a66a59a63a96a114a68a70a101a17a78a94a57a60a69a75a71a79a59a23a90a41a78a81a72a74a95a65a64a143a112a25a97a122a71a79a80a53a59a91a77a17a78a73a68a87a103a87a101a53a61a177a71a153a68a74a96a52a71a79a80a17a59 a57a60a69a70a71a73a59a19a90a70a78a73a72a41a69a53a103a34a192a17a59a23a89a27a72a41a95a92a101a40a72a54a71a79a59a48a103a110a72a54a71a63a71a79a80a53a59a181a61a9a59a23a69a75a71a79a78a73a68a41a57a65a103a5a68a41a96a46a59a48a72a41a61a81a80a110a68a74a96a46a71a73a80a17a59a23a64a79a59a122a64a66a101a53a112a17a77a53a72a74a69a53a59a19a95a65a64a76a72a74a69a53a103a12a71a73a80a17a59a185a72a41a78a79a59a48a72a181a68a41a96 a71a79a80a53a59a99a64a79a101a17a112a17a77a40a72a74a69a17a59a23a95a51a118a25a24a53a78a79a68a70a98 a57a92a69a75a71a73a101a87a71a79a57a60a68a41a69a50a62a153a59a180a112a40a59a23a95a92a57a60a59a19a89a70a59a180a71a79a80a53a72a74a71a149a71a73a80a17a57a60a64a67a64a73a61a81a80a17a59a23a98a180a59a180a57a65a64a149a90a70a68a41a57a60a69a17a90a85a71a73a68a110a90a41a57a60a89a41a59a185a101a40a64a15a72 a98a99a68a41a78a73a59a185a72a41a61a23a61a9a101a17a78a81a72a54a71a73a59a181a72a41a69a53a64a66a62a128a59a19a78a48a118a149a191a117a80a53a59a181a116a75a101a17a59a23a64a66a71a79a57a60a68a41a69a47a192a43a80a53a68a54a62a153a59a23a89a41a59a19a78a48a192a53a57a65a64a15a62a63a80a17a59a9a71a73a80a17a59a19a78a149a71a79a80a53a57a60a64a15a57a65a64a76a71a73a80a17a59a181a112a43a59a23a64a66a71a135a62a128a72a27a97 a71a79a68a67a90a70a68a67a68a41a78a94a72a41a78a79a59a128a71a73a80a17a59a19a78a73a59a76a112a43a59a9a71a66a71a73a59a19a78a94a71a79a59a23a61a81a80a53a69a17a57a60a116a75a101a17a59a48a64a28a27 a88 a72a41a69a180a62a128a59a76a103a87a68a67a71a73a80a17a59a91a64a73a72a74a98a99a59a63a197a25a57a92a69a40a103a99a68a74a96a52a57a92a69a75a71a79a59a23a90a41a78a81a72a54a71a73a57a92a68a70a69a181a68a70a69 a77a53a72a41a69a17a59a19a95a7a203a29a27a125a191a7a68a85a112a43a59a181a72a41a112a17a95a92a59a122a71a79a68a5a72a41a69a53a64a66a62a128a59a19a78a76a71a79a80a17a59a48a64a66a59a185a116a75a101a17a59a48a64a113a71a73a57a92a68a70a69a53a64a63a62a128a59a122a62a63a57a60a95a92a95a7a71a73a72a41a197a41a59a122a72a107a112a17a78a73a57a60a59a9a96a143a95a92a68a25a68a70a197a12a72a74a71a91a80a17a68a54a62 a69a25a101a17a98a99a59a19a78a73a57a60a61a23a72a74a95a52a57a60a69a70a71a73a59a19a90a70a78a73a72a74a71a79a57a60a68a41a69a12a57a60a64a117a77a43a59a19a78a79a96a114a68a41a78a73a98a99a59a23a103a34a192a17a72a181a208a53a59a23a95a60a103a12a68a74a96a143a64a66a71a79a101a53a103a87a97a5a197a75a69a53a68a54a62a63a69a12a72a70a64a30a22a9a116a75a101a53a72a41a103a17a78a73a72a74a71a79a101a17a78a73a59a31a22a74a118 a32 a120 a0 a136a2a1a4a3 a124a122a6a9a8a6a5a128a11a8a7 a0a76a121 a9a10a1a30a136a2a11a110a6a9a11a12a3 a13a117a138a66a139 a140a125a141a143a142a25a144a66a142a146a145a76a147a149a148a91a150a43a151a17a144a113a152a153a148a155a154a157a156a149a156a91a158a17a152a94a141a46a150a40a159 a14a47a161a32a163a43a161a32a163 a164 a37a43a55a32a55a32a37a47a24a41a14a87a20a48a29a32a37a43a28a16a15a12a29a32a16a48a24a74a39a54a22a75a20a23a29a18a17a41a14a87a20a48a29a32a37a43a28a84a37a40a35 a163 a15a184a13a149a31a34a33a30a14a87a20a48a29a32a37a43a28 a165a43a166a75a167a169a168a47a170a20a19 a193a22a21a32a207a24a23a26a25a28a27a16a29 a30a32a31 a21a114a207a34a33a66a207a24a35a36a23a66a196a4a21a114a207a24a35a37a23a39a38a41a40a4a35 a207a43a42a20a44a45a46a33a48a47a6a49 a88 a59a19a69a75a71a73a78a79a68a70a57a60a103a111a61a9a68a41a95a60a95a60a68a25a61a23a72a54a71a73a59a23a103a111a77a17a57a92a59a48a61a9a59a23a62a63a57a60a64a79a59a67a61a9a68a41a69a40a64a113a71a81a72a74a69a75a71a76a64a79a61a81a80a17a59a23a98a99a59 Normalized 1-D Problem Basis Function Approach Collocation Discretization of 1-D Equation () ( ) ( ) 1 0 ,xgxx x dSσ′′′Ψ=" [ ]0,1x ∈ 0 0x = 1nx =1x 1nx ?2x 1σ nσ 1c x 2c x nc x () () 11 , to be evaluated j i jx j x n c j g xx dxSσ ?= Ψ′′=!" !""#""$ a193a22a21a114a207a24a50a52a51a53a23a26a25a55a54a57a56 a2a53a58 a29 a196 a2a59a27a61a60a48a62 a60a62a39a63a65a64 a31 a21a32a207a24a50a52a51a66a33a79a207 a35a23a39a38a41a40 a35 a67 a68a6a69 a70 a204a72a71a8a73a52a74a75a74a77a76a6a78a66a79a81a80a82a78a81a204a72a74a84a83 a172a174a173a34a175a177a176a86a85 a187a162a59a19a71a73a64a117a71a81a72a74a197a70a59a185a72a107a64a66a57a60a98a99a77a17a95a60a59a122a59a9a209a17a72a74a98a99a77a17a95a60a59a67a57a60a69 a210 a190a181a118a53a191a117a80a17a59a185a103a17a68a41a98a107a72a74a57a60a69a102a57a60a64a76a71a79a80a17a59a181a64a66a59a23a90a41a98a99a59a19a69a75a71a87a44a45a46a33a88a47a88a49a47a68a74a96a46a71a73a80a17a59a122a78a73a59a23a72a74a95 a95a60a57a92a69a17a59a70a118a7a186a174a59a63a62a117a72a74a69a75a71a7a71a73a68a122a64a66a68a70a95a92a89a70a59a94a71a73a80a17a59a63a57a60a69a70a71a73a59a19a90a70a78a73a72a41a95a75a59a48a116a75a101a53a72a54a71a73a57a92a68a70a69a181a64a79a80a17a68a54a62a63a69a180a72a74a71a46a71a79a80a17a59a117a71a73a68a41a77a180a68a74a96a43a71a79a80a17a59a63a64a79a95a92a57a65a103a87a59a63a68a41a69a181a71a79a80a53a57a60a64 a103a87a68a70a98a99a72a41a57a92a69a47a118a117a191a117a80a53a59 a4a15a78a79a59a23a59a19a69a162a188 a64a117a96a114a101a17a69a40a61a177a71a79a57a60a68a41a69a179a57a60a64a149a103a87a59a19a69a17a68a41a71a79a59a48a103a110a112a75a97 a31 a21a114a207a89a33a79a207 a35a23a177a118a117a104a105a69a36a72a85a61a19a59a19a69a75a71a79a78a73a68a41a57a65a103a110a61a9a68a70a95a92a95a60a68a87a61a19a72a54a71a73a57a92a68a70a69 a64a73a61a81a80a17a59a19a98a99a59a41a192a30a62a128a59a85a103a87a57a92a89a25a57a65a103a87a59a107a71a73a80a17a59a5a103a87a68a70a98a99a72a41a57a92a69a125a57a60a69a75a71a79a68a36a198a84a64a79a59a19a90a41a98a99a59a23a69a70a71a81a64a90a44a207 a2 a33a66a207 a2a53a91 a29 a49 a9a92a25a93a45a46a33a88a94a48a94a88a94a6a33a79a198a143a192a46a62a63a57a206a71a73a80 a207 a30 a25a32a45a185a72a41a69a53a103a107a207 a56a96a95 a210 a118a46a191a117a80a17a59a149a61a81a80a40a72a74a78a73a90a41a59a76a103a87a59a23a69a53a64a66a57a92a71a113a97a107a196a102a57a65a64a153a72a70a64a79a64a79a101a17a98a99a59a23a103a99a71a73a68a185a112a43a59a91a77a17a57a60a59a23a61a19a59a19a62a63a57a65a64a66a59a15a61a9a68a70a69a53a64a113a71a81a72a74a69a75a71a94a68a70a69 a59a23a72a70a61a81a80a5a68a41a96a7a71a73a80a17a59a23a64a79a59a135a57a92a69a75a71a79a59a23a78a79a89a54a72a41a95a60a64a23a118a143a191a117a80a17a59a135a77a40a68a41a71a79a59a19a69a75a71a73a57a60a72a41a95a51a192a40a193 a57a60a64a117a71a73a80a17a59a19a69a12a59a19a89a54a72a41a95a92a101a53a72a74a71a79a59a48a103a111a72a74a71a128a71a73a80a17a59a122a61a9a59a23a69a75a71a79a78a73a68a41a57a65a103a17a64a128a207a24a50a52a51a177a118 a191a117a80a17a57a65a64a30a78a79a59a48a64a66a101a17a95a92a71a73a64a30a57a60a69a181a198a5a59a23a116a75a101a53a72a54a71a73a57a92a68a70a69a53a64a30a57a92a69a180a198a85a89a54a72a74a78a73a57a65a72a74a112a17a95a60a59a23a64a23a192a27a71a79a80a17a59a117a61a9a68a70a95a92a95a60a68a87a61a19a72a74a71a79a57a60a68a41a69a185a62a128a59a19a57a60a90a41a80a75a71a81a64 a0 a202 a33a79a203a59a25a97a47a98a33a88a94a48a94a88a94a6a33a79a198 a62a63a80a17a57a65a61a81a80a50a61a23a72a74a69a179a112a43a59a99a62a63a78a79a57a92a71a66a71a73a59a19a69a174a57a92a69a174a98a107a72a54a71a73a78a79a57a92a209a102a96a114a68a41a78a73a98a110a118 a19 a101a17a78a149a71a73a72a70a64a66a197a102a57a65a64a15a71a79a68a111a208a53a78a81a64a113a71a67a59a19a89a54a72a41a95a92a101a53a72a74a71a79a59a181a71a79a80a17a59a99a59a23a69a70a71a73a78a79a57a60a59a23a64 a68a74a96a153a71a79a80a17a59a107a98a107a72a74a71a79a78a73a57a206a209a179a72a41a69a53a103a50a64a79a101a17a112a53a64a79a59a23a116a75a101a17a59a23a69a75a71a79a95a60a97a179a64a66a68a70a95a92a89a70a59a185a71a73a80a17a59a85a64a79a59a9a71a67a68a74a96a128a59a23a116a75a101a53a72a54a71a73a57a92a68a70a69a53a64a23a118a22a99a91a68a74a71a79a59a180a71a73a80a53a72a54a71a67a71a73a80a17a59a99a203 a9 a204a32a205 a59a19a69a75a71a73a78a79a97 a68a74a96a67a71a73a80a17a59a179a98a107a72a54a71a73a78a79a57a92a209a84a57a60a64a12a72a74a69a84a57a60a69a75a71a79a59a23a90a41a78a81a72a74a95a91a68a41a96a67a71a79a80a17a59 a4a15a78a79a59a23a59a19a69a162a188 a64a5a96a114a101a17a69a53a61a177a71a73a57a92a68a70a69a162a192a76a59a23a89a27a72a41a95a92a101a40a72a54a71a79a59a48a103a134a72a54a71a111a71a79a80a17a59 a61a9a68a70a95a92a95a60a68a87a61a19a72a74a71a79a57a60a68a41a69a84a77a43a68a41a57a60a69a70a71a102a207 a50a51a81a192a149a68a54a89a70a59a19a78a85a71a73a80a17a59a174a57a92a69a75a71a79a59a23a78a79a89a54a72a41a95a100a44a207 a2a98a33a79a207 a2a53a91 a29 a49a38a192a91a62a63a80a17a57a65a61a81a80a199a57a65a64a111a71a79a80a17a59a174a57a60a69a70a71a73a59a19a78a73a89a54a72a74a95a135a68a54a89a41a59a23a78 a62a63a80a17a57a65a61a81a80a102a71a79a80a53a59a185a112a53a72a70a64a66a57a65a64a63a96a114a101a17a69a40a61a177a71a79a57a60a68a41a69a174a201 a2a185a57a60a64a149a69a17a68a41a69a53a183a19a59a19a78a73a68a43a101a114a78a73a59a23a61a23a72a74a95a60a95a162a71a73a80a53a72a54a71a135a62a153a59a185a80a53a72a27a89a41a59a185a61a81a80a53a68a70a64a79a59a19a69a36a72a85a77a53a57a92a59a48a61a9a59a19a62a63a57a65a64a79a59 a61a9a68a70a69a53a64a66a71a73a72a74a69a75a71a185a72a41a77a17a77a17a78a73a68a27a209a25a57a60a98a107a72a54a71a73a57a92a68a70a69a96a102a177a118a111a104a38a96a113a192a7a80a17a68a54a62a128a59a19a89a70a59a19a78a48a192a47a62a128a59a85a103a87a59a23a61a19a57a60a103a17a59a23a103a174a71a79a68a50a61a81a80a17a68a25a68a70a64a79a59a99a72a102a103a17a57a206a182a52a59a19a78a73a59a19a69a75a71a181a64a79a59a9a71a185a68a41a96 a112a53a72a70a64a66a57a65a64a153a96a114a101a17a69a53a61a177a71a73a57a92a68a70a69a53a64a23a192a25a71a79a80a17a57a65a64a117a57a92a69a75a71a73a59a19a90a41a78a81a72a74a95a52a62a128a68a41a101a17a95a65a103a85a112a43a59a135a69a17a68a41a69a53a183a19a59a19a78a73a68a185a68a70a69a17a95a60a97a85a68a41a69a5a71a79a80a17a59a67a64a79a101a17a77a17a77a43a68a41a78a79a71a117a68a41a96a47a71a73a80a17a59a135a112a53a72a41a64a79a57a60a64 a96a114a101a17a69a53a61a9a71a79a57a60a68a41a69a53a64a103a101a32a68a41a78a48a192a63a71a79a68 a112a40a59a174a77a17a78a73a59a23a61a19a57a60a64a79a59a41a192a76a68a70a69a84a71a79a80a53a59a50a57a60a69a75a71a79a59a19a78a81a64a79a59a23a61a177a71a73a57a92a68a70a69a134a68a74a96a185a71a79a80a17a59a125a64a66a101a53a77a17a77a40a68a70a78a66a71a12a68a74a96a122a71a73a80a17a59a174a112a53a72a41a64a79a57a60a64 a96a114a101a17a69a53a61a9a71a79a57a60a68a41a69a102a72a74a69a40a103a92a44a45a46a33a48a47a6a49a37a102a177a118 a189 a13a117a138a51a137 a0a117a144a2a1 a156a4a3a6a5a8a7 a147a76a141a10a9a76a158a17a141a30a151a17a147a76a158a11a5a12a0a63a150a53a159a13a5a14a1a15a5 a165a43a166a75a167a169a168a47a170a17a16 a27 a29 a30a19a18 a21a114a207 a23a39a38a70a207a21a20 a18 a22 a47 a23a25a24 Normalized 1-D Problem Simple Quadrature Scheme ()fx x0 11 2 Area under the curve is approximated by a rectangle a172a174a173a34a175a177a176a21a26 a24a17a68a70a78a58a71a73a80a17a59a76a71a73a57a92a98a99a59a15a112a40a59a23a57a92a69a17a90a181a95a92a59a19a71a73a64a128a61a9a68a70a69a53a61a9a59a23a69a70a71a73a78a73a72a74a71a79a59a76a68a70a69a99a71a79a80a17a59a149a71a79a68a70a77a17a57a65a61a63a68a74a96a47a103a87a59a23a89a41a59a23a95a92a68a70a77a17a57a92a69a53a90a185a72a122a90a70a68a75a68a87a103a99a69a25a101a17a98a99a59a19a78a73a57a65a61a19a72a74a95 a71a79a59a48a61a81a80a17a69a17a57a65a116a75a101a17a59a76a96a114a68a70a78a153a59a23a89a27a72a41a95a92a101a40a72a54a71a79a57a60a69a17a90a67a71a73a80a17a59a15a57a92a69a75a71a79a59a23a90a41a78a81a72a74a95a40a68a74a96a47a72a67a96a114a101a53a69a53a61a177a71a73a57a92a68a70a69 a18 a21a114a207 a23a94a68a70a69a107a71a79a80a53a59a15a103a87a68a70a98a107a72a74a57a60a69a43a44a45 a33a88a47a6a49a38a118a46a186a50a59 a72a41a64a73a64a79a101a17a98a99a59a149a71a73a80a53a72a54a71a91a71a79a80a17a59a67a57a60a69a75a71a79a59a19a90a70a78a73a72a41a69a53a103a111a57a60a64a76a72a28a27a79a64a79a98a180a68a25a68a41a71a79a80a30a29a15a96a114a101a17a69a53a61a177a71a73a57a92a68a70a69a162a192a17a71a73a80a17a68a41a101a53a90a41a80a110a62a153a59a67a62a63a57a60a95a60a95a162a59a9a209a17a72a74a98a99a57a60a69a17a59a135a71a79a80a53a57a60a64 a72a41a64a73a64a79a101a17a98a99a77a87a71a79a57a60a68a41a69a12a95a60a72a74a71a79a59a23a78a23a118a1a24a30a57a92a78a81a64a66a71a76a62a153a59a67a72a74a78a73a59a135a90a41a68a70a57a92a69a17a90a181a71a79a68a85a103a17a59a19a89a41a59a23a95a92a68a70a77a12a72a180a69a53a72a41a57a92a89a70a59a67a72a74a77a17a77a17a78a73a68a70a72a70a61a81a80a85a96a114a68a70a78a63a68a41a112a87a71a81a72a74a57a60a69a17a57a92a69a53a90 a72a181a90a41a68a25a68a87a103a85a72a74a77a17a77a53a78a79a68a27a209a87a57a60a98a99a72a74a71a79a57a60a68a41a69a85a68a41a96a162a71a79a80a53a59a15a57a60a69a75a71a79a59a23a90a41a78a81a72a74a95a43a68a74a96a47a71a73a80a17a57a65a64a94a96a114a101a53a69a53a61a177a71a73a57a92a68a70a69a111a68a70a69a85a71a79a80a53a57a60a64a117a57a60a69a75a71a79a59a19a78a73a89a54a72a74a95a200a192a25a62a63a80a17a57a60a61a81a80a111a62a153a59 a61a19a72a41a95a92a95a52a71a73a80a17a59 a22a9a64a79a57a60a98a180a77a53a95a92a59a67a116a75a101a53a72a41a103a17a78a73a72a74a71a79a101a17a78a73a59a15a64a73a61a81a80a17a59a23a98a180a59 a22a54a118 a191a117a80a17a59a5a64a79a57a60a98a180a77a53a95a92a59a48a64a113a71a122a71a73a80a17a57a60a69a17a90a179a62a128a59a5a61a19a72a41a69a45a103a87a68a36a57a65a64a67a71a79a68a50a78a79a59a23a77a17a95a65a72a41a61a9a59a99a71a73a80a17a59a85a57a60a69a75a71a79a59a19a90a70a78a73a72a41a95a58a62a63a57a92a71a79a80a45a71a73a80a17a59a85a71a73a80a17a59a85a77a17a78a73a68a87a103a87a101a53a61a177a71 a68a74a96a30a71a79a80a53a59a67a57a92a69a75a71a79a59a23a90a41a78a81a72a74a69a53a103a162a192a87a59a19a89a54a72a74a95a60a101a53a72a54a71a73a59a23a103a12a72a54a71a91a72a180a77a43a68a41a57a60a69a70a71a76a57a60a69a53a64a79a57a60a103a17a59a135a71a79a80a17a59a122a57a60a69a75a71a79a59a23a78a79a89a54a72a74a95a200a192a17a72a74a69a40a103a5a71a73a80a17a59a67a95a60a59a19a69a17a90a41a71a79a80a110a68a74a96a30a71a79a80a17a59 a57a60a69a70a71a73a59a19a78a73a89a54a72a74a95a200a192a87a62a63a80a17a57a60a61a81a80a110a57a60a69a12a71a79a80a17a57a65a64a76a61a19a72a70a64a66a59a135a57a65a64a63a101a17a69a17a57a92a71a113a97a41a118a94a104a38a96a30a62a128a59a122a61a81a80a17a68a25a68a70a64a79a59a91a71a73a80a17a59a122a77a43a68a41a57a60a69a75a71a76a68a74a96a7a59a23a89a54a72a74a95a60a101a53a72a54a71a73a57a92a68a70a69a110a72a70a64a128a71a79a80a17a59 a61a9a59a23a69a75a71a79a78a73a68a41a57a65a103a181a68a41a96a52a71a79a80a17a59a76a57a60a69a75a71a79a59a23a78a79a89a54a72a74a95a200a192a41a57a200a118 a59a41a118a46a207 a25a57a45 a94a31a17a192a41a62a128a59a76a61a19a72a41a95a92a95a53a71a79a80a17a59a91a64a73a61a81a80a17a59a23a98a180a59a25a22a9a98a99a57a60a103a87a77a43a68a41a57a60a69a75a71a153a116a75a101a53a72a41a103a17a78a73a72a74a71a79a101a17a78a73a59a31a22a74a118 a109 a98a99a57a60a103a17a77a40a68a70a57a92a69a75a71a111a116a75a101a53a72a70a103a87a78a73a72a74a71a79a101a17a78a73a59a110a64a79a61a81a80a17a59a23a98a99a59a12a78a73a59a19a77a53a95a60a72a70a61a9a59a23a64a180a71a73a80a17a59a36a72a74a78a73a59a23a72a174a101a17a69a40a103a87a59a19a78a107a71a73a80a17a59a36a61a9a101a53a78a79a89a70a59 a18 a21a32a207a24a23a99a112a25a97 a72 a78a73a59a23a61a177a71a81a72a74a69a17a90a70a95a92a59a63a62a63a80a53a68a70a64a79a59a76a80a17a59a19a57a60a90a41a80a75a71a94a57a65a64a143a71a73a80a17a59a76a96a114a101a53a69a53a61a177a71a73a57a92a68a70a69 a18 a21a114a207 a23a143a59a19a89a54a72a41a95a92a101a53a72a74a71a79a59a48a103a107a72a54a71a153a207 a25a32a45a46a94a32a31a87a118a46a191a117a80a17a59a91a64a73a61a81a80a17a59a23a98a180a59a91a57a60a64 a59a9a209a17a72a70a61a177a71a58a62a63a80a17a59a23a69 a18 a21a114a207 a23a143a57a65a64a94a72a67a61a9a68a70a69a53a64a66a71a73a72a74a69a75a71a48a118a34a33a91a68a54a62a153a59a23a89a41a59a19a78a48a192a74a62a63a80a53a72a74a71a94a57a65a64a58a95a60a59a23a64a73a64a58a68a41a112a25a89a25a57a60a68a41a101a53a64a143a57a65a64a143a71a73a80a53a72a54a71a94a71a79a80a17a59a15a64a79a61a81a80a17a59a23a98a99a59 a57a65a64a149a59a19a209a17a72a41a61a177a71a135a62a63a80a17a59a23a69 a18 a21a114a207 a23a149a57a60a64a135a72a111a95a92a57a60a69a17a59a23a72a41a78a91a96a114a101a53a69a53a61a177a71a73a57a92a68a70a69a50a68a41a96a94a207a45a72a41a64a149a62a153a59a23a95a92a95a200a118a185a191a117a80a17a59a180a98a99a68a70a64a66a71a15a68a70a112a75a89a25a57a60a68a41a101a53a64a149a62a117a72a27a97a12a68a41a96 a64a79a59a19a59a19a57a60a69a17a90a180a71a79a80a53a57a60a64a63a57a65a64a63a112a25a97a85a78a73a59a23a72a41a95a92a57a60a183a19a57a60a69a17a90a180a71a79a80a53a72a74a71a76a62a63a80a17a59a19a69 a18 a21a32a207a24a23a128a57a65a64a76a72a99a64a113a71a73a78a73a72a41a57a92a90a70a80a70a71a63a95a60a57a92a69a53a59a41a192a87a71a73a80a17a59a122a72a74a78a73a59a23a72a181a101a17a69a53a103a87a59a23a78a76a57a92a71a76a57a60a64 a72a111a71a79a78a81a72a74a77a43a59a19a183a23a68a41a57a65a103a34a118a180a191a117a80a53a57a60a64a15a71a73a78a73a72a41a77a40a59a23a183a19a68a41a57a65a103a36a80a40a72a41a64a135a59a9a209a17a72a70a61a177a71a79a95a60a97a110a71a79a80a53a59a85a64a79a72a41a98a180a59a99a72a74a78a73a59a23a72a12a72a41a64a15a71a73a80a17a59a99a78a79a59a48a61a177a71a81a72a74a69a17a90a70a95a92a59a180a62a63a80a17a57a65a61a81a80 a71a79a80a53a57a60a64a76a64a73a61a81a80a17a59a23a98a180a59a135a101a53a64a79a59a23a64a128a71a73a68a107a72a74a77a17a77a17a78a73a68a27a209a87a57a92a98a107a72a74a71a79a59a149a71a79a80a17a59a67a57a60a69a75a71a79a59a19a90a70a78a73a72a41a95a4a101a32a61a19a72a41a69a110a97a41a68a41a101a12a64a66a59a23a59a67a62a63a80a25a97 a27 a102a9a118 a187a162a59a19a71a73a64a76a71a79a78a73a97a111a71a73a68a5a103a87a59a23a78a79a57a60a89a41a59a67a71a73a80a17a57a65a64a91a57a60a69a179a72a85a64a79a95a92a57a60a90a41a80a75a71a73a95a92a97a110a103a87a57a92a182a43a59a23a78a79a59a23a69a70a71a15a62a128a72a27a97a70a118a63a104a105a69a53a64a66a71a79a59a23a72a70a103a110a68a74a96a143a71a73a80a17a59a185a57a60a69a75a71a79a59a19a78a73a89a54a72a74a95a7a112a40a59a23a57a92a69a17a90 a44a45a46a33a88a47a88a49a153a62a128a59a111a61a19a68a41a69a53a64a79a57a60a103a17a59a19a78a181a72a41a69a146a57a60a69a70a71a73a59a19a78a73a89a54a72a74a95 a44a45a46a33a36a35 a49a52a33a37a35a39a38 a45a102a71a73a68a179a112a43a59a111a72a179a112a17a57a206a71a99a98a99a68a41a78a73a59a85a90a41a59a23a69a17a59a19a78a81a72a74a95a200a118a36a186a174a59a111a98a107a72a27a97 a59a9a209a17a72a41a77a17a69a53a103 a18 a21a114a207a24a23a117a72a41a112a40a68a70a101a87a71a63a71a79a80a53a59a122a61a9a59a19a69a75a71a73a78a79a68a70a57a60a103a5a68a74a96a30a71a79a80a17a57a65a64a63a57a60a69a70a71a73a59a19a78a73a89a54a72a74a95a200a192a41a40a207 a25 a205a42 a18 a21a114a207 a23a26a25 a18 a21a43a40a207a24a23a14a44a46a45 a21a114a207 a23 a38 a18 a21a47a40a207 a23 a38a70a207 a44 a45 a21a32a207a24a23 a42 a23a49a48 a38 a42 a18 a21a51a50a98a23 a38a70a207 a42 a18a53a52a55a54a57a56a58a52a55a59a61a60 a50a22a42a103a44a45a46a33a36a35 a49 a62a63a80a17a59a23a78a79a59a62a45 a21a32a207a24a23a61a25 a207a64a63a65a40a207a7a118a132a191a117a80a53a59a50a95a65a72a41a64a66a71a111a71a73a59a19a78a73a98 a57a92a69 a71a79a80a53a59a174a59a19a209a87a77a53a72a74a69a53a64a79a57a60a68a41a69a134a57a65a64a111a71a79a80a17a59a45a191a30a72a27a97a75a95a60a68a41a78a12a64a66a59a23a78a79a57a60a59a23a64 a78a73a59a19a98a107a72a74a57a60a69a53a103a87a59a23a78a23a118a46a187a162a59a19a71a73a64a63a57a60a69a75a71a79a59a19a90a70a78a73a72a74a71a79a59a149a71a79a80a17a57a65a64a63a59a19a209a87a72a41a77a53a72a74a69a40a64a66a57a60a68a41a69a111a68a54a89a41a59a19a78a128a71a73a80a17a59a67a57a92a69a75a71a73a59a19a78a73a89a27a72a41a95a26a44a45 a33a37a35a46a49 a27 a205 a30 a18 a21a114a207 a23a39a38a70a207 a25a66a35 a18 a21a43a40a207a24a23a14a44 a35a11a67 a23a69a68 a38 a42 a18 a21a70a50a98a23 a38a70a207 a42 a33a76a59a23a69a53a61a9a59a15a71a73a80a17a59a67a59a19a78a73a78a79a68a70a78a117a57a92a69a12a71a79a80a53a59a67a98a180a57a65a103a87a77a43a68a41a57a60a69a75a71a91a116a70a101a40a72a41a103a87a78a81a72a54a71a73a101a17a78a79a59a135a72a74a77a53a77a17a78a79a68a27a209a87a57a60a98a107a72a54a71a79a57a60a68a41a69a5a57a60a64 a71 a25 a27 a205 a30 a18 a21a114a207a24a23a84a38a41a207a72a63a62a35 a18 a21a47a40a207 a23 a25 a35a73a67 a23a69a68 a38 a42 a18 a21a70a50a98a23 a38a70a207 a42 a74 a109 a64a66a71a79a78a81a72a74a57a60a90a41a80a75a71a58a95a60a57a92a69a53a59a91a68a41a78a128a72a122a61a9a68a70a69a53a64a66a71a73a72a74a69a75a71a94a62a128a68a41a101a53a95a60a103a99a90a70a59a19a69a17a59a23a78a73a72a74a71a79a59a76a72a67a183a23a59a19a78a73a68a185a72a41a64a58a71a79a80a17a59a15a64a66a59a48a61a9a68a41a69a40a103a107a103a87a59a19a78a73a57a92a89a54a72a74a71a79a57a60a89a41a59a76a72a74a69a40a103 a71a79a80a53a59a19a78a73a59a9a96a114a68a41a78a73a59a63a71a73a80a17a59a15a72a41a112a40a68a54a89a70a59a76a59a9a209a87a77a17a78a73a59a23a64a73a64a79a57a92a68a70a69a99a71a79a59a19a95a60a95a65a64a153a101a40a64a58a71a73a80a53a72a54a71a128a71a79a80a17a59a15a59a19a78a73a78a73a68a41a78a94a57a65a64a153a57a60a103a87a59a23a69a75a71a79a57a65a61a19a72a74a95a60a95a60a97a180a183a23a59a19a78a73a68a122a96a114a68a70a78a153a112a43a68a74a71a73a80 a71a79a80a53a59a23a64a79a59a135a61a23a72a41a64a79a59a23a64a23a118 a99a76a68a54a62a67a192a117a72a74a69a53a68a74a71a79a80a53a59a19a78a99a57a92a98a99a77a43a68a41a78a79a71a73a72a41a69a70a71a85a77a53a57a92a59a48a61a9a59a110a68a74a96a135a57a92a69a17a96a114a68a41a78a73a98a99a72a74a71a79a57a60a68a41a69a21a71a79a80a53a72a74a71a85a62a153a59a102a68a70a112a87a71a73a72a41a57a92a69a21a96a114a78a73a68a41a98 a71a73a80a17a57a60a64a85a95a92a57a92a71a66a71a73a95a92a59 a59a9a209a87a59a23a78a73a61a19a57a60a64a79a59a107a57a60a64a185a71a79a80a40a72a54a71a180a71a79a80a17a59a111a59a19a78a73a78a79a68a70a78a185a64a79a61a23a72a74a95a60a59a23a64a181a72a41a64a67a71a73a80a17a59a110a61a9a101a17a112a43a59a5a68a41a96a63a71a79a80a17a59a12a103a87a68a41a98a107a72a74a57a60a69a146a95a60a59a19a69a17a90a41a71a79a80a162a118a174a191a117a80a17a57a65a64a122a71a73a59a19a95a60a95a60a64 a101a53a64a76a71a79a80a53a72a74a71a76a96a114a68a41a78a91a64a79a98a107a72a74a95a60a95a47a59a23a69a17a68a41a101a17a90a70a80a110a57a92a69a75a71a79a59a23a78a79a89a54a72a41a95a60a64a63a68a70a101a17a78a149a64a79a57a60a98a180a77a53a95a92a59a122a78a73a101a17a95a60a59a67a57a60a64 a19a1a0 a112a53a101a87a71a91a96a114a68a70a78a76a95a60a72a41a78a79a90a70a59a19a78a91a103a17a68a41a98a107a72a74a57a60a69 a95a60a59a19a69a17a90a41a71a79a80a53a64a117a62a128a59a135a80a53a72a27a89a41a59a15a71a79a68a107a64a79a59a19a59a23a197a85a112a40a59a19a71a66a71a73a59a19a78a76a98a99a59a9a71a73a80a17a68a87a103a17a64a23a118 a14a47a161a65a160a162a161a32a163 a2 a42a146a44a7a39a54a37a4a3a162a29a32a28a30a49a174a20a23a26a30a22a6a5a99a24a41a24a70a33a47a39a54a14a17a24a41a18 a165a43a166a75a167a169a168a47a170a8a7 a27 a29 a30 a18 a21a114a207a24a23a84a38a41a207 a20 a47 a23 a18 a22 a47 a68 a24 a44 a47 a23 a18 a22a10a9 a68 a24 Normalized 1-D Problem Simple Quadrature Scheme Area under the curve is approximated by two rectangles Improving the Accuracy ()fx x0 11 4 3 4 1 2 a172a174a173a34a175a177a176a6a11 a19 a69a17a59a58a62a117a72a27a97a15a68a74a96a17a57a60a98a180a77a53a78a79a68a54a89a25a57a60a69a17a90a63a98a107a72a27a97a149a112a43a59a94a71a73a68a91a103a87a57a60a89a25a57a60a103a87a59a94a71a79a80a17a59a128a57a60a69a70a71a73a59a19a78a73a89a54a72a74a95a96a44a45a46a33a48a47a6a49a75a57a92a69a75a71a79a68a149a64a66a101a17a112a53a57a92a69a75a71a79a59a23a78a79a89a54a72a41a95a60a64 a44a45 a33a53a45a46a94a32a31 a49 a72a74a69a40a103a92a44a45a46a94a32a31 a33a48a47a6a49a162a72a74a69a53a103a12a62a63a78a73a57a206a71a73a59a15a71a79a80a17a59a67a57a60a69a75a71a79a59a19a90a70a78a73a72a41a95 a27 a29 a30 a18 a21a32a207a24a23a39a38a70a207 a25 a27 a30a13a12a14 a30 a18 a21a32a207a24a23a39a38a70a207 a44 a27 a29 a30a13a12a14 a18 a21a114a207 a23a39a38a70a207 a72a74a69a40a103 a72a74a77a53a77a17a95a92a97a146a72a174a98a99a57a65a103a87a77a40a68a70a57a92a69a75a71a107a78a73a101a17a95a92a59a111a71a79a68a174a71a73a80a17a59a110a57a92a69a75a71a79a59a23a90a41a78a81a72a74a95a128a68a41a69a21a59a23a72a41a61a81a80 a64a66a101a53a112a17a57a92a69a75a71a73a59a19a78a73a89a27a72a41a95a51a118 a186a50a59a110a68a41a112a17a71a73a72a74a57a60a69a21a72 a64a73a61a81a80a17a59a19a98a99a59a99a64a66a80a17a68a54a62a63a69a174a57a60a69a179a71a73a80a17a59a85a64a66a95a60a57a65a103a87a59a41a118a107a191a117a80a17a59a180a96a32a72a41a61a9a71a79a68a70a78 a29 a42 a72a74a77a17a77a43a59a23a72a41a78a79a57a60a69a17a90a12a57a92a69a50a96a114a78a79a68a70a69a75a71a135a68a41a96 a18 a21 a29 a42 a23a135a72a41a69a53a103 a18 a21 a67 a15 a23 a72a74a78a73a59a17a16a113a101a53a64a113a71a63a71a73a80a17a59a122a103a87a68a70a98a99a72a41a57a92a69a110a95a60a59a19a69a17a90a41a71a79a80a53a64a23a118 a18a20a19a22a21 a176a24a23a26a25a28a27a30a29a41a176a146a178 a88 a72a41a69a179a97a70a68a41a101a179a61a19a68a41a98a99a59a185a101a17a77a50a62a63a57a206a71a73a80a174a72a41a69a36a59a9a209a87a77a17a78a73a59a23a64a73a64a79a57a92a68a70a69a110a96a114a68a41a78a15a71a79a80a17a59a180a59a19a78a73a78a73a68a41a78a91a57a60a69a36a71a73a80a17a57a65a64a15a61a23a72a41a64a79a59 a27 a33a76a68a54a62 a98a185a101a53a61a81a80a110a103a87a68a25a59a23a64a128a71a73a80a17a59a122a72a41a61a23a61a9a101a17a78a81a72a41a61a19a97a107a57a92a98a99a77a17a78a73a68a54a89a41a59 a27 a14a47a161a65a160a162a161a65a160 a82 a22a25a28a30a22a70a39a54a14a17a55a117a28a32a31a34a33a149a37a40a29a51a28a52a20a36a35a7a37a53a39a54a42a174a33a7a55a51a14 a165a43a166a75a167a169a168a47a170a8a37 a27 a29 a30 a18 a21a32a207a24a23a39a38a70a207a61a20 a56 a38 a202a58 a29 a47 a198 a18 a22 a203 a63 a29 a42 a198 a24 SMA-HPC ?1999 MIT Normalized 1-D Problem Simple Quadrature Scheme ()f x x0 1 a39 a0 a59a19a97a111a116a70a101a53a59a23a64a66a71a79a57a60a68a41a69a53a64a63a72a41a112a40a68a70a101a87a71a117a71a79a80a17a59a67a98a99a59a9a71a73a80a17a68a87a103a1a0 a2a102a37a4a3 a35a38a14a17a16a23a20a99a86a30a37a125a20a23a26a30a22a36a22a75a39a27a39a54a37a53a39a54a16a107a86a30a22a75a24a70a14a41a18a5a3a107a29a114a20a48a26 a28a7a6 a5a85a39a54a22a102a20a48a26a7a22a75a39a54a22a9a8a94a22a75a20a23a20a23a22a75a39a85a42a146a22a75a20a23a26a30a37a162a86a30a16a10a6 a172a174a173a34a175a177a176a6a11 a109 a64a91a97a41a68a70a101a110a61a23a72a74a69a36a64a79a59a19a59a185a103a17a57a92a89a25a57a65a103a87a57a92a69a53a90a107a71a79a80a17a59a185a57a92a69a75a71a79a59a23a78a79a89a54a72a41a95a47a57a60a69a70a71a73a68a99a71a113a62a153a68a5a78a79a59a48a103a87a101a53a61a9a59a48a64a117a71a79a80a53a59a185a59a19a78a73a78a73a68a41a78a76a72a74a69a40a103a12a71a79a80a17a59a23a78a79a59a67a57a65a64a76a69a17a68 a78a73a59a23a72a41a64a79a68a41a69a146a71a79a68a45a64a113a71a73a68a41a77a157a72a54a71 a16a113a101a53a64a66a71a107a71a113a62a153a68a125a64a79a101a17a112a17a57a60a69a70a71a73a59a19a78a73a89a54a72a74a95a65a64a99a62a63a80a17a59a19a69 a62a153a59a102a61a23a72a74a69 a80a40a72a27a89a41a59a12a198 a64a79a101a17a112a17a57a60a69a75a71a79a59a19a89a54a72a41a95a60a64a99a72a74a69a40a103 a78a73a59a19a77a43a59a23a72a54a71a143a68a70a101a17a78a143a98a99a57a60a103a87a77a43a68a41a57a60a69a75a71a153a116a75a101a53a72a41a103a17a78a73a72a74a71a79a101a17a78a73a59a153a78a73a101a17a95a60a59a63a68a41a69a99a59a23a72a70a61a81a80a180a64a79a101a17a112a17a57a60a69a75a71a79a59a23a78a79a89a54a72a74a95a200a118a47a186a174a59a76a68a70a112a87a71a73a72a41a57a92a69a181a71a79a80a17a59a149a64a79a61a81a80a17a59a23a98a99a59 a27 a29 a30 a18 a21a114a207 a23a39a38a70a207 a25 a56 a38 a202a58 a29 a47 a198 a67a68a6a69a70 a64a79a101a17a112a17a57a60a69a75a71a79a59a19a78a73a89a54a72a74a95 a95a60a59a19a69a17a90a41a71a79a80 a18 a21a32a207 a50a51a39a23 a62a63a80a17a59a23a78a79a59a67a71a73a80a17a59a99a61a9a59a19a69a75a71a73a78a79a68a70a57a60a103a110a68a41a96a58a71a79a80a17a59a181a203a38a204a32a205a5a64a79a101a17a112a17a57a60a69a75a71a79a59a19a78a73a89a54a72a74a95a30a57a60a64a149a207 a50a51a8a25 a29 a42 a21 a202a12a11 a29 a56 a44 a202 a56 a23a12a25 a202a13a11 a64 a14 a56 a118a67a191a117a80a17a59a23a78a79a59a185a57a60a64 a69a17a68a110a103a87a68a41a101a53a112a87a71a67a71a79a80a53a72a74a71a67a62a153a59a99a90a70a72a41a57a92a69a47a192a34a112a17a101a87a71a67a71a79a80a53a59a180a197a70a59a19a97a179a116a75a101a17a59a48a64a113a71a73a57a92a68a70a69a179a57a65a64a135a112a25a97a36a80a17a68a54a62a194a98a181a101a53a61a81a80 a27 a33a76a68a54a62 a103a87a68a25a59a23a64a15a71a79a80a53a57a60a64 a90a70a72a41a57a92a69a50a64a79a61a23a72a74a95a60a59a181a62a63a57a92a71a79a80a50a71a79a80a17a59a180a69a25a101a17a98a185a112a43a59a19a78a135a68a41a96a153a64a66a101a53a112a17a57a92a69a75a71a73a59a19a78a73a89a27a72a41a95a60a64a15a101a53a64a79a59a23a103 a27 a109 a69a40a103a36a208a53a69a40a72a74a95a60a95a92a97a70a192a162a72a74a78a73a59a122a71a79a80a53a59a19a78a73a59a180a61a19a95a92a59a23a89a41a59a19a78 a62a117a72a27a97a25a64a128a68a41a96a7a68a41a112a17a71a73a72a74a57a60a69a17a57a60a69a17a90a99a112a40a59a19a71a66a71a73a59a19a78a91a72a70a61a19a61a9a101a53a78a73a72a70a61a9a97a107a62a63a57a92a71a79a80a110a95a92a59a48a64a79a64a117a59a19a182a43a68a70a78a66a71 a27 a14a47a161a65a160a162a161a14 a15 a33a30a42a146a22a75a39a54a29a32a24a41a14a53a55a94a13 a13a162a14a17a42a45a44a46a55a32a22 a165a43a166a75a167a169a168a47a170a17a16 a27 a29 a30a19a56 a203a200a198a26a21a32a207a24a23a84a38a41a207a21a20 a56 a38 a202a58 a29 a47 a198 a56 a203a38a198 a22 a203a14a63 a29 a42 a198 a24 SMA-HPC ?1999 MIT Normalized 1-D Problem Simple Quadrature Scheme Numerical Example E r r o r n a172a174a173a34a175a177a176a19a18 a187a162a59a19a71a73a64a111a95a92a68a25a68a70a197a157a72a74a71a12a72a41a69a134a59a9a209a17a72a74a98a99a77a17a95a60a59a179a68a41a96a122a57a92a69a75a71a73a59a19a90a41a78a81a72a54a71a73a57a92a69a53a90 a18 a21a114a207 a23a92a25 a56 a203a200a198a26a21a32a207a24a23a5a68a41a69a134a68a70a101a17a78a110a103a87a68a41a98a107a72a74a57a60a69a162a118 a186a50a59 a68a41a112a17a71a73a72a74a57a60a69a12a77a17a78a79a68a70a90a41a78a73a59a23a64a73a64a66a57a60a89a41a59a23a95a92a97a180a112a43a59a9a71a66a71a73a59a19a78a76a72a41a69a53a64a66a62a128a59a19a78a81a64a153a71a79a68a180a71a79a80a53a59a67a57a92a69a75a71a79a59a23a90a41a78a81a72a74a95a52a112a25a97a85a57a60a69a53a61a9a78a73a59a23a72a70a64a66a57a60a69a17a90a181a71a79a80a17a59a135a69a25a101a17a98a185a112a43a59a19a78a63a68a41a96 a64a79a101a17a112a17a57a60a69a70a71a73a59a19a78a73a89a54a72a74a95a65a64a91a198a143a118a67a191a117a80a17a59a180a59a19a78a73a78a79a68a70a78a91a57a60a69a179a59a23a89a27a72a41a95a92a101a40a72a54a71a79a57a60a69a17a90a5a71a79a80a17a59a99a57a60a69a70a71a73a59a19a90a70a78a73a72a41a95a47a57a65a64a149a77a53a95a92a68a41a71a66a71a79a59a48a103a179a72a70a64a149a72a5a96a114a101a17a69a53a61a9a71a79a57a60a68a41a69a179a68a41a96 a71a79a80a53a59a149a69a25a101a17a98a181a112a40a59a23a78a128a68a41a96a47a64a66a101a53a112a17a57a92a69a75a71a73a59a19a78a73a89a27a72a41a95a60a64 a101a114a198a89a102a9a118a58a191a117a80a17a59a15a59a19a78a73a78a79a68a70a78a153a72a74a77a53a77a40a59a48a72a74a78a81a64a58a71a73a68a181a112a43a59a149a90a70a68a41a57a60a69a17a90a181a103a17a68a54a62a63a69a5a72a41a64a21a20a23a22 a29 a56 a14a4a24 a118 a187a162a59a19a71a73a64a63a64a79a59a19a59a135a62a63a80a25a97a41a118 a24a17a78a73a68a41a98 a62a63a80a53a72a54a71a107a62a128a59a12a80a40a72a27a89a41a59 a16a113a101a53a64a113a71a5a64a79a59a19a59a23a69a162a192a153a71a79a80a17a59a102a59a23a78a79a78a73a68a41a78a180a57a92a69a40a64a66a57a65a103a87a59a110a71a79a80a53a59a12a203a38a204a32a205a45a64a79a101a17a112a17a57a60a69a75a71a79a59a23a78a79a89a54a72a74a95a2a101a32a68a74a96a135a95a92a59a23a69a17a90a74a71a73a80 a25 a35 a25 a29 a56 a188a102a63a57a65a64 a205a27a26a42 a15 a83 a14a29a28a31a30a33a32 a51a13a34 a83 a60 a14 a96a114a68a70a78a91a64a79a68a41a98a99a59a28a50a19a202 a42 a44 a202a13a11 a29 a56 a33 a202 a56 a49a200a118 a33a91a59a19a69a53a61a19a59a41a192a40a96a114a68a41a78a63a71a73a80a17a59a185a59a23a69a75a71a79a57a60a78a79a59a67a57a60a69a75a71a79a59a19a78a73a89a54a72a74a95 a44a45 a33a88a47a6a49 a62a128a59a180a61a23a72a74a69a125a64a66a101a17a98 a71a79a80a17a59a48a64a66a59a180a59a19a78a73a78a79a68a70a78a73a64a15a72a74a69a40a103a179a68a70a112a87a71a73a72a41a57a92a69a179a71a79a80a17a59a99a59a19a78a73a78a73a68a41a78a48a192a53a71 a56 a96a114a68a70a78a135a72a41a69a174a72a41a77a17a77a17a78a73a68a27a209a87a57a92a98a107a72a54a71a73a57a92a68a70a69a102a101a53a64a79a57a92a69a17a90 a35 a188 a198a143a188a17a64a66a101a17a112a53a57a92a69a75a71a79a59a23a78a79a89a54a72a41a95a60a64a63a72a70a64 a71 a56 a25 a198a14a35a73a67 a23a69a68 a0 a47 a198 a56 a38 a202a58 a29 a38 a42 a18 a21a70a50a23a202a18a23 a38a70a207 a42a2a1 a67 a68a6a69 a70 a50 a78a66a79a79a52a204a32a205 a202a4a3a6a5a7a9a8 a104a38a71a143a57a60a64a58a59a23a72a41a64a79a97a135a71a79a68a67a64a79a59a19a59a128a71a79a80a40a72a54a71a58a57a92a96 a18 a21a32a207a24a23a30a57a65a64a143a72a15a61a19a68a41a69a75a71a79a57a60a69a25a101a17a68a41a101a40a64a47a96a114a101a53a69a53a61a177a71a73a57a92a68a70a69a162a192 a25a11a10 a188a46a101a114a112a43a59a19a57a60a69a17a90a135a71a79a80a53a59a63a98a180a59a48a72a74a69a96a102a7a98a181a101a53a64a113a71 a112a43a59a107a112a40a68a70a101a17a69a53a103a87a59a48a103a179a112a25a97a102a71a73a80a17a59a107a98a107a72a54a209a87a57a92a98a181a101a17a98 a72a41a69a53a103a36a71a73a80a17a59a107a98a99a57a92a69a25a101a17a98a181a101a17a98 a68a41a96 a18 a21a32a207a24a23a15a68a41a69a50a71a79a80a53a59a107a57a92a69a75a71a79a59a23a78a79a89a54a72a41a95a12a44a45 a33a88a47a6a49 a72a74a69a40a103a110a80a17a59a19a69a53a61a19a59a41a192a43a71a79a80a53a59a19a78a73a59a122a98a185a101a53a64a66a71a149a59a19a209a87a57a60a64a66a71a15a64a79a68a41a98a99a59a28a50 a42a16a44a45a46a33a88a47a88a49a7a64a79a101a53a61a81a80a102a71a79a80a40a72a54a71 a10 a25a10a38 a42 a18 a21a70a50a98a23a13a12a82a38a70a207 a42 a118a13a33a76a59a23a69a53a61a9a59 a62a128a59a15a68a70a112a87a71a73a72a41a57a92a69a111a71a79a80a53a59a67a59a23a64a66a71a79a57a60a98a99a72a74a71a79a59 a71 a56 a25 a198a14a35a11a67 a23a69a68 a38 a42 a18 a21a70a50a98a23 a38a70a207 a42 a25 a47 a23a69a68 a198 a42 a38 a42 a18 a21a70a50a98a23 a38a70a207 a42 a64a79a57a92a69a53a61a19a59 a35a20a25 a47a14a12a54a198a143a118a122a191a117a80a17a57a60a64a135a59a23a78a79a78a73a68a41a78a149a59a23a64a66a71a79a57a60a98a107a72a54a71a79a59a181a71a79a59a23a95a92a95a65a64a135a101a53a64a149a71a79a80a40a72a54a71a135a71a79a80a17a59a107a64a73a61a81a80a17a59a19a98a99a59a180a57a60a64a135a72a41a90a70a72a74a57a60a69a179a59a19a209a17a72a41a61a177a71a15a96a114a68a41a78 a61a9a68a70a69a53a64a66a71a73a72a74a69a75a71a81a64a181a72a41a69a53a103a146a95a60a57a92a69a53a59a23a72a74a78a181a96a114a101a17a69a53a61a9a71a79a57a60a68a41a69a53a64a180a68a70a69a146a71a79a80a53a59a110a103a17a68a41a98a107a72a74a57a60a69a32a101a114a69a53a68a174a80a53a57a92a90a70a80a17a59a19a78a180a68a41a78a81a103a87a59a19a78a180a77a43a68a41a95a60a97a75a69a53a68a41a98a99a57a60a72a41a95a60a64a16a15a102 a72a74a69a40a103a34a192a25a96a114a68a41a78a76a72a99a64a79a98a180a68a25a68a41a71a79a80a12a96a114a101a17a69a53a61a9a71a79a57a60a68a41a69a162a192a87a71a73a80a17a59a67a59a19a78a73a78a73a68a41a78a117a103a87a59a23a61a23a72a27a97a25a64a63a14a53a55a32a49a43a22 a8a30a39a54a14a17a29a32a24a70a14a17a55a32a55a32a18 a161 a13a117a138a18a13 a17 a5a47a148 a5a162a158a53a141a41a3 a7 a147a76a141a10a9a91a158a17a141a30a151a87a147a91a158 a5a15a0a117a150a40a159 a5a14a1a15a5 a14a47a161a14a47a161a32a163 a82 a22a25a28a30a22a70a39a54a14a17a55 a163 a15 a35a47a37a40a39a54a42 a165a43a166a75a167a169a168a47a170a19a18 a27 a29 a30 a18 a21a114a207 a23a39a38a70a207a61a20 a56 a38 a202a58 a29a21a20 a202 a67a68a6a69a70 a22 a74 a202a24a23 a205a23a204 a18 a21a32a207a40a202a52a23 a67 a68a69 a70 a25 a76a6a78a66a79a80a78a177a204 a202 a71 a56a27a26 a71 a202 a56 a204 a24a17a78a73a59a19a59a15a71a73a68a180a77a53a57a60a61a81a197a5a71a79a80a17a59a122a22 a3a17a14a17a55a51a33a7a14a87a20a48a29a32a37a43a28 a44a94a37a40a29a51a28a52a20a23a16a70a118 a24a17a78a73a59a19a59a15a71a73a68a180a77a53a57a60a61a81a197a5a71a79a80a17a59 a3a185a22a25a29a51a49a40a26a52a20a48a16a15a96a114a68a41a78a63a59a23a72a70a61a81a80a5a77a43a68a41a57a60a69a75a71a23a118 a5a5a28a84a28a32a31a79a44a153a37a40a29a51a28a52a20a180a35a38a37a40a39a27a42a174a33a30a55a32a14a179a26a30a14a17a16 a160 a28a157a86a7a22a25a49a40a39a27a22a25a22a25a16a99a37a40a35a91a35a51a39a54a22a25a22a25a86a7a37a43a42a29a28 a172a174a173a34a175a177a176a31a30 a109 a96a169a71a73a59a19a78a117a72a74a95a60a95a53a71a73a80a17a59a91a80a40a72a74a78a81a103a99a62a153a68a70a78a79a197a180a62a153a59a149a103a87a57a60a103a111a103a87a57a60a89a75a57a65a103a87a57a60a69a17a90a122a71a73a80a17a59a149a103a17a68a41a98a99a57a60a72a41a69a107a57a92a69a75a71a73a68a181a64a79a101a17a112a17a57a60a69a70a71a73a59a19a78a73a89a54a72a74a95a65a64a19a192a75a62a153a59a91a78a79a59a48a72a74a95a60a57a92a183a23a59 a71a79a80a40a72a54a71a135a62a153a59a99a61a19a72a41a69a17a69a17a68a74a71a135a59a19a89a70a59a19a69a36a57a60a69a75a71a79a59a23a90a41a78a81a72a54a71a79a59a180a72a5a77a40a72a74a78a81a72a74a112a43a68a41a95a65a72a107a59a9a209a17a72a41a61a9a71a79a95a60a97a102a68a41a69a179a71a79a80a17a59a99a103a87a68a41a98a107a72a41a57a92a69a162a118a122a191a117a80a53a59a19a78a73a59a181a98a181a101a53a64a113a71 a112a43a59a185a64a66a68a70a98a99a59a9a71a79a80a53a57a92a69a17a90a85a71a79a80a53a72a74a71a91a62a153a59a181a61a19a72a74a69a36a103a17a68a107a71a79a68a85a57a60a98a99a77a17a78a73a68a54a89a41a59a15a71a79a80a17a57a65a64a15a64a79a61a81a80a17a59a23a98a99a59a41a118a117a186a50a59a185a90a41a68a85a112a53a72a70a61a81a197a111a72a41a69a53a103a110a95a92a68a25a68a70a197a111a72a54a71 a71a79a80a53a59a185a90a41a59a23a69a17a59a19a78a81a72a74a95a162a96a114a68a41a78a73a98 a68a41a96a30a71a73a80a17a59a181a116a75a101a53a72a70a103a87a78a81a72a54a71a79a101a53a78a79a59a122a72a41a77a17a77a17a78a73a68a27a209a25a57a60a98a107a72a54a71a73a57a92a68a70a69a12a64a73a61a81a80a17a59a19a98a99a59a70a118 a109 a95a60a95a7a62a128a59a185a72a74a78a73a59a122a103a87a68a70a57a92a69a17a90a5a57a60a64 a72a74a77a53a77a17a78a79a68a27a209a87a57a60a98a107a72a54a71a79a57a60a69a17a90a99a72a74a69a110a57a60a69a70a71a73a59a19a90a70a78a73a72a41a95a162a112a25a97a111a72a99a62a128a59a19a57a60a90a41a80a75a71a79a59a48a103a12a64a79a101a17a98 a68a74a96a46a96a114a101a17a69a53a61a9a71a79a57a60a68a41a69a102a59a19a89a54a72a41a95a92a101a53a72a74a71a79a57a60a68a41a69a53a64a63a72a70a64a76a64a79a80a17a68a54a62a63a69 a57a60a69 a71a73a80a17a57a60a64a5a64a66a95a60a57a60a103a17a59a41a118 a115 a68a174a96a32a72a74a78a107a62a128a59a12a80a40a72a27a89a41a59a12a112a40a59a23a59a19a69a157a61a81a80a17a68a25a68a70a64a79a57a92a69a17a90a50a71a79a80a53a59a23a64a79a59a12a62a128a59a19a57a60a90a41a80a75a71a81a64a107a72a41a64a181a71a79a80a17a59a36a64a79a101a17a112a17a57a60a69a75a71a79a59a23a78a79a89a54a72a74a95 a95a60a59a19a69a17a90a41a71a79a80a53a64a23a118a46a186a50a59a67a80a53a72a27a89a70a59a15a72a41a95a60a64a79a68a180a112a40a59a23a59a19a69a110a61a81a80a17a68a25a68a70a64a79a57a92a69a17a90a99a72a74a95a60a95a52a71a79a80a53a59a67a59a19a89a54a72a74a95a60a101a53a72a54a71a73a57a92a68a70a69a5a77a43a68a41a57a60a69a75a71a73a64a23a118a58a191a117a80a17a59a135a62a153a59a23a57a92a90a70a80a70a71a81a64a117a72a74a78a73a59 a16a113a101a53a64a66a71a122a64a66a68a70a98a180a59a180a69a17a68a70a78a79a98a107a72a74a95a60a57a60a183a19a57a60a69a17a90a5a96a32a72a70a61a177a71a73a68a41a78a81a64a91a62a63a80a17a57a65a61a81a80a174a59a19a69a53a64a79a101a17a78a73a59a185a71a73a80a53a72a54a71a135a71a79a80a53a59a99a72a41a77a17a77a17a78a73a68a27a209a87a57a92a98a107a72a54a71a73a57a92a68a70a69a36a57a65a64a149a59a19a209a17a72a41a61a177a71a67a57a92a96 a18 a21a114a207 a23a26a25a97a47a128a72a74a69a53a103a122a71a79a80a17a59a128a59a23a116a75a101a53a72a41a95a92a57a92a71a113a97a135a68a74a96a53a72a41a78a79a59a48a72a41a64a162a68a41a96a87a71a79a78a81a72a74a77a43a59a19a183a23a68a41a57a65a103a17a64a7a72a41a69a53a103a67a78a73a59a23a61a9a71a73a72a74a69a53a90a41a95a60a59a23a64a47a71a79a80a53a72a74a71a7a62a128a59a128a103a87a57a65a64a73a61a9a101a53a64a73a64a66a59a48a103 a90a41a57a60a89a41a59a48a64a76a101a53a64a149a71a73a80a17a59a181a59a19a209a25a71a79a78a81a72a85a77a43a68a41a95a60a97a75a69a53a68a41a98a99a57a60a72a41a95a30a72a70a61a19a61a19a101a17a78a73a72a70a61a9a97a111a68a74a96a94a112a43a59a19a57a60a69a17a90a110a72a74a112a17a95a60a59a122a71a73a68a111a68a70a112a87a71a73a72a41a57a92a69a102a71a73a80a17a59a99a72a74a78a73a59a23a72a107a101a17a69a40a103a87a59a19a78 a72a102a64a66a71a79a78a81a72a74a57a60a69a17a90a102a95a92a57a60a69a17a59a85a59a19a209a17a72a41a61a177a71a73a95a92a97a70a118 a115 a68a53a192a46a62a63a80a53a72a54a71a122a62a128a68a41a101a53a95a60a103a125a80a53a72a41a77a17a77a43a59a19a69a45a57a92a96a117a62a128a59a85a62a128a59a19a78a73a59a99a71a79a68a179a61a81a80a53a68a75a68a75a64a66a59a99a112a43a68a74a71a73a80a45a71a79a80a17a59 a57a60a69a70a71a73a59a19a90a70a78a73a72a74a71a79a57a60a68a41a69a122a77a43a68a41a57a60a69a75a71a73a64a143a72a74a69a40a103a122a71a73a80a17a59a117a62a153a59a23a57a92a90a70a80a70a71a81a64a30a57a92a69a75a71a79a59a23a95a92a95a60a57a60a90a41a59a19a69a75a71a73a95a92a97 a27a18a24a17a68a70a78a46a72a41a69a185a198a52a195a38a77a40a68a70a57a92a69a75a71a46a96a114a68a41a78a73a98a185a101a53a95a60a72a149a62a153a59a128a80a53a72a27a89a70a59 a32 a25 a198a143a188a34a62a153a59a23a57a92a90a70a80a70a71a81a64a67a72a74a69a53a103 a25 a198a143a188a162a59a19a89a54a72a74a95a60a101a53a72a74a71a79a57a60a68a41a69a174a77a40a68a70a57a92a69a75a71a73a64a15a71a73a68a102a61a81a80a17a68a25a68a70a64a79a59a41a118a107a191a117a80a53a72a54a71a67a90a70a57a92a89a70a59a23a64a15a101a53a64 a25 a23 a198a143a188a162a103a87a59a23a90a41a78a73a59a19a59a48a64a149a68a41a96 a96a114a78a73a59a19a59a23a103a17a68a41a98a110a118 a33a76a59a19a69a40a61a9a59a107a62a128a59a180a98a181a101a53a64a66a71a122a112a40a59a85a72a41a112a17a95a92a59a180a71a73a68a12a57a60a69a75a71a79a59a23a90a41a78a81a72a54a71a79a59a107a72a12a77a43a68a41a95a60a97a75a69a53a68a41a98a99a57a60a72a41a95a143a68a74a96a117a103a87a59a19a90a70a78a79a59a23a59a99a72a74a71a67a98a180a68a75a64a113a71 a25a21 a23 a198a62a63a10a47a23a9a188a63a118 a191a117a80a53a57a60a64a107a64a79a57a92a98a99a77a17a95a60a59a110a57a60a103a17a59a23a72a50a90a41a57a60a89a41a59a23a64a180a78a73a57a60a64a79a59a111a71a79a68a174a71a79a80a17a59a102a64a79a68a174a61a23a72a74a95a60a95a92a59a48a103 a22 a4a135a72a74a101a53a64a73a64a99a116a70a101a40a72a41a103a87a78a81a72a54a71a73a101a17a78a79a59 a22 a64a73a61a81a80a17a59a19a98a99a59a41a118 a14a47a161a14a47a161a65a160 a33a149a37a40a29a51a28a52a20 a31a1a0 a22a75a29a51a49a40a26a52a20a3a2a47a22a25a55a32a22a25a24a74a20a48a29a32a37a43a28 a164 a39a54a29a114a20a48a22a70a39a54a29a51a14 a165a43a166a75a167a169a168a47a170a5a4 a6a76a59a23a64a79a101a17a95a206a71a91a64a79a80a17a68a70a101a17a95a60a103a12a112a43a59a67a59a9a209a17a72a41a61a9a71a63a57a206a96 a18 a21a114a207 a23a128a57a60a64a63a72a99a77a43a68a41a95a60a97a25a69a17a68a41a98a99a57a65a72a74a95 a18 a21a114a207 a23 a25a8a7 a30 a44a9a7 a29 a207a57a44a9a7 a42 a207 a42 a44a11a10a12a10a12a10a55a44a13a7 a79 a207 a79 a25a13a14 a79 a21a32a207a24a23 a115 a59a23a95a92a59a48a61a177a71a76a207a40a202a79a188 a64a63a72a74a69a53a103 a20 a202a79a188 a64a63a64a79a101a53a61a81a80a5a71a73a80a53a72a54a71 a27 a29 a30 a14 a79 a21a32a207a24a23a84a38a41a207a86a25 a56 a38 a202a58 a29 a20 a202a15a14 a79 a21a114a207a40a202a77a23 a96a114a68a41a78 a109 a99a17a16a194a77a43a68a41a95a60a97a25a69a17a68a41a98a99a57a65a72a74a95a34a101a17a77a87a71a73a68a43a101a32a72a41a69a53a103a5a57a60a69a53a61a19a95a92a101a53a103a17a57a92a69a17a90a41a102a19a18a114a204a32a205a99a68a70a78a73a103a17a59a19a78 a0 a29a114a20a48a26 a23 a198 a86a30a22a75a49a40a39a54a22a25a22a75a16a99a37a43a35a76a35a200a39a54a22a75a22a25a86a30a37a40a42a11a20a21a18a89a25 a23 a198a72a63 a47 a172a174a173a34a175a177a176a23a22 a187a162a59a19a71a24a14 a79 a21a114a207a24a23a149a103a87a59a23a69a17a68a74a71a73a59a181a72a107a77a43a68a41a95a60a97a25a69a17a68a41a98a99a57a65a72a74a95a30a68a74a96a58a103a87a59a19a90a70a78a79a59a23a59a25a18a58a57a60a69a102a71a73a80a17a59a185a89a54a72a74a78a73a57a65a72a74a112a17a95a60a59a122a207 a101a26a7 a79a28a27 a25 a45a102a177a118a91a186a50a59a180a62a128a72a41a69a75a71 a71a79a68a107a64a79a59a19a95a60a59a23a61a9a71a63a71a79a80a17a59a67a62a128a59a19a57a60a90a41a80a75a71a73a64a117a72a41a69a53a103a12a57a92a69a75a71a79a59a23a90a41a78a81a72a54a71a73a57a92a68a70a69a5a77a43a68a41a57a60a69a70a71a81a64a63a64a66a101a40a61a81a80a5a71a73a80a53a72a54a71a63a71a73a80a17a59a135a96a114a68a41a78a73a98a185a101a17a95a65a72 a27 a29 a30 a18 a21a114a207 a23a39a38a70a207 a25 a56 a38 a202a58 a29 a20 a202 a14 a79 a21a114a207 a202a23 a57a65a64a135a59a9a209a17a72a41a61a9a71a15a96a114a68a70a78a122a72a74a95a60a95a58a77a43a68a41a95a60a97a75a69a53a68a41a98a99a57a60a72a41a95a60a64a15a68a41a96a117a103a87a59a23a90a41a78a73a59a19a59a180a101a17a77a87a71a73a68a61a101a32a72a74a69a40a103a50a57a60a69a53a61a19a95a92a101a53a103a17a57a92a69a17a90a41a102a29a18a113a118 a19 a112a25a89a25a57a92a68a70a101a53a64a66a95a60a97a41a192a47a62a63a57a206a71a73a80 a23 a198a174a103a87a59a19a90a70a78a79a59a23a59a23a64a128a68a74a96a7a96a114a78a79a59a23a59a23a103a87a68a70a98a110a192a87a71a79a80a17a59a67a112a43a59a23a64a66a71a76a62a128a59a135a61a23a72a74a69a110a103a87a68a99a57a60a64a30a18 a25 a23 a198a21a63 a47a41a118 a14a47a161a14a47a161a14 a0 a26a52a18 a20a23a26a30a22a102a13a10a13a162a14a17a24a41a20a23a28a30a22a25a16a23a16 a164 a39a54a29a114a20a48a22a70a39a54a29a51a14 a6 a165a43a166a75a167a169a168a47a170a21a171a32a31 a88 a68a70a69a53a64a79a57a60a103a87a59a23a78a117a71a79a80a17a59a122a191a30a72a27a97a25a95a92a68a70a78a63a64a66a59a23a78a79a57a60a59a23a64a153a96a114a68a41a78 a18 a21a32a207a24a23 a18 a21a114a207a24a23 a25 a18 a21a18a45a23a25a44a34a33 a18 a21a18a45a23 a33 a207 a207 a44a35a10a36a10a12a10a58a44 a47 a18 a48 a33 a79 a18 a21a72a45a23 a33 a207 a79 a207 a79 a44a9a37 a79 a91 a29 a37 a79 a91 a29 a57a65a64a128a71a79a80a17a59a67a39a54a22a25a42a146a14a17a29a51a28a7a86a30a22a75a39 a37 a79 a91 a29 a25 a47 a21a38a18a11a44 a47a23 a48 a33 a79 a91 a29 a18 a21a12a39a207a24a23 a33 a207 a79 a91 a29 a207 a79 a91 a29 a62a63a80a17a59a23a78a79a59a40a39a207a86a42 a44a45 a33a66a207 a49 a172a174a173a34a175a177a176a42a41 a19 a96a143a72a74a95a60a95a52a96a114a101a17a69a53a61a177a71a73a57a92a68a70a69a53a64a23a192a17a62a63a80a25a97a111a72a74a78a73a59a15a62a153a59a67a57a60a69a75a71a79a59a19a78a73a59a23a64a66a71a79a59a48a103a5a57a60a69a110a57a92a69a75a71a73a59a19a90a41a78a81a72a54a71a73a57a92a69a53a90a181a77a43a68a41a95a60a97a25a69a17a68a41a98a99a57a65a72a74a95a65a64a28a27a102a191a117a80a53a59a67a78a79a59a48a72a41a64a79a68a41a69 a61a9a68a70a98a99a59a23a64a67a96a114a78a73a68a41a98 a71a79a80a53a59a111a64a66a71a79a78a73a101a53a61a177a71a73a101a17a78a79a59a5a68a74a96a76a191a7a72a27a97a25a95a60a68a41a78a48a188 a64a185a64a66a59a23a78a79a57a60a59a23a64a67a59a19a209a87a77a53a72a74a69a53a64a79a57a60a68a41a69a162a118a36a191a117a80a17a59a111a191a7a72a27a97a25a95a60a68a41a78a122a59a19a209a25a77a40a72a74a69a53a64a79a57a92a68a70a69 a43 a68a74a96a58a72a99a96a114a101a53a69a53a61a177a71a73a57a92a68a70a69a36a57a92a69a36a72a85a95a92a68a87a61a19a72a41a95a47a69a17a59a23a57a92a90a70a80a25a112a40a68a70a78a79a80a17a68a25a68a87a103a12a68a74a96a94a72a99a77a43a68a41a57a60a69a75a71a22a101a114a80a53a59a19a78a73a59a67a71a79a80a17a57a65a64a149a77a40a68a70a57a92a69a75a71a91a57a65a64a91a61a81a80a17a68a75a64a66a59a23a69a36a72a41a64a1a0 a62a63a57a92a71a79a80a17a68a70a101a87a71a76a95a92a68a75a64a79a64a128a68a41a96a30a90a70a59a19a69a17a59a23a78a73a72a41a95a92a57a92a71a113a97a46a102a58a57a65a64a117a69a17a68a74a71a73a80a17a57a60a69a17a90a107a112a17a101a87a71a91a72a180a77a43a68a54a62a153a59a23a78a117a64a79a59a19a78a73a57a60a59a23a64a117a59a9a209a87a77a53a72a41a69a53a64a79a57a92a68a70a69 a15a30a191a117a80a53a59a67a80a17a57a92a90a70a80a17a59a19a78 a71a79a80a53a59a122a68a41a78a81a103a87a59a19a78a63a68a41a96a30a77a43a68a41a95a60a97a25a69a17a68a41a98a99a57a65a72a74a95a65a64a117a71a79a80a40a72a54a71a76a68a41a101a53a78a91a64a79a61a81a80a53a59a19a98a99a59a122a61a19a72a41a69a12a57a60a69a75a71a79a59a23a90a41a78a81a72a54a71a79a59a135a71a73a80a17a59a67a80a17a57a60a90a41a80a17a59a23a78a63a71a79a80a53a59a122a68a41a78a81a103a87a59a19a78a63a68a41a96 a71a79a80a53a59a111a78a73a59a19a98a107a72a74a57a60a69a53a103a87a59a23a78a180a71a79a59a19a78a73a98 a57a60a69 a71a73a80a17a59a111a59a19a209a87a77a53a72a74a69a53a64a79a57a60a68a41a69a162a118 a191a117a80a53a59a111a57a60a69a75a71a79a59a19a90a70a78a73a72a41a95a153a68a41a96a76a71a79a80a53a59a12a78a73a59a19a98a107a72a41a57a92a69a17a59a23a78a181a68a54a89a70a59a19a78a181a71a79a80a17a59 a103a87a68a70a98a99a72a41a57a92a69a12a57a65a64a117a77a17a78a79a59a48a61a9a57a65a64a66a59a23a95a92a97a85a71a79a80a17a59a67a59a23a78a79a78a73a68a41a78a128a57a92a69a110a69a25a101a17a98a99a59a19a78a73a57a65a61a19a72a74a95a34a57a60a69a75a71a79a59a23a90a41a78a81a72a54a71a79a57a60a68a41a69a47a118 a14a47a161a14a47a161a3a2 a13a15a16a19a20a48a29a32a42a146a14a17a20a23a29a32a28a30a49a50a20a23a26a30a22a102a13a91a39a27a39a54a37a53a39 a165a43a166a75a167a169a168a47a170a21a171a17a171 a4a91a64a79a57a60a69a17a90a180a71a79a80a17a59a122a191a30a72a27a97a25a95a92a68a70a78a117a64a79a59a19a78a73a57a60a59a23a64a128a78a79a59a48a64a66a101a53a95a206a71a81a64a63a72a74a69a53a103a111a71a79a80a53a59a67a59a9a209a17a72a41a61a9a71a79a69a17a59a48a64a79a64a117a61a9a78a73a57a92a71a79a59a19a78a73a57a65a72 a5 a29a30 a18 a21a32a207a24a23a39a38a70a207a21a63 a54a57a56 a202a58 a29 a20 a202 a18 a21a32a207a40a202a52a23 a25 a47 a21a38a18a11a44 a47a23 a48 a27a16a29 a30 a33 a79 a91 a29 a18 a21a12a39a207 a21a114a207a24a23a53a23 a33 a207 a79 a91 a29 a207 a79 a91 a29 a38a41a207 a67 a68a69 a70 a6 a74a8a7a8a78 a202 a56 a83a74a8a9 a109 a64a73a64a66a101a53a98a180a57a60a69a17a90a5a103a17a59a19a78a73a57a92a89a54a72a54a195 a71a79a57a60a89a41a59a48a64a117a68a74a96 a18 a21a114a207a24a23a76a72a74a78a73a59a149a112a43a68a41a101a17a69a40a103a87a59a23a103a12a68a41a69a103a44a45a46a33a48a47a6a49 a10 a10 a10 a10 a10 a27 a29 a30a19a18 a21a114a207a24a23a84a38a41a207a72a63 a56 a38 a202a58 a29 a20 a202 a18 a21a114a207 a202 a23 a10 a10 a10 a10 a10a12a11 a13 a21a26a18a11a44a57a47a23 a48 a88 a68a41a69a25a89a70a59a19a78a73a90a41a59a19a69a40a61a9a59a91a57a65a64 a3a40a22a75a39a27a18a125a96a32a72a41a64a66a71a16a15a4a15 a172a174a173a34a175a177a176a179a178a15a14 a109 a64a73a64a79a101a17a98a99a59a122a71a79a80a40a72a54a71a135a68a41a101a17a78a135a64a79a61a81a80a53a59a19a98a99a59a185a57a65a64a91a59a19a209a17a72a41a61a177a71a135a101a17a77a87a71a73a68a111a72a5a77a40a68a70a95a92a97a25a69a17a68a70a98a180a57a65a72a74a95a30a68a41a78a81a103a87a59a23a78 a25a18a113a188a60a118a135a191a117a80a40a72a54a71a149a98a99a59a48a72a74a69a53a64a149a62a153a59 a61a19a72a41a69a85a57a60a69a70a71a73a59a19a90a70a78a73a72a74a71a79a59a76a71a73a80a17a59a149a208a53a78a73a64a66a71a2a21a26a18 a44a16a47a23a46a71a73a59a19a78a73a98a99a64a128a57a60a69a85a71a73a80a17a57a60a64a117a191a30a72a27a97a25a95a92a68a70a78a128a64a79a59a19a78a73a57a92a59a48a64a58a59a19a209a25a77a40a72a74a69a53a64a79a57a92a68a70a69a85a59a9a209a17a72a70a61a177a71a79a95a60a97a41a118a143a191a117a80a17a59 a59a19a78a73a78a73a68a41a78a128a57a92a69a110a69a25a101a17a98a99a59a19a78a73a57a60a61a23a72a74a95a34a57a60a69a75a71a79a59a19a90a70a78a73a72a74a71a79a57a60a68a41a69 a71 a25 a27 a29 a30 a18 a21a114a207 a23a39a38a70a207 a63 a56 a38 a202a58 a29 a20 a202 a18 a21a114a207 a202a23a26a25 a47 a21a26a18a11a44a57a47a23 a48 a27 a29 a30 a33 a79 a91 a29 a18 a21a36a39a207a34a21a114a207 a23a84a23 a33 a207 a79 a91 a29 a207 a79 a91 a29 a38a70a207a89a94 a14a47a161a14a47a161a17a16 a11a194a22a25a22a75a20a23a29a51a28a7a49a174a20a48a26a7a22a102a13 a13a162a14a17a24a74a20a48a28a7a22a25a16a48a16 a164 a39a27a29a32a20a23a22a75a39a54a29a32a14 a165a43a166a75a167a169a168a47a170a21a171 a15 a93a143a209a17a72a41a61a9a71a79a69a17a59a48a64a79a64a63a61a19a68a41a69a53a103a17a57a206a71a73a57a92a68a70a69a111a78a73a59a23a116a75a101a17a57a60a78a79a59a48a64 a27 a29 a30 a14 a79 a21a114a207 a23a39a38a70a207 a25 a27 a29 a30 a21a38a7 a30 a44a13a7 a29 a207a57a44a13a7 a42 a207 a42 a44a11a10a12a10a36a10a55a44a9a7 a79 a207 a79a23a39a38a70a207 a25 a56 a38 a202a58 a29 a20 a202 a14 a79 a21a32a207 a202 a23 a96a114a68a41a78a76a72a41a69a25a97a5a64a66a59a19a71a63a68a74a96 a18a11a44 a47a67a61a9a68a25a59a19a18a107a61a19a57a92a59a23a69a75a71a73a64a17a7 a30 a33 a7 a29 a33a88a94a48a94a88a94a88a33 a7 a79 a13a15a31a34a33a30a29 a3a17a14a17a55a51a22a75a28a52a20a48a55a114a18 a27 a29 a30 a7 a30 a38a41a207a57a44 a27 a29 a30 a7 a29 a207a96a38a70a207 a44 a27 a29 a30 a7 a42 a207 a42 a38a70a207a57a44a35a10a36a10a12a10a58a44 a27 a29 a30 a7 a79 a207 a79a38a70a207 a25 a56 a38 a202a58 a29 a20 a202 a14 a79 a21a32a207 a202 a23 a172a174a173a34a175a177a176a179a178a43a178 a191a117a80a17a57a65a64a63a64a66a95a60a57a65a103a87a59a67a69a17a59a19a59a48a103a17a64a117a95a60a57a206a71a79a71a79a95a60a59a122a61a9a95a65a72a74a78a73a57a206a208a40a61a23a72a54a71a73a57a92a68a70a69a162a118 a19 a101a17a78a63a59a9a209a17a72a41a61a9a71a79a69a17a59a48a64a79a64a117a61a19a78a79a57a92a71a79a59a23a78a79a57a60a68a41a69a12a57a60a64 a27 a29 a30 a14 a79 a21a114a207 a23a39a38a70a207 a25 a27 a29 a30 a21a38a7 a30 a44a13a7 a29 a207a57a44a13a7 a42 a207 a42 a44a11a10a12a10a36a10a55a44a9a7 a79 a207 a79a23a39a38a70a207 a25 a56 a38 a202a58 a29 a20 a202 a14 a79 a21a32a207 a202 a23 a20 a62a63a80a17a57a65a61a81a80a111a57a65a64a117a71a73a80a17a59a122a64a79a72a41a98a99a59a135a72a70a64 a7 a30 a27 a29 a30 a38a41a207 a44 a7 a29 a27 a29 a30 a207a96a38a70a207a53a44 a10a12a10a36a10 a44 a7 a79 a27 a29 a30 a207 a79a38a70a207 a25a8a7 a30 a56 a38 a202a58 a29 a20 a202 a44 a7 a29 a56 a38 a202a58 a29 a20 a202 a207 a202 a44 a10a12a10a36a10a32a44 a7 a79 a56 a38 a202a58 a29 a20 a202 a207 a202a79 a24a17a68a70a78a117a71a79a80a17a57a65a64a128a71a79a68a107a112a43a59a122a72a74a69a12a57a65a103a87a59a19a69a75a71a79a57a92a71a113a97a107a96a114a68a70a78a117a71a79a80a17a59 a21a38a18a11a44 a47a23a153a72a41a78a79a112a53a57a206a71a73a78a73a72a41a78a79a97a85a61a19a68a75a59a19a18a85a61a9a57a60a59a19a69a75a71a73a64a19a7 a202 a192a53a62a153a59a135a98a185a101a40a64a113a71a63a80a53a72a27a89a70a59 a71a79a80a53a59 a21a26a18a11a44a57a47a23a128a61a9a68a70a69a53a103a87a57a92a71a79a57a60a68a41a69a53a64 a56 a38 a202a58 a29 a20 a202a51a207a40a202 a2 a25a28a27 a29 a30 a207 a2 a38a70a207 a18a53a52a55a54 a9 a25a32a45a46a33a48a47a98a33a48a94a88a94a88a94a33 a18 a14a47a161a14a47a161a1a0 a11a194a22a25a22a75a20a23a29a51a28a7a49a174a20a48a26a7a22a102a13 a13a162a14a17a24a74a20a48a28a7a22a25a16a48a16 a164 a39a27a29a32a20a23a22a75a39a54a29a32a14 a165a43a166a75a167a169a168a47a170a21a171a98a19 a93a143a209a17a72a41a61a9a71a79a69a17a59a48a64a79a64a63a61a19a68a41a69a53a103a17a57a206a71a73a57a92a68a70a69a111a62a63a57a60a95a92a95a47a112a40a59a122a64a73a72a54a71a73a57a60a64a66a208a53a59a23a103a111a57a92a96a46a72a41a69a53a103a5a68a70a69a17a95a60a97a5a57a92a96 a27a16a29 a30 a38a70a207 a25 a56 a38 a202a58 a29 a20 a202 a10 a47 a27 a29 a30 a207a24a38a41a207 a25 a56 a38 a202a58 a29 a20 a202 a10a9a207 a202 a118 a118 a118 a27 a29 a30 a207 a79a38a41a207 a25 a56 a38 a202a58 a29 a20 a202 a10a9a207 a202a79 a14a47a161a14a47a161a3a2 a11a194a22a25a22a75a20a23a29a51a28a7a49a174a20a48a26a7a22a102a13 a13a162a14a17a24a74a20a48a28a7a22a25a16a48a16 a164 a39a27a29a32a20a23a22a75a39a54a29a32a14 a165a43a166a75a167a169a168a47a170a21a171a69a16 a6a63a59a23a68a41a78a73a90a70a72a74a69a53a57a92a183a23a57a92a69a17a90a181a59a9a209a17a72a41a61a9a71a79a69a17a59a48a64a79a64a117a59a48a116a75a101a53a72a54a71a73a57a92a68a70a69a53a64 SMA-HPC ?2000 MIT Normalized 1-D Problem General Quadrature Scheme Meeting the exactness criteria 1 12 2 1 12 0 111 1 0n lll l nn w xx x w xdxxx xw #$#$#$ %&%&%& %& %&?= %&%&%& %&’(’( %&’(" % & % & &&’&& % a15 a37a43a28a7a55a51a29a32a28a30a22a75a14a17a39a75a192a87a64a79a57a92a69a40a61a9a59a67a207a40a202a79a188 a64a63a72a74a69a53a103 a20 a202a66a188 a64a63a72a74a78a73a59a135a101a17a69a17a197a25a69a17a68a54a62a63a69a53a64 a172a174a173a34a175a177a176a179a178 a17 a99a76a68a54a62 a62a63a80a40a72a54a71a135a57a60a64a135a72a85a77a53a78a73a72a70a61a177a71a79a57a65a61a19a72a41a95a7a62a117a72a27a97a110a68a74a96a153a61a19a68a41a98a99a77a17a101a87a71a73a57a92a69a53a90a5a71a73a80a17a59a181a59a23a89a54a72a74a95a60a101a53a72a54a71a73a57a92a68a70a69a36a77a40a68a70a57a92a69a75a71a81a64a15a72a41a69a53a103a36a62a128a59a19a57a60a90a41a80a75a71a81a64a28a27 a186a174a59a76a61a19a72a74a69a180a62a63a78a73a57a206a71a73a59a117a71a79a80a17a59a63a59a19a209a87a72a70a61a177a71a73a69a17a59a23a64a73a64a46a61a9a78a73a57a206a71a73a59a19a78a73a57a60a72a15a57a60a69a70a71a73a68a67a72a149a98a107a72a74a71a79a78a73a57a206a209a122a96a114a78a73a68a41a98a110a118a46a191a117a80a17a59a76a64a79a97a25a64a66a71a79a59a23a98a194a68a74a96a52a59a23a116a75a101a53a72a74a71a79a57a60a68a41a69a53a64 a57a65a64a117a69a17a68a74a71a76a59a48a72a41a64a79a97a99a71a79a68a85a64a79a68a41a95a60a89a41a59a135a64a66a57a60a69a53a61a9a59a67a207 a202 a188 a64a63a72a74a69a40a103 a20 a202 a188 a64a63a72a41a78a79a59a135a101a17a69a53a197a75a69a53a68a54a62a63a69a53a64a19a118 a210 a0 a14a47a161a14a47a161a1a0 a164 a37a43a42a45a44a30a33a7a20a23a29a51a28a7a49a50a20a23a26a30a22a6a33a149a37a40a29a51a28a52a20a23a16a99a14a17a28a30a86a35a0 a22a25a29a32a49a43a26a52a20a23a16 a165a43a166a75a167a169a168a47a170a21a171 a7 a88 a68a41a101a53a95a60a103a111a101a53a64a79a59 a15 a22 a3a99a20a48a37a40a28a3a2a92a16 a11a194a22a75a20a23a26a30a37a162a86 a4a90a21a6a5a46a23a26a25a57a45a8a7a10a9a12a11 a22 a5a14a13 a24 a22 a5a14a13 a91 a29 a63a15a5a14a13 a24 a25a8a63a16a4 a22 a5a14a13 a24 a191a117a80a17a59a135a69a17a68a70a69a17a95a92a57a60a69a17a59a48a72a74a78a117a96a114a101a17a69a53a61a9a71a79a57a60a68a41a69a12a96a114a68a41a78a12a99a91a59a19a62a117a71a73a68a41a69a110a57a60a64a128a71a73a80a17a59a19a69 SMA-HPC ?2000 MIT Normalized 1-D Problem General Quadrature Scheme Computing the points and weights 1 12 2 1 12 0 111 1 0n lll l nn w xx x wwF x xdxxx xw #$#$#$ %&%&%& %&)* %&=?=+, %&-.%&%& %&’(’( %&’(" % & % & &&’&& % a172a174a173a34a175a177a176a179a178 a85 a99a76a59a23a62a117a71a79a68a70a69a162a188 a64a30a98a99a59a9a71a73a80a17a68a87a103a122a57a65a64a46a72a74a69a122a57a92a71a79a59a23a78a73a72a74a71a79a57a60a89a41a59a58a71a73a59a23a61a81a80a17a69a53a57a60a116a75a101a17a59a153a96a114a68a41a78a7a208a53a69a53a103a87a57a60a69a17a90a67a72a76a89a54a72a41a95a92a101a17a59a17a5a67a64a66a101a40a61a81a80a122a71a79a80a40a72a54a71a18a4a90a21a19a5a46a23a26a25 a45a17a118a46a191a117a80a17a59a76a98a99a59a19a71a79a80a17a68a87a103a107a57a65a64a58a112a53a72a70a64a66a59a48a103a180a68a41a69a107a95a60a57a92a69a17a59a48a72a74a78a73a57a92a183a23a57a92a69a53a90a135a71a79a80a17a59a76a77a53a78a79a68a70a112a17a95a92a59a23a98 a72a41a112a40a68a70a101a87a71a94a72a67a90a70a101a17a59a23a64a73a64a58a72a54a71a3a5a43a192a25a72a74a69a40a103a180a71a79a80a17a59a23a69 a101a17a77a52a103a17a72a54a71a73a57a92a69a53a90a102a71a79a80a17a59a111a89a54a72a74a95a60a101a17a59a85a68a74a96a20a5a125a112a25a97a174a64a79a68a41a95a60a89a25a57a92a69a53a90a110a71a79a80a17a59a85a95a60a57a60a69a17a59a23a72a41a78a79a57a60a183a19a59a48a103a125a77a17a78a79a68a70a112a17a95a60a59a19a98a110a118a102a104a105a69a45a77a53a72a74a78a79a71a79a57a65a61a9a101a53a95a60a72a41a78a23a192a7a71a79a80a17a59 a57a92a71a79a59a19a78a81a72a54a71a73a59a21a5 a13 a91 a29 a57a60a64a91a103a87a59a9a71a73a59a19a78a73a98a180a57a60a69a17a59a48a103a85a96a114a78a73a68a41a98a22a5 a13 a112a25a97a5a64a79a68a41a95a60a89a75a57a60a69a17a90a181a71a79a80a17a59a67a95a60a57a60a69a17a59a23a72a41a78a76a64a66a97a87a64a66a71a79a59a19a98a132a68a74a96a46a59a23a116a75a101a53a72a54a71a73a57a92a68a70a69a53a64 a4a90a21a6a5 a13 a23 a44a23a9 a11 a21a6a5 a13 a23a7a22a24a5 a13 a91 a29 a63a15a5 a13 a24 a25a57a45 a62a63a80a17a59a23a78a79a59a25a9a26a11 a21a19a5 a13 a23a107a57a65a64a107a71a73a80a17a59a28a27a75a72a41a61a19a68a41a112a17a57a65a72a74a69a10a101a114a98a181a101a17a95a206a71a73a57a60a103a17a57a92a98a99a59a19a69a40a64a66a57a60a68a41a69a53a72a41a95a15a103a17a59a19a78a73a57a92a89a54a72a54a71a73a57a92a89a70a59 a102a180a68a41a96a67a71a79a80a17a59a36a69a53a68a41a69a17a95a60a57a92a69a53a59a23a72a74a78 a96a114a101a17a69a53a61a9a71a79a57a60a68a41a69a29a4a90a21a19a5a46a23a9a118a134a191a117a80a17a59a102a57a206a71a73a59a19a78a81a72a54a71a73a57a92a68a70a69a21a57a65a64a107a61a9a68a70a69a70a71a73a57a92a69a25a101a17a59a48a103 a101a53a69a70a71a73a57a92a95a76a71a79a80a17a59a102a101a17a77a52a103a17a72a74a71a79a59a23a103a30a5a21a57a60a64a107a64a79a101 a18a85a61a9a57a60a59a19a69a75a71a79a95a60a97 a61a9a95a60a68a70a64a79a59a135a71a79a68a85a71a79a80a17a59a122a59a19a209a17a72a41a61a177a71a15a64a66a68a70a95a92a101a87a71a73a57a92a68a70a69a162a192a53a72a85a61a19a78a79a57a92a71a79a59a23a78a79a57a60a68a41a69a12a71a73a80a53a72a54a71a149a61a23a72a74a69a102a112a43a59a185a103a87a57a3a18a85a61a9a101a17a95a92a71a91a71a73a68a85a89a41a59a23a78a79a57a92a96a114a97a41a118a75a99a76a59a23a62a117a71a79a68a41a69a47a188 a64 a98a99a59a9a71a73a80a17a68a87a103a181a103a87a68a25a59a48a64a7a69a53a68a74a71a58a72a74a95a60a62a128a72a27a97a87a64a7a61a9a68a41a69a25a89a70a59a19a78a73a90a41a59a41a192a54a72a149a77a17a80a17a59a23a69a17a68a41a98a99a59a19a69a53a68a41a69a122a71a73a80a53a72a54a71a143a57a65a64a46a98a99a68a41a78a73a59a153a95a60a57a92a197a70a59a19a95a60a97a67a62a63a80a17a59a19a69a31a9 a11 a21a6a5 a23 a57a65a64a47a69a17a59a48a72a74a78a73a95a92a97a135a64a79a57a92a69a17a90a70a101a17a95a65a72a74a78a48a118 a24a17a68a70a78a47a98a99a68a41a78a73a59a153a72a41a112a40a68a70a101a87a71 a99a91a59a19a62a117a71a79a68a70a69a162a188 a64a47a98a99a59a9a71a73a80a17a68a87a103a34a192a74a64a79a59a19a59a58a71a73a80a17a59 a35 a118 a189a70a189 a35a33a32a25a210 a35 a118 a20a31a32 a0 a32 a32 a118 a0 a20a31a35 a61a9a68a70a101a17a78a81a64a66a59a15a69a17a68a74a71a73a59a23a64a23a118 a14a47a161a14a47a161a1a34 a164 a37a43a42a45a44a30a33a7a20a23a29a51a28a7a49a50a20a23a26a30a22a6a33a149a37a40a29a51a28a52a20a23a16a99a14a17a28a30a86a35a0 a22a25a29a32a49a43a26a52a20a23a16 a165a43a166a75a167a169a168a47a170a21a171 a37 a99a76a59a23a62a117a71a79a68a70a69a12a100a102a59a19a71a79a80a17a68a87a103a35a27a75a72a41a61a9a68a70a112a17a57a65a72a74a69a5a78a73a59a19a89a70a59a23a72a41a95a60a64a128a77a17a78a73a68a41a112a17a95a60a59a19a98 Normalized 1-D Problem General Quadrature Scheme Computing the points and weights 12 1 2 11 12 11 1 11 1 1 1 1 nn F ll l l l xx x w w wwJ x xx xlwx lwx?? #$ %&)* =+, -. ’( %% &&’& & &’ & %%% Columns become linearly dependent for high order 2n 00 0 a172a174a173a34a175a177a176a179a178a49a26 a187a162a68a25a68a70a197a75a57a60a69a17a90a102a72a54a71a67a71a73a80a17a59a36a27a75a72a41a61a19a68a41a112a17a57a65a72a74a69a50a68a74a96a117a71a79a80a17a59a107a77a53a78a79a68a70a112a17a95a92a59a23a98a110a192a7a62a128a59a107a78a79a59a48a72a74a95a60a57a92a183a23a59a180a71a79a80a53a72a74a71a67a71a79a80a17a59a107a208a40a78a73a64a66a71 a25 a198a143a188a7a61a9a68a41a95a60a101a17a98a99a69a53a64 a112a43a59a23a61a9a68a70a98a99a59a143a57a92a69a53a61a19a78a79a59a48a72a41a64a79a57a92a69a53a90a41a95a60a97a76a95a92a57a60a69a17a59a48a72a74a78a73a95a92a97a15a103a87a59a19a77a43a59a19a69a53a103a17a59a19a69a75a71a162a96a114a68a41a78a47a95a60a72a41a78a79a90a70a59 a18a66a118a47a191a117a80a17a57a65a64a162a57a65a64a34a112a43a68a41a101a17a69a53a103a15a71a73a68a76a80a53a72a41a77a17a77a40a59a23a69a67a64a66a57a60a69a53a61a9a59 a210a41a210 a62a128a59a185a72a74a78a73a59a67a95a92a68a25a68a70a197a75a57a60a69a17a90a107a57a60a69a75a71a79a68a107a71a73a80a17a59a185a64a79a77a53a72a41a61a19a59 a56 a14 a7a75a198a1a0 a47a33a66a207a89a33a48a94a88a94a48a94a33a66a207 a79a3a2a70a118a63a191a117a80a17a57a65a64a76a112a53a72a70a64a66a57a65a64a91a72a41a95a92a62a117a72a27a97a87a64a117a112a40a59a48a61a9a68a41a98a99a59a48a64a76a57a60a95a92a95 a61a9a68a70a69a53a103a87a57a92a71a79a57a60a68a41a69a17a59a48a103a107a62a63a57a92a71a79a80a85a57a60a69a53a61a9a78a73a59a23a72a70a64a66a57a60a69a17a90a25a18a66a118a30a191a117a80a53a59a149a64a79a68a41a95a60a101a87a71a79a57a60a68a41a69a85a57a65a64a94a71a79a68a185a68a70a112a87a71a73a72a41a57a92a69a85a72a185a77a40a68a70a95a92a97a25a69a17a68a70a98a99a57a60a72a41a95a40a112a53a72a70a64a66a57a65a64a143a71a73a80a53a72a54a71 a57a65a64a30a22a177a69a17a68a70a78a79a98a107a72a74a95a60a57a60a183a19a59a23a103 a22a76a57a60a69a12a64a79a68a41a98a99a59a67a64a79a59a19a69a53a64a79a59a122a64a66a68a181a71a79a80a40a72a54a71a76a57a92a71a76a57a60a64a63a77a53a78a79a68a70a77a40a59a23a78a79a95a60a97a85a61a9a68a70a69a53a103a87a57a92a71a79a57a60a68a41a69a17a59a48a103a34a118 a14a47a161a14a47a161a32a163a5a4 a6 a16a48a22a43a15a12a29a8a7a94a22a75a39a54a22a25a28a52a20 a33a149a37a40a55a32a18a162a28a7a37a43a42a146a29a51a14a17a55a32a16 a165a43a166a75a167a169a168a47a170a21a171 a16 a93a143a209a17a72a41a61a9a71a79a69a17a59a48a64a79a64a63a61a19a78a79a57a92a71a79a59a23a78a79a57a65a72a181a62a63a57a92a95a60a95a162a112a43a59a122a64a79a72a74a71a79a57a65a64a113a208a53a59a48a103a5a57a92a96a143a72a74a69a53a103a12a68a41a69a53a95a92a97a85a57a92a96 Normalized 1-D Problem General Quadrature Scheme Use different polynomials () ( ) 1 00 10 n ii i cxdx wcx = =!" () ( ) 1 11 10 n ii i cxdx wcx = =!" () ( ) 1 10 n lili i cxdx wcx = =!" Each polynomial mustic Contain an termix Be linearly independent BUT a172a174a173a34a175a177a176a179a178 a11 a191a117a80a17a59a99a68a41a69a53a95a92a97a36a103a87a57a92a182a52a59a19a78a73a59a19a69a53a61a19a59a185a96a114a78a73a68a41a98 a71a73a80a17a59a99a77a17a78a73a59a19a89a25a57a92a68a70a101a53a64a135a64a66a59a19a71a15a68a41a96a128a61a19a78a79a57a92a71a79a59a23a78a79a57a65a72a5a57a65a64a149a71a79a80a53a72a74a71a135a71a79a80a17a59a48a64a66a59a99a77a43a68a41a95a60a97a75a69a53a68a41a98a99a57a60a72a41a95a60a64 a80a53a72a27a89a70a59a149a112a43a59a9a71a79a71a79a59a19a78a91a77a17a78a79a68a70a77a40a59a23a78a66a71a73a57a92a59a48a64a128a71a79a80a53a72a41a69a5a71a73a80a17a59a67a68a41a69a53a59a23a64a117a62a128a59a135a61a81a80a53a68a70a64a79a59a15a112a43a59a9a96a114a68a70a78a79a59a70a118 a14a47a161a14a47a161a32a163a43a163 a9 a39a27a20a23a26a30a37a43a49a40a37a43a28a7a14a53a55a17a33a149a37a40a55a32a18a162a28a7a37a43a42a146a29a32a14a53a55a32a16 a165a43a166a75a167a169a168a47a170a21a171 a18 a24a17a68a70a78a117a71a79a80a17a59a67a69a53a68a41a78a73a98a99a72a41a95a92a57a60a183a19a59a48a103a85a57a60a69a75a71a79a59a19a90a70a78a73a72a41a95a51a192a75a71a113a62a128a68a99a77a40a68a54a97a25a69a17a68a70a98a180a57a65a72a74a95a65a64a117a72a74a78a73a59a67a64a79a72a41a57a60a103a5a71a79a68a99a112a43a59a122a37a53a39a27a20a23a26a30a37a43a49a40a37a43a28a7a14a53a55a7a57a92a96 a27 a29 a30a11a10 a202 a21a114a207 a23 a10 a2 a21a32a207a24a23a84a38a41a207 a25a32a45 a18a53a52a55a54 a9 a27 a25a134a203 a191a117a80a17a59a110a72a74a112a43a68a54a89a41a59a111a57a92a69a75a71a79a59a23a90a41a78a81a72a74a95a153a57a60a64a180a68a74a96a169a71a73a59a19a69 a78a79a59a19a96a114a59a19a78a73a78a79a59a48a103a174a71a79a68a125a72a41a64a99a72a74a69a21a57a92a69a53a69a17a59a19a78a99a77a17a78a73a68a87a103a87a101a53a61a177a71a107a72a74a69a40a103a146a72a41a64a73a61a9a78a73a57a92a112a43a59a23a103a146a71a79a80a17a59 a69a17a68a41a71a73a72a54a71a73a57a92a68a70a69 a21 a10 a202a53a33 a10 a2 a23a26a25a28a27a16a29 a30 a10 a202a84a21a32a207a24a23 a10 a2a65a21a114a207 a23a39a38a70a207 a191a117a80a17a59a99a61a9a68a70a69a17a69a17a59a23a61a9a71a79a57a60a68a41a69a50a112a40a59a19a71a113a62a153a59a23a59a19a69a179a77a43a68a41a95a60a97a25a69a17a68a41a98a99a57a65a72a74a95a30a57a92a69a53a69a17a59a19a78a135a77a17a78a73a68a87a103a87a101a53a61a177a71a81a64a135a72a74a69a53a103a179a89a41a59a23a61a9a71a79a68a70a78a149a57a60a69a17a69a17a59a23a78a15a77a53a78a79a68a87a103a87a101a53a61a9a71a73a64 a61a19a72a41a69a111a112a43a59a122a64a79a59a19a59a19a69a12a112a25a97a5a64a73a72a74a98a99a77a17a95a60a57a92a69a17a90a40a118 a14a47a161a14a47a161a32a163a17a160 a13a10a13a34a44a46a55a32a37a43a29a114a20a48a29a32a28a30a49a179a20a48a26a7a22 a15a111a29a12a7a94a22a75a39a27a22a25a28a52a20a36a33a91a37a43a55a114a18a162a28a30a37a40a42 a29a32a14a17a55a51a16 a165a43a166a75a167a169a168a47a170a21a171 a4 a88 a68a41a69a40a64a66a57a65a103a87a59a19a78a117a78a73a59a19a62a63a78a73a57a206a71a73a57a92a69a53a90a181a71a73a80a17a59a67a59a9a209a17a72a70a61a177a71a79a69a53a59a23a64a73a64a128a61a19a78a79a57a92a71a79a59a23a78a79a57a65a72 Normalized 1-D Problem General Quadrature Scheme Exploiting the different polynomials () ( ) 1 00 10 n ii i cxdx wcx = =!" () ( ) 1 11 10 Low order terms n nii i cxdx wcx?? = !" !"""""#"""""$ () ( ) 1 10 n nini i cxdx wcx = =!" () ( ) 1 21 21 10 High Order Terms n nini i cxdx wcx?? = ! !"""""#"""""$ a210 a32 a6a63a59a23a61a23a72a74a95a60a95a52a71a79a80a53a72a74a71a30a18a84a21a1a0 a14 a52 a18a6a5 a56 a23a4a25 a23 a198a21a63a16a47a65a21a1a0 a52 a18 a10 a52a69a60a69a18 a18 a203 a10 a203 a60 a198a3a2 a56 a23 a172a174a173a34a175a177a176a179a178 a18 a186a174a59a17a16a113a101a40a64a113a71a76a61a19a72a41a95a92a95a43a71a73a80a17a59a15a208a53a78a73a64a66a71a2a21a114a198 a63 a47a23a153a61a9a68a41a69a40a103a87a57a206a71a73a57a92a68a70a69a53a64a117a72a41a64a153a71a79a80a53a59 a22a9a95a92a68a54a62a128a59a19a78a153a68a41a78a81a103a87a59a19a78a153a71a79a59a23a78a79a98a107a64 a22a91a72a41a69a53a103a5a71a79a80a17a59a15a95a60a72a70a64a113a71 a198a50a61a19a68a41a69a53a103a17a57a206a71a73a57a92a68a70a69a53a64a63a72a41a64a128a71a73a80a17a59 a22a177a80a17a57a60a90a41a80a53a59a19a78a117a68a41a78a81a103a87a59a19a78a117a71a73a59a19a78a73a98a107a64 a22a74a118 a14a47a161a14a47a161a32a163 a14 a13a10a13a34a44a46a55a32a37a43a29a114a20a48a29a32a28a30a49a179a20a48a26a7a22 a15a111a29a12a7a94a22a75a39a27a22a25a28a52a20a36a33a91a37a43a55a114a18a162a28a30a37a40a42 a29a32a14a17a55a51a16 a165a43a166a75a167a169a168a47a170a16a15 a31 a88 a72a74a69a12a62a63a78a73a57a206a71a73a59a149a71a73a80a17a59a67a80a17a57a60a90a41a80a17a59a23a78a63a68a41a78a81a103a87a59a19a78a128a71a73a59a19a78a73a98a99a64a91a103a87a57a206a182a52a59a19a78a73a59a19a69a75a71a73a95a92a97 a5 a29a30 a10 a56 a21a32a207a24a23a84a38a41a207a86a25 a54 a56 a202a58 a29 a20 a202 a10 a56 a21a32a207a40a202a52a23 a7 a5 a29a30 a10 a56 a21a114a207a24a23 a10 a30 a21a114a207 a23a39a38a70a207 a25 a54 a56 a202a58 a29 a20 a202 a10 a56 a21a114a207a40a202a77a23 a10 a30 a21a32a207a40a202a77a23 a118 a118 a118 a118 a118 a118 a5 a29a30 a10 a42 a56 a11 a29 a21a32a207a24a23a39a38a70a207 a25 a54 a56 a202a58 a29 a20 a202 a10 a42 a56 a11 a29 a21a114a207a43a202a23 a7 a5 a29a30 a10 a56 a21a114a207 a23 a10 a56 a11 a29 a21a114a207a24a23a84a38a41a207a86a25 a54a57a56 a202a58 a29 a20 a202 a10 a56 a21a114a207a43a202a23 a10 a56 a11 a29 a21a32a207a40a202a52a23 a191a117a80a17a59a135a77a17a78a73a68a87a103a87a101a53a61a177a71a81a64 a10 a56 a21a32a207a24a23 a10 a2a65a21a114a207 a23a117a72a74a78a73a59a15a95a60a57a60a69a17a59a23a72a41a78a79a95a60a97a85a57a60a69a53a103a87a59a19a77a43a59a19a69a40a103a87a59a19a69a75a71a16a15 a172a174a173a34a175a177a176a179a178 a30 a104a105a69a102a71a73a80a17a57a60a64a15a64a79a95a92a57a65a103a87a59a180a62a153a59a122a59a19a209a87a77a17a78a79a59a48a64a79a64a76a71a79a80a17a59 a22a177a80a53a57a92a90a70a59a19a78a76a68a70a78a73a103a87a59a23a78a76a71a73a59a19a78a73a98a99a64 a22a122a72a70a64a91a61a19a68a41a69a53a103a87a57a92a71a79a57a60a68a41a69a40a64a91a57a60a69a25a89a41a68a41a95a60a89a25a57a92a69a53a90 a22a177a95a60a68a54a62a153a59a23a78 a68a41a78a81a103a87a59a23a78a153a71a73a59a19a78a73a98a107a64 a22a74a118 a14a47a161a14a47a161a32a163 a2 a6 a16a48a29a32a28a30a49 a9 a39a27a20a23a26a30a37a43a49a40a37a43a28a7a14a53a55a32a29a114a20a81a18 a14a17a28a7a86a5a4a111a37a162a37a40a20a23a16 a165a43a166a75a167a169a168a47a170a16a15a34a171 a4a91a64a79a59a135a68a41a78a79a71a79a80a17a68a70a90a41a68a70a69a53a72a74a95a52a77a43a68a41a95a60a97a75a69a53a68a41a98a99a57a60a72a41a95a60a64 SMA-HPC ?2000 MIT Normalized 1-D Problem General Quadrature Scheme Using orthogonality and roots () () ( ) ( ) 1 11 10 n nn inini i cxc xdx wcxc x?? = !" () () ( ) ( ) 1 00 10 n nii i cxcxdx wcxcx = =!" 0 0 0 0 a56a58a57a65a61a81a197a107a71a79a80a53a59a67a207a40a202a79a188 a64a117a71a73a68a99a112a40a59a67a198a174a78a79a68a25a68a41a71a73a64a128a68a74a96 a10 a56 a21a114a207 a23 a191a117a80a17a59a135a80a17a57a60a90a41a80a17a59a23a78a63a68a41a78a81a103a87a59a19a78a63a61a19a68a41a69a53a64a66a71a79a78a81a72a74a57a60a69a75a71a73a64a117a72a74a78a73a59a135a59a9a209a17a72a41a61a9a71a79a95a60a97a5a64a79a72a74a71a79a57a65a64a113a208a40a59a23a103 a15 a172a174a173a34a175a177a176a179a178 a22 a191a117a80a17a57a65a64a99a59a19a95a60a59a19a90a75a72a74a69a75a71a99a64a66a71a79a59a23a77a21a78a73a59a19a95a60a57a60a59a23a64a99a68a41a69 a77a43a68a41a95a60a97a75a69a53a68a41a98a99a57a60a72a41a95a117a68a41a78a79a71a79a80a17a68a70a90a41a68a41a69a40a72a74a95a60a57a206a71a113a97a70a118 a104a38a96a149a62a128a59a102a61a81a80a17a68a25a68a70a64a79a59a5a71a73a80a17a59a102a77a43a68a41a95a60a97a70a195 a69a17a68a70a98a180a57a65a72a74a95 a10 a56 a21a114a207 a23a107a64a66a101a40a61a81a80a21a71a73a80a53a72a54a71a111a57a206a71a5a57a65a64a107a68a70a78a66a71a73a80a17a68a41a90a70a68a41a69a53a72a41a95a153a71a73a68 a72a74a95a60a95a76a77a43a68a41a95a60a97a25a69a17a68a41a98a99a57a65a72a74a95a65a64a99a68a74a96a67a57a92a69a17a96a114a59a19a78a73a57a92a68a70a78a85a103a87a59a19a90a70a78a79a59a23a59 a101a114a57a200a118 a59a41a118 a10 a30 a21a32a207a24a23a6a33 a10 a29 a21a114a207 a23a33a88a94a48a94a88a94a33 a10 a56 a11 a29 a21a32a207a24a23a53a102a122a72a74a69a40a103a174a71a79a80a17a59a85a207 a202 a188 a64a185a72a41a78a79a59a107a78a73a68a25a68a74a71a81a64a67a68a74a96a117a71a79a80a53a57a60a64a185a77a40a68a70a95a92a97a25a69a17a68a70a98a180a57a65a72a74a95a200a192a7a71a79a80a17a59a23a69 a210 a189 a71a79a80a53a59a12a80a53a57a92a90a70a80a17a59a19a78a99a68a41a78a81a103a87a59a23a78a180a198 a61a9a68a41a69a40a103a87a57a206a71a73a57a92a68a70a69a53a64a107a72a74a78a73a59a110a72a74a101a87a71a73a68a41a98a107a72a54a71a73a57a60a61a23a72a74a95a60a95a92a97a45a64a73a72a54a71a79a57a65a64a66a208a53a59a23a103a34a118 a99a76a68a41a71a79a59a12a71a79a80a53a72a74a71a99a96a114a68a41a78a99a71a79a80a53a57a60a64 a103a87a59a23a78a79a57a60a89a54a72a54a71a79a57a60a68a41a69a180a62a128a59a63a101a53a64a66a59a48a103a180a77a40a68a70a95a92a97a25a69a17a68a70a98a180a57a65a72a74a95a65a64a46a62a63a80a17a57a65a61a81a80a107a72a74a78a73a59a153a68a70a78a66a71a73a80a17a68a41a90a70a68a41a69a53a72a41a95a75a68a70a69a180a71a79a80a17a59a63a57a60a69a75a71a79a59a19a78a73a89a54a72a74a95a89a44a45a46a33a48a47a6a49a38a118 a115 a101a40a61a81a80 a77a43a68a41a95a60a97a75a69a53a68a41a98a99a57a60a72a41a95a60a64a117a72a41a78a79a59a135a64a66a80a53a57a206a96a169a71a73a59a23a103a110a72a74a69a53a103a12a64a73a61a19a72a74a95a60a59a23a103a5a89a41a59a23a78a73a64a79a57a92a68a70a69a53a64a153a68a74a96a30a71a79a80a17a59a67a61a9a95a65a72a41a64a73a64a79a57a60a61a23a72a74a95a34a187a162a59a23a90a41a59a19a69a40a103a87a78a79a59a15a77a43a68a41a95a60a97a75a69a53a68a41a98a99a57a206a195 a72a74a95a65a64a23a192a87a62a63a80a17a57a60a61a81a80a110a72a74a78a73a59a135a68a41a78a79a71a79a80a17a68a70a90a41a68a70a69a53a72a74a95a43a68a41a69a111a71a79a80a17a59a67a57a60a69a75a71a79a59a23a78a79a89a54a72a74a95a26a44a63a2a47a98a33a48a47a6a49a38a118 a14a47a161a14a47a161a32a163 a16 a2a162a14a17a20a23a29a32a16a48a35a200a18a162a29a32a28a30a49a174a20a23a26a30a22a1a0a94a37a4a3a185a22a75a39 a9 a39a54a86a7a22a75a39 a164 a37a43a28a7a16a23a20a19a39a54a14a17a29a51a28a52a20a23a16 a165a43a166a75a167a169a168a47a170a16a15 a15 a109 a69a102a72a41a112a17a112a17a78a73a59a19a89a25a57a60a72a74a71a79a59a48a103a85a59a9a209a17a72a70a61a177a71a79a69a53a59a23a64a73a64a117a59a23a116a75a101a53a72a54a71a73a57a92a68a70a69 a2a3 a4 a10 a30 a21a114a207 a29 a23 a10a12a10a36a10 a10 a30 a21a114a207 a56 a23 a118 a118 a118 a118 a118 a118 a118 a118 a118 a10 a56 a11 a29 a21a32a207 a29 a23 a10a12a10a36a10 a10 a56 a11 a29 a21a114a207 a56 a23 a5a7a6 a8 a2a3 a4 a20 a29 a118 a118 a118 a20 a56 a5a7a6 a8 a25 a2a3 a4 a5 a29a30 a10 a30 a21a114a207 a23a39a38a70a207 a118 a118 a118 a5 a29a30 a10 a56 a11 a29 a21a32a207a24a23a39a38a70a207 a5a7a6 a8 a99a76a68a54a62 a55a32a29a51a28a7a22a25a14a87a39a25a192a87a207a43a202a66a188 a64a76a72a74a78a73a59a15a197a75a69a53a68a54a62a63a69 a6a63a68a54a62a76a64a117a72a74a78a73a59a67a64a79a72a41a98a99a77a17a95a92a59a48a103a5a68a70a78a66a71a73a80a17a68a41a90a70a68a41a69a53a72a41a95a40a77a43a68a41a95a60a97a25a69a17a68a41a98a99a57a65a72a74a95a65a64a19a118 a172a174a173a34a175a177a176a179a178 a41 a9a128a97a50a101a40a64a66a57a60a69a17a90a36a71a73a80a17a59a5a78a73a68a25a68a74a71a81a64a122a68a74a96 a10 a56 a21a114a207a24a23a122a96a114a68a41a78a181a71a79a80a53a59a5a207a40a202a79a188 a64a19a192a46a71a73a80a17a59a5a80a17a57a60a90a41a80a53a59a19a78a180a68a41a78a81a103a87a59a19a78a185a61a19a68a41a69a53a64a66a71a79a78a81a72a74a57a60a69a75a71a73a64a181a72a74a78a73a59a5a72a74a101a87a195 a71a79a68a70a98a107a72a54a71a79a57a65a61a19a72a41a95a92a95a60a97a36a64a73a72a54a71a79a57a65a64a66a208a53a59a23a103a34a118 a115 a57a60a69a53a61a19a59a181a71a73a80a17a59a107a207a40a202a79a188 a64a122a72a74a78a73a59a181a69a53a68a54a62 a197a25a69a17a68a54a62a63a69a162a192a34a68a70a69a17a95a92a97a36a71a79a80a17a59a107a62a128a59a19a57a60a90a41a80a75a71a73a64a135a72a41a78a79a59a107a64a66a71a79a57a60a95a92a95 a101a17a69a17a197a25a69a17a68a54a62a63a69a47a118a91a191a117a80a17a59a122a95a60a68a54a62a153a59a23a78a91a198a45a61a9a68a70a69a53a64a66a71a79a78a81a72a74a57a60a69a70a71a81a64a91a61a23a72a74a69a102a112a43a59a185a101a53a64a79a59a23a103a110a71a73a68a111a103a87a59a19a71a79a59a23a78a79a98a99a57a60a69a17a59a67a71a79a80a53a59a185a62a128a59a19a57a60a90a41a80a75a71a73a64a23a192a40a90a70a59a19a69a87a195 a59a19a78a81a72a54a71a73a57a92a69a53a90a180a72a180a95a60a57a92a69a17a59a48a72a74a78a76a64a79a97a87a64a113a71a73a59a19a98a110a118 a13a117a138a11a10 a17a84a141a46a147a149a142a25a142a75a144a66a141a143a148 a7a146a147a91a141a10a9a76a158a17a141a30a151a17a147a76a158 a5a15a0a63a147a13a1 a1 a141a46a158a13a12 a14a47a161a3a2a7a161a32a163 a5a111a55a32a49a43a37a53a39a54a29a114a20a48a26a7a42 a2a34a20a48a22a70a44a46a16 a165a43a166a75a167a169a168a47a170a16a15a46a19 a210 a118 a88 a68a70a69a53a64a66a71a79a78a73a101a53a61a177a71a63a198 a44 a47a135a68a70a78a66a71a73a80a17a68a41a90a70a68a41a69a53a72a41a95a40a77a43a68a41a95a60a97a25a69a17a68a41a98a99a57a65a72a74a95a65a64 a27 a29 a30 a10 a202a84a21a32a207a24a23 a10 a2a65a21a114a207 a23a39a38a70a207 a25a32a45 a18a53a52a55a54 a9 a27 a25a134a203 a32 a118 a88 a68a70a98a99a77a17a101a87a71a79a59a91a198a111a39a54a37a162a37a40a20a23a16a70a192a74a207a40a202a53a33a58a203 a25 a47a98a33a48a94a88a94a48a94a6a33a66a198a12a68a74a96a43a71a79a80a17a59a91a198a34a204a32a205a135a68a41a78a81a103a87a59a19a78a143a68a70a78a66a71a73a80a17a68a41a90a70a68a41a69a53a72a41a95a25a77a40a68a70a95a92a97a25a69a17a68a70a98a180a57a65a72a74a95 a64a79a101a53a61a81a80a12a71a79a80a40a72a54a71 a10 a56 a21a114a207a43a202a77a23a26a25a57a45 a189a53a118 a115 a68a41a95a60a89a41a59a67a72a180a95a92a57a60a69a17a59a23a72a41a78a76a64a66a97a87a64a66a71a79a59a23a98 a96a114a68a41a78a117a71a73a80a17a59 a3a185a22a25a29a32a49a43a26a52a20a23a16 a20 a202 a74 a118 a109 a77a17a77a17a78a73a68a27a209a87a57a92a98a107a72a54a71a73a59a15a71a79a80a17a59a67a57a60a69a75a71a79a59a19a90a70a78a73a72a41a95a34a72a41a64a63a72a99a64a66a101a53a98 a5 a29a30 a18 a21a32a207a24a23a39a38a70a207 a25a55a54 a56 a202a58 a29 a20 a202 a18 a21a114a207 a202a23 a172a174a173a34a175a177a176a18a17 a14 a191a117a80a17a57a65a64a67a64a66a95a60a57a65a103a87a59a107a64a66a101a17a98a99a98a107a72a74a78a73a57a60a183a19a59a23a64a15a71a73a80a17a59a180a71a79a59a48a61a81a80a17a69a17a57a65a116a70a101a53a59a107a68a74a96a94a208a40a69a53a103a87a57a60a69a17a90a102a62a153a59a23a57a92a90a70a80a75a71a73a64a135a72a74a69a53a103a50a57a60a69a70a71a73a59a19a90a70a78a73a72a74a71a79a57a60a68a41a69a179a77a43a68a41a57a60a69a75a71a73a64 a96a114a68a41a78a10a4a135a72a41a101a53a64a73a64a63a116a70a101a40a72a41a103a87a78a81a72a54a71a73a101a17a78a79a59a70a118 a210 a74 a14a47a161a3a2a7a161a65a160 a5a99a24a41a24a70a33a47a39a54a14a17a24a41a18 a4a5a22a25a16a48a33a7a55a32a20 a165a43a166a75a167a169a168a47a170a16a15a30a16 a27 a29 a30 a18 a21a32a207a24a23a84a38a41a207 a20 a56 a38 a202a58 a29 a20 a202 a18 a21a32a207a40a202a52a23 a0 a59a19a97a85a77a53a78a79a68a70a77a40a59a23a78a66a71a73a57a92a59a48a64a117a68a74a96a47a71a73a80a17a59a67a98a99a59a9a71a73a80a17a68a87a103 a0 a109 a69a99a198a52a195a200a77a43a68a41a57a60a69a75a71 a4a135a72a41a101a53a64a73a64a143a116a70a101a40a72a41a103a87a78a81a72a54a71a73a101a17a78a79a59a153a78a79a101a53a95a92a59a117a57a65a64a58a22a31a13a162a14a17a24a74a20a117a96a114a68a70a78a143a77a40a68a70a95a92a97a25a69a17a68a70a98a180a57a65a72a74a95a65a64a30a68a74a96a52a68a41a78a81a103a87a59a19a78 a23 a198a4a63 a47 a0 a93a94a78a79a78a73a68a41a78a117a57a65a64a63a77a17a78a73a68a41a77a43a68a41a78a79a71a79a57a60a68a41a69a53a72a41a95a40a71a73a68 a22 a29 a42 a56 a24 a42 a56 a13a117a138a2a1 a0a117a144a2a1 a156a4a3a6a5a8a7 a147a76a141a10a9a76a158a17a141a30a151a17a147a76a158a11a5a12a0a63a150a53a159a13a5a14a1a15a5 a14a47a161a17a16a162a161a32a163 a82 a22a25a28a30a22a70a39a54a14a17a55a117a28a32a31a34a33a149a37a40a29a51a28a52a20a36a35a7a37a53a39a54a42a174a33a7a55a51a14 a165a43a166a75a167a169a168a47a170a16a15 a7 a27 a29 a30 a18 a21a32a207a24a23a39a38a70a207a61a20 a56 a38 a202a58 a29 a47 a198 a18 a22 a203 a63 a29 a42 a198 a24 SMA-HPC ?1999 MIT Normalized 1-D Problem Simple Quadrature Scheme ()f x x0 1 a0 a59a19a97a5a77a17a78a79a68a70a77a40a59a23a78a66a71a113a97a107a68a41a96a7a71a73a80a17a59a67a98a99a59a9a71a73a80a17a68a87a103 a0 a93a94a78a79a78a73a68a41a78a117a57a65a64a63a77a17a78a73a68a41a77a43a68a41a78a79a71a79a57a60a68a41a69a53a72a41a95a40a71a73a68 a29 a56 a14 a18a20a19a22a21 a176a24a23a26a25a28a27a30a29a41a176 a17a50a190a91a68a12a97a41a68a70a101a50a64a79a59a19a59a122a71a73a80a53a72a54a71a135a71a73a80a17a59a180a64a79a57a60a98a180a77a53a95a92a59a99a116a75a101a53a72a41a103a87a78a81a72a54a71a73a101a17a78a73a59a181a64a73a61a81a80a17a59a19a98a99a59a180a57a60a64a135a72a5a64a79a77a43a59a23a61a9a57a65a72a74a95a143a61a23a72a41a64a79a59 a68a74a96a1a4a135a72a74a101a53a64a73a64a63a116a75a101a53a72a41a103a17a78a73a72a74a71a79a101a17a78a73a59 a27 a13a117a138a2a3 a4a50a152 a1 a156a149a141a30a158a40a144a66a148a6a5 a0a117a144a2a1 a156a4a3a6a5a19a7a146a147a91a141a10a9 a141a143a148a13a9 a17a84a141a46a147a149a142a25a142 a7 a147a76a141a10a9 a165a43a166a75a167a169a168a47a170a16a15 a37 SMA-HPC ?1999 MIT Gauss-Quad Simple Quad Number of Points E r r o r Approximately Integrating 1 1 cos2 xdxpi ? " a172a174a173a34a175a177a176a18a17a47a178 a99a76a68a41a71a79a57a65a61a9a59a15a71a79a80a40a72a54a71a63a96a114a68a41a78a76a72a99a64a79a98a99a68a75a68a41a71a79a80a12a96a114a101a17a69a53a61a9a71a79a57a60a68a41a69 a18 a21a114a207 a23 a25 a10 a52 a56 a21 a23a8a7 a207 a23a177a192a53a62a63a80a53a57a60a61a81a80a110a57a65a64a117a57a92a69a87a208a40a69a17a57a206a71a73a59a19a95a60a97a111a103a87a57a92a182a52a59a19a78a73a59a19a69a87a195 a71a79a57a65a72a74a112a53a95a92a59a70a192 a4a135a72a74a101a53a64a73a64a63a116a75a101a53a72a41a103a87a78a81a72a54a71a73a101a17a78a73a59a91a96a32a72a41a78a63a68a41a101a87a71a73a77a40a59a23a78a66a96a114a68a70a78a79a98a107a64a128a71a79a80a53a59a122a64a66a57a60a98a99a77a17a95a92a59a122a116a75a101a53a72a70a103a87a78a81a72a54a71a79a101a53a78a79a59a135a64a79a61a81a80a53a59a19a98a99a59 a210a28a39 a13a117a138a1a0 a4a50a152 a1 a156a149a141a30a158a40a144a66a148a6a5 a0a117a144a2a1 a156a4a3a6a5a19a7a146a147a91a141a10a9 a141a143a148a13a9 a17a84a141a46a147a149a142a25a142 a7 a147a76a141a10a9 a14a47a161a3a2a162a161a32a163 a13a22a3a17a14a17a55a32a33a30a14a87a20a48a29a32a37a43a28 a33a91a37a43a29a32a28a52a20 a33a67a55a32a14a17a24a70a22a25a42a146a22a25a28a52a20 a165a43a166a75a167a169a168a47a170a16a15 a16 SMA-HPC ?1999 MIT Symmetrically Normalized 1-D Problem Comparing Quadratures Evaluation Point Placement -1 1 Simple-Quadrature Points 0 Gauss-Quadrature Points 0-1 1 a99a76a68a41a71a79a57a65a61a9a59a15a71a79a80a53a59a122a61a9a95a60a101a53a64a113a71a73a59a19a78a73a57a92a69a53a90a99a72a74a71a63a57a92a69a75a71a79a59a23a78a79a89a54a72a41a95a34a59a19a69a53a103a17a64 a172a174a173a34a175a177a176a18a17 a17 a104a105a69a146a71a73a80a17a59 a4a135a72a74a101a53a64a73a64a99a116a75a101a53a72a41a103a87a78a81a72a54a71a73a101a17a78a73a59a111a64a73a61a81a80a17a59a19a98a99a59a85a71a73a80a17a59a110a59a19a89a54a72a74a95a60a101a53a72a74a71a79a57a60a68a41a69 a77a43a68a41a57a60a69a75a71a73a64a99a72a74a78a73a59a12a78a79a68a25a68a74a71a81a64a185a68a41a96a149a187a162a59a23a90a41a59a23a69a53a103a87a78a73a59 a77a43a68a41a95a60a97a75a69a53a68a41a98a99a57a60a72a41a95a60a64a128a62a63a80a17a57a65a61a81a80a102a72a74a78a73a59a15a61a19a95a92a101a40a64a113a71a73a59a19a78a73a59a23a103a110a72a54a71a117a71a79a80a53a59a67a59a19a69a53a103a17a64a117a68a41a96a7a71a73a80a17a59a67a57a60a69a70a71a73a59a19a78a73a89a54a72a74a95a200a118 a2 a3a5a4 a11a7a6a107a8a9a10a9a8a185a2a111a6a9a124 a1a11a10 a11a75a1a30a10a12a11a63a6 a9a10a1a46a136a2a11a12a6a19a11a12a3 a10a76a138a66a139 a13a13a12a15a14a122a141a46a156a4a3a66a141a143a150a53a5a17a16a19a18a63a141a10a1 a156a4a3a6a5 a2a7a161a32a163a43a161a32a163 a164 a14a17a55a51a24a41a33a30a55a32a14a87a20a48a29a32a28a30a49a125a20a23a26a30a22a21a20 a2a47a22a25a55a32a35 a31 a106 a22a75a39a54a42a22a20 a165a43a166a75a167a169a168a47a170a16a15 a18 3-D Laplace’s Equation Panel i ic x Collocationpoint , 0 ii ii cc Panel AreaA xx?≈"#$#%One point quadrature Approximation x y z , is an integrable singularity 1 i ii panel i cxx AdS′?=! Basis Function Approach , 1 i ii pa cnel i xx A dS′′?= ! a172a174a173a34a175a177a176a18a17 a85 a99a76a68a54a62a199a95a92a59a19a71a73a64a117a90a41a68a181a112a53a72a41a61a81a197a99a71a73a68a99a68a41a101a17a78a117a77a17a78a73a68a41a112a17a95a60a59a19a98 a68a41a96a7a64a79a68a41a95a60a89a75a57a60a69a17a90a99a187a47a72a41a77a17a95a60a72a70a61a9a59a70a188 a64a153a59a23a116a75a101a53a72a74a71a79a57a60a68a41a69a5a68a70a69a12a72a180a189a41a190 a103a17a68a41a98a107a72a74a57a60a69 a101a53a64a79a57a92a69a53a90a146a112a40a68a70a101a17a69a53a103a17a72a41a78a79a97 a57a60a69a70a71a73a59a19a90a70a78a73a72a41a95a91a78a73a59a19a77a17a78a73a59a23a64a79a59a19a69a75a71a73a72a74a71a79a57a60a68a41a69a162a118 a186a174a59a179a78a73a59a23a72a41a95a92a57a60a183a19a59a36a71a79a80a53a72a74a71a111a62a128a59a179a69a17a68a54a62 a80a53a72a27a89a70a59a179a64a66a68a70a98a99a59 a64a79a68a41a77a17a80a17a57a65a64a66a71a79a57a65a61a19a72a54a71a73a59a23a103a111a71a79a68a25a68a70a95a60a64a76a71a79a68a85a80a53a72a41a69a53a103a87a95a60a59a122a57a60a69a70a71a73a59a19a90a70a78a73a72a41a95a60a64a117a68a41a96a30a96a114a101a53a69a53a61a177a71a73a57a92a68a70a69a53a64a63a71a73a80a53a72a54a71a135a72a74a78a73a59a122a64a66a98a99a68a25a68a74a71a73a80a162a118 a9a153a101a17a71a91a62a63a80a53a72a54a71 a72a74a112a43a68a41a101a17a71a128a71a73a80a17a59a67a57a60a69a70a71a73a59a19a90a70a78a73a72a41a95a43a68a70a69a5a71a73a80a17a59a67a77a53a72a41a69a17a59a19a95a162a62a63a80a53a59a19a78a73a59a91a71a73a80a17a59a122a61a9a59a23a69a70a71a73a78a79a68a70a57a60a103a111a207a24a50a52a51a153a57a65a64a63a95a92a68a87a61a23a72a54a71a79a59a48a103 a27a110a191a117a80a17a59a25a4a15a78a73a59a19a59a23a69a162a188 a64 a96a114a101a17a69a53a61a9a71a79a57a60a68a41a69a102a112a17a95a60a68a54a62a76a64a63a101a17a77a50a72a54a71a76a71a73a80a17a59a185a61a9a59a23a69a75a71a79a78a73a68a41a57a65a103a34a118 a33a76a68a54a62a128a59a19a89a41a59a23a78a23a192a25a71a73a80a17a59a67a96a114a101a17a69a40a61a177a71a79a57a60a68a41a69a179a57a60a64a91a57a92a69a75a71a79a59a23a90a41a78a81a72a74a112a17a95a60a59a67a112a43a59a23a61a19a72a41a101a53a64a79a59 a71a79a80a53a59a179a57a60a69a70a71a73a59a19a90a70a78a73a72a41a69a53a103a157a112a17a95a60a68a54a62a76a64a85a101a17a77a199a72a74a71a12a72a45a78a81a72a54a71a73a59a36a71a73a80a53a72a54a71a12a57a60a64a12a64a66a95a60a68a54a62a128a59a19a78a107a71a73a80a53a72a74a69a84a71a73a80a17a59a179a78a81a72a54a71a73a59a179a72a54a71a12a62a63a80a17a57a65a61a81a80a84a71a79a80a17a59 a64a79a101a17a78a66a96a32a72a70a61a9a59a149a98a180a59a48a72a41a64a79a101a17a78a73a59a91a90a41a68a25a59a48a64a94a71a73a68a181a183a23a59a19a78a73a68a185a57a60a69a5a71a73a80a17a59a135a89a25a57a60a61a19a57a92a69a17a57a92a71a113a97a107a68a74a96a47a71a79a80a17a59a67a64a79a57a92a69a17a90a70a101a17a95a65a72a74a78a73a57a206a71a113a97a70a118 a115 a68a180a62a128a59a149a197a25a69a17a68a54a62a84a71a73a80a53a72a54a71 a210 a35 a71a79a80a53a59a122a57a92a69a75a71a73a59a19a90a41a78a81a72a74a95a34a59a19a209a87a57a60a64a66a71a73a64a76a72a41a69a53a103a12a57a60a64a63a208a53a69a17a57a92a71a79a59a41a192a40a112a17a101a87a71a91a57a65a64 a4a135a72a41a101a53a64a73a64a63a116a70a101a40a72a41a103a87a78a81a72a54a71a73a101a17a78a79a59a67a61a23a72a74a77a53a72a41a112a17a95a92a59a67a68a41a96a30a77a43a59a19a78a79a96a114a68a41a78a73a98a99a57a92a69a53a90 a62a128a59a19a95a60a95a34a57a92a69a12a71a79a80a53a59a67a77a17a78a79a59a48a64a66a59a23a69a53a61a9a59a135a68a41a96a47a71a79a80a53a57a60a64a76a64a79a57a60a69a17a90a41a101a17a95a65a72a74a78a73a57a92a71a113a97 a27 a10a76a138a51a137 a0 a12 a1 a1a15a5a34a151a87a158a53a144a1a0 a5 a9 a139 a12 a16a19a18a117a141a10a1 a156a4a3 a5 a2a7a161a65a160a162a161a32a163 a13 a13a162a14a17a42a45a44a46a55a32a22 a165a43a166a75a167a169a168a47a170a16a15 a4 a5 a29 a11 a29 a29 a2 a3 a60 a3 a38a41a207a61a20 a54 a56 a202a58 a29 a20 a202 a29 a2 a3 a60 a51 a3 The Singular Kernel Problem Symmetrized 1-D Exmaple Note no 0ix =1 x Quad Point x a172a174a173a34a175a177a176a18a17 a26 a104a105a69 a210 a190 a62a128a59a67a95a60a68a75a68a70a197a111a72a74a71a91a72a180a96a114a101a17a69a53a61a177a71a73a57a92a68a70a69 a18 a21a114a207 a23 a25 a29 a2 a3 a60 a3 a62a63a80a17a57a60a61a81a80a110a57a65a64a63a57a60a69a70a71a73a59a19a90a70a78a73a72a41a112a17a95a92a59a135a68a70a69a103a44a63a2a47a98a33a48a47a6a49a162a112a17a101a17a71a91a80a53a72a41a64 a72a99a64a66a57a60a69a17a90a70a101a17a95a60a72a41a78a79a57a92a71a113a97a85a72a54a71a63a207 a25a32a45a17a118 a2a7a161a65a160a162a161a65a160 a2 a28a52a20a48a22a25a49a53a39a54a14 a8a46a55a32a22a36a14a53a28a7a86 a15 a37a43a28a30a29a32a28a52a20a23a22a25a49a40a39a54a14 a8a30a55a51a22 a2a162a29a51a28a7a49a43a33a7a55a51a14a87a39a54a29a114a20a48a29a32a22a25a16 a165a43a166a75a167a169a168a47a170a20a19 a31 ∞→ Areafinite Area → x 1f(x) = x 1f(x) = a10a76a138a18a13 a4a50a152 a1 a156a149a141a30a158a40a144a66a148a6a5 a7 a147a76a141a10a9a76a158a53a141a7a151a17a147a91a158 a5 a0a117a150a53a159a13a5a14a1a15a5a47a142 a165a43a166a75a167a169a168a47a170a20a19a34a171 Gauss-Quad Simple Quad Number of Points E r r o r Approximately Integrating 1 1 1 dx x?" Large errors even with many points! a210 a32 a172a174a173a34a175a177a176a18a17 a11 a186a174a59a76a68a41a112a53a64a79a59a19a78a73a89a41a59a153a71a79a80a40a72a54a71 a4a135a72a41a101a53a64a79a64a94a116a70a101a40a72a41a103a87a78a81a72a54a71a73a101a17a78a79a59a117a57a65a64a143a69a17a68a41a71a58a89a41a59a19a78a73a97a185a90a70a68a25a68a25a103a180a72a54a71a153a57a92a69a75a71a79a59a23a90a41a78a81a72a54a71a73a57a92a69a17a90a15a71a73a80a17a57a60a64a58a96a114a101a17a69a53a61a9a71a79a57a60a68a41a69a162a118 a191a117a80a17a59a135a61a9a68a70a69a75a89a70a59a19a78a73a90a41a59a23a69a53a61a9a59a149a57a60a64a128a78a81a72a54a71a79a80a53a59a19a78a117a77a40a68a25a68a70a78a23a118 a109 a64a117a72a180a98a107a72a54a71a79a71a79a59a19a78a117a68a41a96a47a96a32a72a41a61a9a71a23a192a87a57a92a71a63a57a60a64a117a98a99a68a41a78a73a59a91a57a60a69a53a72a70a61a19a61a9a101a53a78a73a72a74a71a79a59a149a71a79a80a53a72a41a69 a71a79a80a53a59a180a64a79a57a92a98a99a77a17a95a60a59a107a116a75a101a53a72a41a103a87a78a81a72a54a71a73a101a17a78a73a59a181a64a73a61a81a80a17a59a19a98a99a59a41a118a67a104a105a69a50a71a79a80a53a59a181a69a17a59a19a209a25a71a15a96a114a59a23a62 a64a79a95a92a57a65a103a87a59a23a64a15a62a128a59a181a77a53a78a79a59a48a64a66a59a23a69a70a71a67a64a79a59a19a89a41a59a23a78a73a72a41a95a162a71a73a59a23a61a81a80a87a195 a69a17a57a65a116a75a101a17a59a23a64a117a68a41a96a30a80a53a72a41a69a53a103a87a95a60a57a92a69a53a90a180a57a60a69a75a71a79a59a23a90a41a78a81a72a74a95a65a64a153a62a63a57a92a71a79a80a102a64a79a57a92a69a53a90a41a101a17a95a65a72a74a78a73a57a206a71a73a57a92a59a48a64 a101a114a62a63a80a17a57a65a61a81a80a102a72a74a78a73a59a15a57a60a69a75a71a79a59a23a90a41a78a81a72a74a112a17a95a60a59a41a192a25a68a74a96a46a61a9a68a70a101a17a78a73a64a79a59 a102 a0 a115 a101a17a112a17a57a60a69a70a71a73a59a19a78a73a89a54a72a74a95 a101a51a72a41a103a17a72a41a77a87a71a79a57a60a89a41a59 a102a153a116a75a101a53a72a70a103a87a78a81a72a54a71a79a101a53a78a79a59 a0 a88 a80a40a72a74a69a17a90a70a59a15a68a41a96a30a89a54a72a41a78a79a57a65a72a74a112a17a95a60a59a23a64a128a68a41a96a7a57a60a69a70a71a73a59a19a90a70a78a73a72a74a71a79a57a60a68a41a69 a0 a115 a57a92a69a17a90a70a101a17a95a65a72a74a78 a101a4a135a72a41a101a53a64a73a64a66a57a65a72a74a69a96a102a128a116a75a101a53a72a70a103a87a78a73a72a74a71a79a101a17a78a73a59 a10a76a138a11a10 a0 a1 a156a149a158a17a152a2a1 a5a14a9a4a3 a5a7a150a53a159a91a148a149a144a6a5a117a147 a5a7a142 a2a7a161a3a2a7a161a32a163 a2a47a33 a8a46a29a32a28a52a20a48a22a75a39a28a3a17a14a53a55a8a7 a5a99a86a7a14a87a44a30a20a23a29a30a3a40a22a10a9a107a83a107a33a30a14a17a86a7a39a54a14a87a20a23a33a7a39a54a22 a165a43a166a75a167a169a168a47a170a20a19 a15 SMA-HPC ?1999 MIT Improved Techniques x 1 x a115 a101a17a112a52a103a87a57a60a89a25a57a60a103a87a59a135a71a73a80a17a59a67a57a60a69a70a71a73a59a19a90a70a78a73a72a74a71a79a57a60a68a41a69a5a57a60a69a75a71a79a59a23a78a79a89a54a72a74a95 a5 a29 a11 a29 a29 a2 a3 a60 a3 a38a41a207 a25 a5 a11 a30a13a12 a29 a11 a29 a29 a2 a3 a60 a3 a38a70a207a41a44 a5 a30 a11 a30 a12 a29 a29 a2 a3 a60 a3 a38a41a207a10a44 a5 a30a13a12 a29a30 a29 a2 a3 a60 a3 a38a70a207a41a44 a5 a29a30a13a12 a29 a29 a2 a3 a60 a3 a38a70a207 a4a91a64a79a59 a4a135a72a41a101a53a64a79a64 a116a75a101a53a72a41a103a17a78a73a72a74a71a79a101a17a78a73a59a149a57a60a69a110a59a23a72a70a61a81a80a111a64a79a101a17a112a17a57a60a69a70a71a73a59a19a78a73a89a54a72a74a95 a56a46a68a41a95a60a97a25a69a17a68a41a98a99a57a65a72a74a95a65a64a94a208a53a71a91a64a66a101a17a112a53a57a92a69a75a71a79a59a23a78a79a89a54a72a41a95a60a64a117a112a43a59a9a71a79a71a79a59a19a78 a93a143a209a87a77a43a59a19a69a53a64a79a57a92a89a70a59a135a57a206a96a46a98a107a72a74a69a25a97a5a64a79a101a17a112a17a57a60a69a70a71a73a59a19a78a73a89a54a72a74a95a65a64a117a101a53a64a66a59a48a103a34a118 a2a7a161a3a2a7a161a65a160 a164 a26a30a14a17a28a7a49a43a22a36a37a43a35a12a11a67a14a87a39a54a29a32a14 a8a46a55a32a22a25a16 a31a122a35a38a37a53a39 a2a47a29a32a42a146a44a30a55a51a22 a164 a14a17a16a48a22a25a16 a165a43a166a75a167a169a168a47a170a20a19a46a19 a88 a80a53a72a41a69a17a90a41a59a15a89a54a72a74a78a73a57a60a72a41a112a17a95a92a59a48a64a128a71a79a68a99a59a19a95a60a57a60a98a180a57a60a69a53a72a74a71a79a59a122a64a79a57a92a69a17a90a70a101a17a95a65a72a74a78a73a57a206a71a113a97 a5 a42 a25a157a207 a7 a23 a5 a38 a5a100a25a57a38a70a207 a5 a29 a11 a29 a29 a2 a3 a60 a3 a38a41a207 a25 a23 a5 a29a30 a29 a2 a3a13 a14 a3 a23 a5 a38 a5a100a25 a23 a5 a29a30 a23 a38a26a5 a109 a77a17a77a17a95a60a97a18a4a135a72a41a101a53a64a79a64a63a116a75a101a53a72a70a103a87a78a81a72a54a71a79a101a53a78a79a59a15a68a41a69a102a103a87a59a48a64a66a57a60a69a17a90a41a101a53a95a60a72a41a78a79a57a60a183a19a59a48a103a85a57a60a69a70a71a73a59a19a78a73a90a41a78a81a72a74a69a53a103 a210a36a43 a2a7a161a3a2a7a161a14 a2a47a29a32a28a30a49a40a33a30a55a32a14a17a39a12a83a107a33a30a14a17a86a7a39a27a14a17a20a23a33a7a39a27a22a6a31 a164 a37a43a42a45a44a46a55a32a29a32a24a70a14a87a20a23a22a25a86 a164 a14a17a16a23a22a25a16 a165a43a166a75a167a169a168a47a170a20a19a30a16 a104a105a69a75a71a79a59a23a90a41a78a81a72a74a69a53a103a5a80a53a72a41a64a117a197a25a69a17a68a54a62a63a69a110a64a66a57a60a69a17a90a70a101a17a95a60a72a41a78a79a57a92a71a113a97 a56 a21a114a207 a23 a5 a29 a11 a29 a18 a21a114a207 a23 a56 a21a114a207 a23a39a38a70a207a36a62a63a80a17a59a23a78a79a59 a18 a21a114a207 a23a117a57a65a64a63a64a66a98a99a68a25a68a74a71a73a80 a190a149a59a19a89a70a59a19a95a60a68a41a77a110a72a107a116a75a101a53a72a41a103a87a78a81a72a54a71a73a101a17a78a73a59a91a96a114a68a70a78a79a98a181a101a17a95a60a72a181a59a9a209a17a72a70a61a177a71a117a96a114a68a41a78 a5 a29 a11 a29 a14 a79 a21a32a207a24a23 a56 a21a32a207a24a23a84a38a41a207a50a62a63a80a17a59a19a78a73a59a19a14 a79 a21a114a207 a23a117a57a65a64a63a77a40a68a70a95a92a97a25a69a17a68a70a98a99a57a60a72a41a95a43a68a41a96a30a68a70a78a73a103a17a59a19a78a17a18 a88 a72a41a95a60a61a19a101a17a95a65a72a54a71a79a59a135a62a128a59a19a57a60a90a41a80a75a71a73a64a117a95a60a57a60a197a41a59a30a4a135a72a74a101a53a64a73a64a63a116a75a101a53a72a41a103a17a78a73a72a74a71a79a101a17a78a73a59 a172a174a173a34a175a177a176a18a17 a18 a104a38a71a76a57a65a64a117a77a40a68a75a64a79a64a79a112a17a57a60a95a92a59a15a71a79a68a107a90a70a59a19a69a17a59a23a78a73a72a74a71a79a59 a4a135a72a74a101a53a64a73a64a79a57a60a72a41a69a12a116a75a101a53a72a70a103a87a78a73a72a74a71a79a101a17a78a73a59a135a64a79a61a81a80a17a59a23a98a99a59a23a64a128a68a74a96a30a71a79a80a17a59a135a96a114a68a70a78a79a98 a27 a29 a11 a29 a56 a21a114a207a24a23 a18 a21a114a207 a23a39a38a70a207a86a25 a56 a38 a202a58 a29 a20 a202 a18 a21a114a207a43a202a23 a96a114a68a41a78a94a96a114a101a17a69a53a61a9a71a79a57a60a68a41a69a53a64a94a62a63a80a17a57a60a61a81a80a107a80a40a72a27a89a41a59a91a72a135a197a25a69a17a68a54a62a63a69a85a64a79a57a92a69a53a90a41a101a17a95a65a72a74a78a73a57a206a71a113a97 a56 a21a114a207a24a23 a38 a45a53a118a30a191a117a80a17a59a149a116a70a101a40a72a41a103a87a78a81a72a54a71a73a101a17a78a79a59a117a96a114a68a70a78a79a98a181a101a17a95a65a72 a69a17a59a23a59a23a103a17a64a7a71a79a68a15a112a43a59a128a59a19a209a17a72a41a61a177a71a46a62a63a80a17a59a23a69 a18 a21a114a207 a23a7a57a60a64a143a72a149a77a40a68a70a95a92a97a25a69a17a68a70a98a99a57a60a72a41a95a70a68a41a96a53a68a41a78a81a103a87a59a19a78a46a72a54a71a143a98a99a68a70a64a66a71 a18a113a118a34a99a91a68a74a71a58a64a79a101a17a78a79a77a53a78a79a57a65a64a66a57a60a69a17a90a70a95a92a97a70a192 a57a92a71a67a71a79a101a17a78a73a69a53a64a135a68a41a101a17a71a15a71a73a80a53a72a54a71a67a71a73a80a17a59a99a57a92a69a75a71a73a59a19a90a41a78a81a72a54a71a73a57a92a68a70a69a179a77a43a68a41a57a60a69a75a71a73a64a67a72a74a78a73a59a180a71a79a80a17a59a99a198 a78a79a68a25a68a41a71a73a64a15a68a74a96a117a72a111a77a43a68a41a95a60a97a25a69a17a68a41a98a99a57a65a72a74a95 a10 a56 a21a32a207a24a23 a68a74a96a143a103a87a59a23a90a41a78a73a59a19a59a135a198 a25 a21a38a18a11a44a32a47a23a12 a23 a62a63a80a17a57a60a61a81a80a110a57a65a64a63a68a41a78a73a80a17a68a70a90a41a68a41a69a40a72a74a95a43a71a79a68a85a72a41a95a92a95a162a77a43a68a41a95a60a97a25a69a17a68a41a98a99a57a65a72a74a95a65a64a117a68a74a96a46a57a92a69a17a96a114a59a19a78a73a57a92a68a70a78a76a103a87a59a19a90a70a78a79a59a23a59 a62a63a57a92a71a79a80a110a78a79a59a48a64a66a77a43a59a23a61a9a71a117a71a79a68a99a71a79a80a53a59a15a62a128a59a19a57a60a90a41a80a75a71a76a64 a101a169a209 a102a177a192a17a57a200a118 a59a41a118 a27 a29 a11 a29 a56 a21a32a207a24a23 a10 a56 a21a114a207 a23 a10 a2a41a21a114a207a24a23 a25a57a45 a18a53a52a55a54 a9 a25a32a45a46a33a48a47a98a33a48a94a88a94a88a94a88a33a48a21a114a198a21a63 a47a23a94 a109 a69a110a59a19a209a87a72a41a98a99a77a17a95a92a59a135a57a65a64a128a71a79a80a17a59a122a64a79a57a60a69a17a90a41a101a17a95a65a72a74a78a63a57a60a69a75a71a79a59a23a90a41a78a81a72a74a95 a0a22a25 a27 a29 a11 a29 a18 a21a32a207a24a23 a2 a47 a63a50a207 a42 a33a76a59a23a78a79a59a70a192 a56 a21a114a207 a23 a25a10a47 a12 a2 a47 a63a50a207 a42 a72a74a69a53a103a185a71a79a80a17a59a63a61a19a68a41a78a73a78a79a59a48a64a66a77a43a68a41a69a53a103a17a57a92a69a17a90a149a68a41a78a79a71a79a80a17a68a70a90a41a68a70a69a53a72a74a95a75a77a40a68a70a95a92a97a25a69a17a68a70a98a99a57a60a72a41a95a60a64a7a71a79a101a53a78a79a69a180a68a41a101a87a71 a71a79a68a179a112a40a59a107a71a73a80a17a59 a88 a80a17a59a23a112a25a97a25a64a79a80a17a59a23a89a50a77a43a68a41a95a60a97a25a69a17a68a41a98a99a57a65a72a74a95a65a64a19a118a110a191a117a80a17a59a85a57a60a69a75a71a79a59a19a90a70a78a73a72a74a71a79a57a60a68a41a69a174a77a40a68a70a57a92a69a75a71a81a64a185a72a74a78a73a59a107a90a41a57a60a89a41a59a23a69a45a57a92a69a146a61a19a95a92a68a75a64a66a59a48a103 a96a114a68a41a78a73a98a3a112a25a97 a207 a202 a25 a10 a52 a56 a22 a7 a23 a203a25a63 a47 a23 a198 a24 a72a74a69a40a103a85a71a73a80a17a59a122a61a9a68a70a78a79a78a73a59a23a64a79a77a40a68a70a69a53a103a87a57a60a69a17a90a180a62a153a59a23a57a92a90a70a80a70a71a81a64a117a72a74a78a73a59 a20 a202a34a25 a7 a12a27a198a143a118 a210 a20 a2a7a161a3a2a7a161a3a2 a2a47a29a32a28a30a49a40a33a30a55a32a14a17a39a12a83a107a33a30a14a17a86a7a39a27a14a17a20a23a33a7a39a27a22a42a0 a22a25a29a32a49a43a26a52a20a23a16 a165a43a166a75a167a169a168a47a170a20a19 a7 a2a3 a4 a10 a30 a21a32a207 a29 a23 a10a12a10a36a10 a10 a30 a21a114a207 a56 a23 a118 a118 a118 a118 a118 a118 a118 a118 a118 a10 a56 a11 a29 a21a114a207 a29 a23 a10a12a10a36a10 a10 a56 a11 a29 a21a32a207 a56 a23 a5a7a6 a8 a2a3 a4 a20 a29 a118 a118 a118 a20 a56 a5a7a6 a8 a25 a2a3 a3 a4 a5 a29 a11 a29 a10 a30 a21a32a207a24a23 a56 a21a32a207a24a23a84a38a41a207 a118 a118 a118 a5 a29 a11 a29 a10 a56 a11 a29 a21a114a207a24a23 a56 a21a114a207a24a23a84a38a41a207 a5 a6 a6 a8 a99a76a59a23a59a23a103a103a101a32a72a41a69a53a72a74a95a60a97a75a71a79a57a65a61a48a102a94a96a114a68a41a78a73a98a185a101a53a95a60a72a70a64a94a96a114a68a70a78a63a57a92a69a75a71a73a59a19a90a41a78a81a72a74a95a65a64a128a68a74a96 a10 a21a114a207 a23 a56 a21a114a207 a23 a172a174a173a34a175a177a176a18a17 a30 a191a117a80a17a59a135a95a60a68a54a62a153a59a23a78a128a68a70a78a73a103a17a59a19a78a63a61a9a68a70a69a53a64a113a71a73a78a73a72a41a57a92a69a75a71a81a64a63a61a19a72a74a69a12a112a43a59a67a101a53a64a66a59a48a103a85a71a73a68a85a61a9a68a70a98a180a77a53a101a87a71a79a59a15a71a79a80a53a59a67a57a92a69a75a71a79a59a23a90a41a78a81a72a54a71a73a57a92a68a70a69a5a62a128a59a19a57a60a90a41a80a75a71a73a64a23a118 a0 a6a107a2 a3 a3 a124a2a1a2a1 a165a43a166a75a167a169a168a47a170a20a19 a37 a13a15a14a17a16a19a18a21a20a23a22a25a24a27a26a7a28a30a29a32a31a34a33a30a22a36a35a38a37a40a39a5a24a41a37a43a42a45a44a46a33a47a20a48a29a32a28a30a49a50a29a51a28a52a20a23a22a25a49a53a39a54a14a17a55a51a16 a56a58a57a60a59a23a61a19a59a19a62a63a57a65a64a66a59a67a61a9a68a70a69a53a64a66a71a73a72a74a69a75a71a76a72a74a77a53a77a17a78a79a68a75a72a41a61a81a80 a82 a14a17a33a30a16a23a16a48a29a32a14a53a28 a31a34a33a30a14a17a86a7a39a54a14a87a20a23a33a7a39a54a22 a24a53a72a70a64a113a71a73a59a19a78a76a61a19a68a41a69a25a89a41a59a23a78a79a90a70a59a19a69a53a61a19a59 a93a94a64a73a64a66a59a23a69a75a71a79a57a65a72a74a95a34a78a73a68a41a95a60a59a15a68a41a96a7a68a70a78a66a71a73a80a17a68a41a90a70a68a41a69a53a72a41a95a40a77a43a68a41a95a60a97a25a69a17a68a41a98a99a57a65a72a74a95a65a64 a106 a22a75a24a27a26a30a28a7a29a51a31a34a33a7a22a25a16a107a35a38a37a40a39a5a16a23a29a51a28a7a49a43a33a30a55a32a14a87a39a107a108a43a22a70a39a54a28a30a22a25a55a32a16 a109 a103a17a72a41a77a87a71a73a72a74a71a79a57a60a68a41a69a110a72a74a69a53a103a4a3a58a72a74a78a73a57a60a72a41a112a17a95a60a59a135a191a7a78a81a72a74a69a53a64a66a96a114a68a41a78a73a98a107a72a54a71a79a57a60a68a41a69 a115 a57a60a69a17a90a41a101a53a95a60a72a41a78a117a116a75a101a53a72a70a103a87a78a81a72a54a71a79a101a53a78a79a59 a0 a26a7a14a17a20a99a14 a8a94a37a40a33a7a20a107a42a50a33a30a55a114a20a48a29a114a44a46a55a32a22a12a86a7a29a51a42a146a22a25a28a7a16a48a29a32a37a43a28a30a16a10a6 a32 a0