Solve BVP with Arrays as Coefficients

Can bvp4c and/or bvp5c be used to solve a boundary value problem where the coefficients (or terms) of the BVP are provided by an array, versus symbolic (e.g. x, x^2,...)?

For example, I need to solve:

y'' + a*y' + b*y = f. where, say x, is defined between 0 and 1. A, b and f are vectors/arrays at certain values of x. Can Matlab's BVP solvers interpolate these arrays as needed?

Thanks!

Answer 1

[Torsten]

If you have functions depending on x for a, b and f, you are done.

If not, you can try the following:

Use "interp1" to interpolate a, b and f in the x-value where the coefficients are needed by bvp4c.

Or fit functions to a, b and f with x as independent variable before calling bvp4c and use these function to supply a, b and f at the relevant x-positions.