需要编写一个matlab的m文件，用于解决多项式

制作 一个m文件，可以达到输入一个多项式，得出多项式的极值，单调性，凹凸位置等等。
问题在于（多项式案例：y=1/5（x^4-6x^2+8x+7)）

1.最后的答案不对（程序运行后的结果）
函数在区间（-inf,-3.2028])上单调减少.
函数在区间[-3.2028,+inf])上单调增加.
函数在x=-3.2028处取得极大值5.0106
函数在x=-3.2028处取得极小值5.0106
函数在区间(-inf,-1]上为凹.
函数在区间[-1,1]上为凸.
2.显示我的程序出错，但是找不到问题。
? Attempted to access r2(3); index out of bounds because numel(r2)=2.

Error in ==> polyp at 117
disp(['函数在区间[',num2str(r2(k-1)),',',num2str(r2(k)),']上为凹.']);

我的m文件