Define ‘a’ in your workspace, then use this bvp4c call:
sol = bvp4c(@(x,y)twoode(x,y,a),@(ya,yb)twobc(ya,yb,a),solinit);
You are passing the extra parameter correctly, however bvp4c does not need to know about it, so only show bvp4c the arguments it wants, thus:
The anonymous function construction (that you wrote correctly) will pick up ‘a’ from your workspace and pass it to your functions.
With this change (and specifying ‘a’), your code ran without error with this bvp4c call when I tested it.