Section: Optimization and Curve Fitting
pint = polyint(p,k)
 or, for a default k = 0,
pint = polyint(p);
 where p is a vector of polynomial coefficients assumed to be in
 decreasing degree and k is the integration constant.
 Contributed by Paulo Xavier Candeias under GPL
polyint.
--> polyint([2,3,4])
ans = 
    0.6667    1.5000    4.0000         0 
And
--> polyint([2,3,4],5)
ans = 
    0.6667    1.5000    4.0000    5.0000