system of coupled ordinary differential equations

i need matlab code by using bvp4c for the above equations.....

(1+1/beta)f'''(eta)-2*(f'+g')f'+(f+g)f''-M^2f'=0

(1+1/beta)g'''(eta)-2(f'+g')g'+(f+g)g''-M^2g'=0

teta''+pr[(f+g)theta'-A(f'+g')theta+S*theta]=0

phi''+Sc[(f+g)phi'-B(f'+g')phi-r*phi]=0

boundary conditions

f(0)=0 ,f'(0)=1 ,g(0)=0, g'(0)=alpha ,theta(0)=1+r1theta'(0) ,phi(0)=1+r2phi'(0)

f'(infinity)=0 g'(infinity)=0,theta(infinity)=0,phi(infinity)=0

[Stephan]

What have you tried so far?

[Priya M]

xlow=0;

xhigh=6;

solyinit=bvpinit(linspace(xlow,xhigh,2),[0,1.01,11,1,0,1,0,12,11.05,1.005]);

M=0.5;

beta=0.02;

B=0.5;

Pr=0.7;

A=0.5;

Sc=1.0;

S=0.1;

r=0.3;

r1=0.1;

r2=0.1;

dydx=@(x,y)[y(2);

y(3);

(((M^2*y(2))-2*y(2)^2-2*(y(2)*y(5)-y(1)*y(3)-y(4)*y(3))/(1+1/beta)));

y(5);

y(6);

(((M^2*y(5))-2*(y(2)*y(5))-2*y(5)^2-y(1)*y(6)-y(4)*y(6))/(1+1/beta));

y(8);

((1/Pr)*(A*(y(2)*y(7)-y(1)*y(8)-y(4)*y(8)+A*(y(5)*y(7))-S*y(7))));

y(10);

((1/Sc)*(B*(y(2)*y(9))+B*(y(5)*y(9))-y(1)*y(10)-y(4)*y(10)+r*y(7)))];

res=@(ya,yb)[ya(1); ya(2)-1; ya(4); ya(5)-alpha; ya(7)-r1*ya(8); ya(9)-r2*ya(10); yb(2); yb(5); yb(7); yb(9)];

fsol=bvp4c(dydx,res,solyinit);

X=fsol.x;

Y=fsol.y;

plot(X,Y)