% Example 3.1
time = [ -.24:0.06:1.8 ]';
x = 0.5;
T_Duhamel =Duhamel_q( x, time, 'q_triangle');
T_exact = 0.6*T_1pulse( time, 0, 1.2, x )'; 
fprintf(' % 6.4f  %10.8f  %10.8f\n',[ time T_exact T_Duhamel ]'  )
figure(96); clf; plot( time, T_exact, 'k--', ...
                        time, T_Duhamel, 'ro', 'linewidth', 1.0 )
xlabel('time')
ylabel('Temperature')
    
fig=figure(96);
print(fig,fullfile(pwd, 'Figure3_2.tif'),'-dtiff','-r600')
