<gh_stars>0
#%%
#!/usr/bin/env python

# 定積分の求め方
# 参考　https://atarimae.biz/archives/22555

import sympy as sym

x, y = sym.symbols("x y")


# 積分
def IndefiniteIntegral(y):
    return sym.Integral(y).doit()


# 定積分（密度関数） コード上 cは無視する
# a to b: integral(f(x)dx + c = g(x))
# 解 g(a) - g(b)
def DefiniteIntegral(formula, a, b):
    y = formula
    indefiniteIntegral = IndefiniteIntegral(y)
    a = indefiniteIntegral.subs(x, a)
    b = indefiniteIntegral.subs(x, b)

    return b - a


def main():
    formula = 4 * x**3
    definiteIntegral = DefiniteIntegral(formula, 1, 2)

    print(definiteIntegral)


if __name__ == "__main__":
    main()

# %%

# %%

#%%
#!/usr/bin/env python

# 定積分の求め方
# 参考　https://atarimae.biz/archives/22555

import sympy as sym

x, y = sym.symbols("x y")


# 積分
def IndefiniteIntegral(y):
    return sym.Integral(y).doit()


# 定積分（密度関数） コード上 cは無視する
# a to b: integral(f(x)dx + c = g(x))
# 解 g(a) - g(b)
def DefiniteIntegral(formula, a, b):
    y = formula
    indefiniteIntegral = IndefiniteIntegral(y)
    a = indefiniteIntegral.subs(x, a)
    b = indefiniteIntegral.subs(x, b)

    return b - a


def main():
    formula = 4 * x**3
    definiteIntegral = DefiniteIntegral(formula, 1, 2)

    print(definiteIntegral)


if __name__ == "__main__":
    main()

# %%

# %%
#%% #!/usr/bin/env python # 定積分の求め方 # 参考　https://atarimae.biz/archives/22555 # 積分 # 定積分（密度関数） コード上 cは無視する # a to b: integral(f(x)dx + c = g(x)) # 解 g(a) - g(b) # %% # %%