Skip to main content

j0_maclaurin_series_fast

Function j0_maclaurin_series_fast 

Source
pub(crate) fn j0_maclaurin_series_fast(x: f64) -> f64
Expand description

Generated by SageMath:

mp.prec = 180
def print_expansion_at_0():
    print(f"const J0_MACLAURIN_SERIES: [(u64, u64); 12] = [")
    from mpmath import mp, j0, taylor
    poly = taylor(lambda val: j0(val), 0, 24)
    real_i = 0
    for i in range(0, 24, 2):
        print_double_double("", DD(poly[i]))
        real_i = real_i + 1
    print("];")
    print(poly)

print_expansion_at_0()