为什么解析`dsolve`与`dsolve`数字结果有很大的不同？

Question

这是Windows上的Maple 2016.2。

同样的ODE。从分析dsolve绘制的解决方案看起来与Maple的数值解决方案有很大不同。差异显示在域的右侧。我怀疑数值解是正确的，无论是绘制解析解，还是给出的解析解都不正确？怎么解释呢？

restart;
eps:=1/100:
ode:=eps*diff(y(x),x$2)+(x^2+1)*diff(y(x),x)-x^3*y(x)=0;
sol:=dsolve({ode,y(0)=1,y(1)=1},y(x)):
plot(rhs(sol),x=0..1,y=0..1);

这是数值解图

sol:=dsolve({ode,y(0)=1,y(1)=1},y(x),range=0..1,'numeric',maxmesh=500);
plots:-odeplot(sol, x=0..1,view=[0..1,0..1]);

顺便说一句，Maple的数值解也同意Mathematica的数值解法。

问题：哪个是正确的情节？从分析或数值解决方案？ Maple分析解决方案不正确吗？或者我在某个地方犯了一些错误？