Hi, I'm solving a system of ODEs with BCs, with bvp4c function.
My simplified script is as follows:
rmesh = ...
solinit = ...
sol = bvp4c(@bvpfcn,@bcs,solinit);
function dydr = bvpfcn(r,y)
function res = bcs(ya,yb,C)
In this case, I solved this system with prefixed values of BCs, reported inside of the second function with the vector C. Thus, I have only one solution, which is related to those BCs values.
However, I'd like to solve this problem with a for loop, through different vectors as input in BCs function in order to have multiple distict solutions, one for every BCs I inserted.
Is this the correct way?
sol = ??
Obviously, I should pass one BC at a time, but I'm interested to do it all together.