>> f

f =
    @(t,x,v)[4*x(1)+v(1)*x(1)*x*(2);-2*x(2)+v(2)*x(1)*x(2)]

>> g

g =
    @(ya,yb,v)[ya(1)-2;ya(2)-1;yb(1)-4;yb(2)-2]


>> sol=bvp5c(f,g,solinit);

??? Error using ==> bvparguments at 121
Error in calling BVP5C(ODEFUN,BCFUN,SOLINIT):
  The derivative function ODEFUN should return a column
  vector of length 2.

Error in ==> bvp5c at 137
[neqn,nparam,nregions,atol,rtol,Nmax,xyVectorized,printstats]
= ...

>> whos

  Name         Size            Bytes  Class              Attributes

  f            1x1                16  function_handle              
  g            1x1                16  function_handle              
  solinit      1x1              1146  struct                       
  x1           2x1                16  double                       
  x2           2x1                16  double         

dydt=@(x,y)[y(2); -6*y(2)+x];
bc=@(ya,yb)[ya(1)+1; yb(1)-4];
yinit=@(x)[exp(-x); -exp(-x)];
solinit=bvpinit(linspace(0,1,20),yinit);
sol=bvp4c(dydt,bc,solinit);
% "sol" callback you three part:x,y,y'
plot(sol.x,sol.y(1,:))