f(x)=x^3
_m=EvalFunc(f( x ),x,-5,5,30)
PlotData(_m[0],_m[1])
f(x,y0,y1)=y1
g(x,y0,y1)=-y0
m=MMiddlePoint([f(x,y0,y1);g(x,y0,y1)],[1,0],0,6,30)