এসসিপি solve_bvp

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])