此微分方程可以用数值方法求得。如变步长四、五阶龙格-库塔法，实现过程：
1、自定义微分方程求解函数
function f= func(x,y)f(1)=y(2);f(2)=0.5*sqrt(1+y(1)^2)/(1-x);f= f(:);end
2、在窗口命令输入
x0=[0,0];[x,y] = ode45(@func,[0 10],x0)
plot(x,y)
3、运行结果
