Hi,

I need some help.  I have a system of ODE's subject to a system of Boundary Conditions.  I can't figure out the sequence of commands to get MAPLE to generate a solution.

ODE System:

SYS1 := {diff((C[10])(z), `$`(z, 2)) = (BR[1]+BR[2])*(C[10])(z)-BR[3]-(BR[4]+BR[5])*C[50], diff((C[51])(z), `$`(z, 2)) = Ped*u[0]*C[50*x]-BR[5]*(C[10])(z)+BR[3]*C[50]}

Boundary Conditions:  I am going to be evaluating from z=s to z=h

BC1 := {(C[10])(h) = 1, diff((C[10])(s), s) = 0, diff((C[51])(h), h) = 0, diff((C[51])(s), s) = Da[1]*(C[50]-J3)}

Given a cubic, x^3 - c*x +16=0, I need to find values of c so that the real parts of the complex roots are negative.  Any suggestions?