        set(0, 'DefaultAxesLineWidth', 2.0)
        set(0, 'DefaultLineLineWidth', 1.0)
        set(0, 'DefaultTextFontSize', 15)
        %set(0, 'DefaultTextFontWeight', 'bold')
        set(0, 'DefaultAxesFontSize', 15)
        %set(0, 'DefaultAxesFontWeight', 'bold')
        %set(0, 'DefaultLineMarkerSize', 10)

x=avsread('../huber/hc/fbp.fld');
image(256-x'/0.13*256); axis image; axis off; drawnow; print fbphc.ps

ideal=avsread('../huber/phtmtbl.fld');
x1=avsread('../huber/hc/out11.fld');
image(256-x1(:,:,11)'/0.13*256); axis image; axis off; drawnow;
print hb_img_hc.ps

load ../data4/myout4; x2=myout4;
image(256-x2(:,:,11)'/0.13*256); axis image; axis off; drawnow;
print my_img_hc.ps

x=avsread('../huber/phtmtbl.fld');
p1=plot(x(:,65)); hold on; p2=plot(x1(:,65,11),'o');
p3=plot(x2(:,65,11),'*'); hold off; axis([0 128 0 0.13])
xlabel('pixel index'); ylabel('estimate of \mu');
legend([p3 p2 p1],'Proposed penalty','Huber penalty','Ideal'); drawnow;
print profilecmphc.ps

ex=avsread('fbphcfbp.fld'); image(256-ex'/4e4*256);
axis image; axis off; drawnow;
print emfbp.ps

ex=avsread('fbphcmine3.fld'); image(256-ex(:,:,11)'/4e4*256);
axis image; axis off; drawnow;
print emmine.ps

ex=avsread('fbphchuber.fld'); image(256-ex(:,:,11)'/4e4*256);
axis image; axis off; drawnow;
print emhuber.ps

% ****************************************************************
x=avsread('../huber/lc/fbp.fld');
image(256-x'/0.13*256); axis image; axis off; drawnow; print fbplc.ps

ideal=avsread('../huber/phtmtbl.fld');
x1=avsread('../huber/lc/out11.fld');
image(256-x1(:,:,33)'/0.13*256); axis image; axis off; drawnow;
print hb_img_lc.ps

load ../myout4; x2=myout4;
image(256-x2(:,:,33)'/0.13*256); axis image; axis off; drawnow;
print my_img_lc.ps

x=avsread('../huber/phtmtbl.fld');
p1=plot(x(:,65)); hold on; p2=plot(x1(:,65,33),'o');
p3=plot(x2(:,65,33),'*'); hold off; axis([0 128 0 0.13])
legend([p3 p2 p1],'Proposed penalty','Huber penalty','Ideal'); drawnow;
xlabel('pixel index'); ylabel('estimate of \mu');
print profilecmplc.ps

ex=avsread('fbplcfbp.fld'); image(256-ex'/4e4*256);
axis image; axis off; drawnow;
print emfbplc.ps

ex=avsread('fbplcmine3.fld'); image(256-ex(:,:,33)'/4e4*256);
axis image; axis off; drawnow;
print emminelc.ps

ex=avsread('fbplchuber.fld'); image(256-ex(:,:,33)'/4e4*256);
axis image; axis off; drawnow;
print emhuberlc.ps

