Question: Boundary conditions not applied correctly in numerical solution of ODE system?

Hi everyone,

I am trying to solve a system of coupled ODEs numerically. My worksheet runs and produces results, but when I look at the graphs, it seems that the boundary conditions are not being satisfied correctly. Could anyone help me identify the issue and  fix my implementation?

BCs_help.mw

Please Wait...