Hi everyone,
I'm trying to understand how bvodeS works and whether I can use it to
solve numerically my problem.
A bit of context:
I have a 1D coupled-mode problem, where I model the propagation along
one direction (let's call it y) of two modes (one in +y-direction and
one in the
in zeta you will give the locations (relative to y) of your boundary
conditions
in your case something like [0 0 0 , L L L ]
if you have 6 states 3 with initial constraints and 3 with final contraints
The constraints should then be defined by the gsub and dgsub functions
Serge
Le 12/11/2015
