1.降阶
z1=y'
z2=y
==>
z1'=y''=-y^2=-z2^2
z2'=y'=z1
2.求解
function  dz = confun2(x,z)
dz(1)=-z(2)^2;
dz(2)=z(1);

[x,num_y]=ode23(@confun2,[0 20],[0 1])
3.画图
subplot(2,1,1) ,plot(x,num_y(:,1))
title('1st der of y''),xlabel('x'),grid
subplot(2,1,2) ,plot(x,num_y(:,2))
title('y'),xlabel('x'),grid