matlab怎么对数值积分结果再次积分?

您好,我需要求一个分段函数的数值积分,求完后对数值积分结果再次求积分(类似于对加速度曲线积分得到速度,再对速度曲线积分得到位移曲线),代码该怎么写呢。
下面的代码是求第一次数值积分,已经运行成功(用的是您的方法)。接下来怎么对第一次数值积分得到的函数p再求积分呢?
g=@(t)(t).*(t<1)+1.*(1<t&t<2)+(-t+3).*(t>2&t<3);
p=@(z)integral(g,0,z);
figure;
fplot(p,[0 3]);hold on;