Problem-solving strategies for boundary value problems in odes pdes
Use this skill when working on boundary-value-problems problems in odes pdes.
Problem Classification
Shooting Method
scipy.integrate.solve_ivp + root findingFinite Difference Method
sympy_compute.py linsolve "tridiagonal_matrix" "boundary_vector"Collocation/BVP Solver
scipy.integrate.solve_bvp(ode, bc, x, y_init)Eigenvalue Problems
sympy_compute.py eigenvalues "sturm_liouville_matrix"uv run python -c "from scipy.integrate import solve_bvp; import numpy as np; ode = lambda x, y: [y[1], -y[0]]; bc = lambda ya, yb: [ya[0], yb[0]-1]; x = np.linspace(0, np.pi, 10); y = np.zeros((2, 10)); sol = solve_bvp(ode, bc, x, y); print('Solution at pi/2:', sol.sol(np.pi/2)[0])"
uv run python -m runtime.harness scripts/sympy_compute.py linsolve "tridiagonal_matrix" "boundary_vector"
uv run python -m runtime.harness scripts/z3_solve.py prove "eigenvalue_real"
From indexed textbooks:
See .claude/skills/math-mode/SKILL.md for full tool documentation.