[御景风yjf]

要求解的微分方程是y''-2*10^(-6)/y^3-1.56*10^(-11)/y=0 ,用ode15和ode45求解图像不对，四阶RK法应该是可以求解的但是改了半天还是有问题，求大神帮忙求解下。初值条件为y（0）=0.1，y（0）'=0