Multi initial gusset at bvp

Hello I want to solve ode system at two domain by bvp4c at matlab and I need two initially gusset for two domain because any domain have different value and I must writ two initially gusset for bvpinit and I don't known what do I do about this problem.

Answer 1

[Torsten]

Use an if-statement depending on the spatial coordinate:

Example:

solinit = bvpinit(linspace(0,pi,10),@mat4init)

function yinit = mat4init(x)

if x<pi/2
   yinit = 0;
else
   yinit = 1;
end

Best wishes

Torsten.

Answer 2

[Mostafa Mostafa]

Thank you. Your answer is very useful I have another question

I want solve two systems equation by bvp4c in matlab at one domain. the two systems couple together and I must be solved simultaneously and one system is nonlinear algebraic equation system and another system is ordinary diffrantioal equation system.how do I define nonlinear algebraic equation system at bvp4c in matlab ???? Please write an example about my problem I wish best you and your company (like)

[Torsten]

As far as I know, bvp4c can not solve a mixture of differential and algebraic equations.

Maybe you could post your equations to see if there is another possibility for solving.

Best wishes

Torsten.