Solving Boundry value problem with fsolve

[Walter Roberson]

http://www.mathworks.com/matlabcentral/answers/6200-tutorial-how-to-ask-a-question-on-answers-and-get-a-fast-answer

[Alex Sha]

Hi, Justin, try to solve your BVP problems by using 1stOpt:

Code (Case-1):

ODEStep = 0.05;
Constant B=-0.056, Peh=2, Pem=2, r=17.6;
ConstStr Rct=3.36*c*exp(r-r/T);
SubjectTo c'[0]=(c[0]-1)*Pem, T'[0]=(T[0]-1)*Peh;
Variable x,c',T';
ODEFunction c''=(c'+Rct)*Pem;
            T''=(T'+B*Rct)*Peh;
Data;
0 NAN NAN
1 0 0

Result (Case-1):

Root of Mean Square Error (RMSE): 3.98105846789366E-10
Sum of Squared Residual: 6.33953060991512E-19

Parameter          	Best Estimate
--------------------	-------------
c' Initial Value 	-1.15593943548473
t' Initial Value 	0.0647326083792722
c Initial Value 	0.422030282257633
t Initial Value 	1.03236630418964

SubjectTo:
t'[initialvalue]-((t[initialvalue]-1)*2): -2.77555756156289E-16
c'[initialvalue]-((c[initialvalue]-1)*2): 3.88578058618805E-15

chart of x-c

chart of x-t

chart of x-dc/dx

chart of x-dt/dx