Solve a BVP with Python (e.g., 'solve_bvp') with arbitrary precision (e.g. 'mpmath')

Question

I am trying to solve a BVP using \'solve_bvp\'. Below is the code:

import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import solve_bvp
# Inp