h = 2.^[ -1 : -1 : -6 ];
err = [ 1.147e-1, 2.840e-2, 7.084e-3, 1.770e-3, 4.425e-4, 1.106e-4 ];
loglog( h, err, '-ok', h, h, '-.b', h, h.^2, '--k', h, h.^3, ':r' );
xlabel( 'h' ); ylabel( 'err' ); grid on
legend( 'err', 'h', 'h^2', 'h^3', 'Location', 'SouthEast' );