এসসিপি

from scipy import integrate

result, error = integrate.quad(
lambda x: x ** 2,
    0,
    3
)

print(result)
print(error)

from scipy import integrate

result, error = integrate.dblquad(
    lambda y, x: x * y,
    0,
    2,
    lambda x: 0,
    lambda x: x
)

print(result)
print(error)

from scipy import integrate

result, error = integrate.tplquad(
    lambda z, y, x: x * y * z,
    0,
    1,
    lambda x: 0,
    lambda x: x,
    lambda x, y: 0,
    lambda x, y: x + y
)

print(result)
print(error)

from scipy import integrate

result, error = integrate.nquad(
    lambda x, y: x ** 2 + y ** 2,
    [[0, 1], [0, 2]]
)

print(result)
print(error)

from scipy import integrate

result, error = integrate.fixed_quad(
    lambda x: x ** 2,
    0,
    3,
    n=5
)

print(result)

from scipy import integrate

data = [0, 1, 4, 9, 16, 25, 36]

result = integrate.simpson(
    data,
    x=[0, 1, 2, 3, 4, 5, 6]
)

print(result)

from scipy import integrate

data = [0, 1, 4, 9, 16, 25, 36]

result = integrate.trapezoid(
    data,
    x=[0, 1, 2, 3, 4, 5, 6]
)

print(result)

from scipy import integrate

data = [0, 1, 4, 9, 16, 25, 36]

result = integrate.cumulative_trapezoid(
    data,
    x=[0, 1, 2, 3, 4, 5, 6],
    initial=0
)

print(result)

from scipy import integrate

def equation(t, y):
    return [-0.5 * y[0]]

result = integrate.solve_ivp(
    equation,
    [0, 10],
    [10],
    t_eval=[0, 2, 4, 6, 8, 10]
)

print(result.t)
print(result.y[0])

from scipy import integrate

def equation(y, t):
    return -0.5 * y

time = [0, 2, 4, 6, 8, 10]

result = integrate.odeint(
    equation,
    10,
    time
)

print(result)

from scipy import integrate

def equation(x, y):
    return [
        y[1],
        -y[0]
    ]

def boundary(ya, yb):
    return [
        ya[0],
        yb[0]
    ]

x = [0, 1, 2, 3]
y = [
    [0, 0, 0, 0],
    [0, 0, 0, 0]
]

result = integrate.solve_bvp(
    equation,
    boundary,
    x,
    y
)

print(result.success)
print(result.y[0])