 imagesc(zz,xx,v_t);xlabel('X [m]');ylabel('Z [m]');caxis(a);colorbar
  title('Result w/ the true source wavelet')

  load('../results/Result_estsrc.mat');
  m_e = Result.mf;
  v_e = reshape(1./sqrt(m_e),n);
  info_e = Result.info;
  info_e = info_e.obj;

  subplot(1,2,2)
  imagesc(zz,xx,v_e);xlabel('X [m]');ylabel('Z [m]');caxis(a);colorbar
  title('Result w/ source estimation')
  x0=10;
  y0=10;
  width=900;
  height=400;
  set(gcf,'units','points','position',[x0,y0,width,height])

%The comparison of the 