le('True model')

  subplot(1,2,2)
  imagesc(zz,xx,velini);xlabel('X [m]');ylabel('Z [m]');caxis(a);colorbar
  title('Initial model')
  x0=10;
  y0=10;
  width=900;
  height=400;
  set(gcf,'units','points','position',[x0,y0,width,height])

%% Inversion result

%The inversion results w/ the true source wavelet and w/ the source estimation are:

  load('../results/Result_truesrc.mat');
  m_t = Result.mf;
  v_t = reshape(1./sqrt(m_t),n);
  info_t = Result.info;
  info_t = info_t.obj;

  figure;subplot(1,2,1)
 