function dydx = bvpfcn(x,y) % equation to solve
dydx = zeros(2,1);
dydx(1)=0.2-5*y(1)-0.1*(y(1)^2)-2*y(1)*y(2)
dydx(2)=0.22-5*y(2)-0.12*(y(2)^2)-2.5*y(1)*y(2)
end