1、建立自定義函數func（）
function
f
=
func(t,x)
%x''+0.7x'+0.8x'|x'|+25.6x-25.6x?=0
f(1)=x(2);
f(2)=25.6*x(1)^3-25.6*x(1)-0.8*x(2)*abs(x(2))-0.7*x(2);
f=f(:);
2、建立龍格庫塔算法函數runge_kutta（）
調用格式：[t,x]
=
runge_kutta(@(t,x)func(t,x),x0,h,a,b);
3、然後根據x和x'數據，繪制出x(t)、x′(t)的圖形。
plot（x（：，1）,x（：，2））