 = m0 + dv; 
figure(1);imagesc(x,z,m);caxis([2250,2750]);xlabel('x [m]');ylabel('z [m]');title('true model [m/s]');colorbar;
figure(2);imagesc(x,z,m0);caxis([2250,2750]);xlabel('x [m]');ylabel('z [m]');title('initial model [m/s]');colorbar;

%%

% reconstructions model

% reflection example, source and receivers are on the top of the model
load FWIresult_camenbert.mat
figure(3);imagesc(x,z,results(:,:,end));caxis([2250,2750]);xlabel('x [m]');ylabel('z [m]');title('inverted result from relection data');color