round(0.9 / dt) + 1, :);
u_slice_2 = u_slice_1(:);
figure; hold on;
title(sprintf('t = %.2f', 0.9));
plot(x_, u_slice_1, 'r-', x_, u_slice_2, 'b-');
legend('dt = 0.001', 'dt = 0.1');
xlabel('x'); ylabel('u');
fprintf('max error = %e\n', max(abs(u_slice_1  - u_slice_2)));
##### SOURCE END #####
--></body></html>