bvp4c problem, unbounded shear layer

Answer 1

[darova]

You need 2 initial conditions (you have 2th order ODE)

function yinit = mat4init(x)
yinit = [ cosh(x)+sinh(x) ];
end