Shooting method with a constraint

Hi,

I want to solve a system of ODE. It's a BVP. I changed it to IVP and used Shooting method to solve this problem. For the method I defined

    X_i_new = X_i - fX_i.*((X_i - X_i_1)/(fX_i - fX_i_1));

X_i and X_i_1 are matrixes like [0.018 0.0258 0.9562] and [0.0150 0.0259 0.9591]. The sum of X_i elements is unity, and also for X_i_1.

The question is X_i_new calculates the new X for the 3rd loop. But the sum of X_i_new does not equal unity. I don't know how to put this criterion in the shooting method.

I am not sure but I changed the shooting method like:

    X_i_new = [X_i_new(1) X_i_new(2) 1-(X_i_new(1)+X_i_new(2))];

was that the correct change?

Thank you, Regards

[Torsten]

It would be helpful if you presented the BVP you try to solve.