原函数（求积分）是计算不了的，因为不唯一。以下代码为花本身函数及其导函数的方法。其中红线为原函数，绿线为导函数p=[2 3 0 4];q1=polyder(p);x=-3:0.1:3;plot(x,sum([x'.^3 x'.^2 x' ones(size(x'))]*p',2),'r-',x,sum([ x'.^2 x' ones(size(x'))]*q1',2),'g-');
