% Example 3.2
time = [ -0.2400 -0.0600  0.1200  0.3000  0.4800  0.6600  0.8400 ...
          1.0200  1.2000  1.3800  1.5600  1.7400 ]';
Y = [  0.0000  0.0000  0.0034  0.0337  0.0964  0.1907  0.2790 ...
       0.3351  0.3588  0.3600  0.3600  0.3600 ]';
dt = time(2) - time(1);
time_fine = [ -.24:0.006:1.8 ]';
qhat = Stolz( dt, x, Y );
figure(96); clf; 
plot( time - dt/2, qhat, 'ro', time_fine, q_triangle(time_fine), ... 
      'k--', 'linewidth', 1.0)

xlabel('time')
ylabel('Heat Flux')

% fig=figure(96);
% print(fig,fullfile(pwd, 'Figure3_3.tif'),'-dtiff','-r600')