Added error norm calculation function and example usage in run_boundary_layer.m