a = imread( 'cameraman.tif' );
    a = im2double(a);
    [L, Lx, Ly, Lxx, Lxy, Lyy] = g2Jet( a, 3 );
    subplot( 3,3,1 ); imshow(L, []);   title('L');
    subplot( 3,3,4 ); imshow(Lx, []);  title('Lx');
    subplot( 3,3,5 ); imshow(Ly, []);  title('Ly');
    subplot( 3,3,7 ); imshow(Lxx, []); title('Lxx');
    subplot( 3,3,8 ); imshow(Lxy, []); title('Lxy');
    subplot( 3,3,9 ); imshow(Lyy, []); title('Lyy');