[apciat]

x=0.5 : 1 : p+0.5;
y=D(1) : D(p+1);
plot(x,y,'g-x')

其中p是上面程序中算出来的自然数；D是上面程序的计算结果，1*p+1的矩阵

程序运行正常，唯一的问题就是最后应该画拆线图，出来的却是散点图，在下是初学者，不明何故，求高手高见^_^