e2d03438-d39e-4017-af33-b70cac58023c 5 6 12 decimal

&[DATE] &[TIME] - &[FILENAME]

n 3 :

δ 0.1 :

n κ b ω H B b 1 n range el vectorize 11 line : κ 2^ω n^1 - (n 1 - 2 / k -^B 2 k * el * ω 2 k * 1 -^* (k 1 n 1 - 2 / sum - (2^1 - (n 1 - 2 / k -^B 2 k * 1 + el * ω 2 k *^* (k 0 n 1 - 2 / sum (2^+ / 21 line :

n κ b ω H κ 2^ω 3^1 - (31 - 2 / 1 -^b 2 el * ω 21 * 1 -^* - (2^1 - (31 - 2 / 0 -^b 1 el * ω 20 *^* 1 - (31 - 2 / 1 -^b 3 el * ω 21 *^* + (2^+ /

μ b ω expand μ 2^ω 3^1 - (31 - 2 / 1 -^b 2 el * ω 21 * 1 -^* - (2^1 - (31 - 2 / 0 -^b 1 el * ω 20 *^* 1 - (31 - 2 / 1 -^b 3 el * ω 21 *^* + (2^+ / :

b 12313 mat transpose :

μ 1.23456789 :

x f μ b x expand :

h 0.001 :

x h f D1 x 2 h * + f - 8 x h + f * + 8 x h - f * - x 2 h * - f + 12 h * / eval :

μ b x expand x h f D1 2 1 sys

x f x 2 - 2 int 2.862