matlab解二阶微分方程>>

我解的方程是光在非均匀介质中的传播,介质折射率为位置的函数n=2-(x^2+y^2)

用费马定理得到二阶微分方程y''+2(y+xy')(1+y'^2)/(2-x^2-y^2)=0

我把它转换为两个一阶微分方程,t=x,y(1)=y,y(2)=y',

所以y(2)'=2(y(1)+ty(2))(1+y(2)^2)/((t^2+y(1)^2)-2)]

其中边界条件为x=-1,y=0

从上面的思路,

需要解的就是这个二阶微分方程,下面是matlab语句,报错

>>clear;close;

>>sinit=bvpinit(-1:1,[1;0]);

>>odefun=inline('[y(2);2*(y(1)+t*y(2))*(1+y(2)^2)/((t^2+y(1)^2)-2)]','t','y');

>>bcfun=inline('[ya(1)+1;yb(1)]','ya','yb');

>>sol=bvp4c(odefun,bcfun,sinit)

运行结果报错：

Errorusing==>bvp4cat203

Unabletosolvethecollocationequations--asingularJacobianencountered

这个二阶微分方程该怎么解呢,

请给正确的语句,