%mm2103.mx = (0:.1:1)*2*pi;y = sin(x); % create rough datapp = spline(x,y); % pp-form fitting rough datappi = mmppint(pp,0); % pp-form of integralxi = linspace(0,2*pi); % finer points for interpolationyi = ppval(pp,xi); % evaluate curveyyi = ppval(ppi,xi); % evaluate integralplot(x,y,'o',xi,yi,xi,yyi,'--') % plot resultstitle('Figure 21.3: Spline Integration')