    % Experimental Setup
    % clear all;

    init;

    guideExpSetupTrans; % set transmitters positions
    guideExpSetupRece;  % set receivers positions
    guideSetGrid; % set a grid

    seti.model = 'helmholtz2D';
    seti.k = 250; % wave number

    seti = setIncField(seti); % Evaluate incident fields on ROI
    seti = setMeasKer(seti); % set up measurement kernels (such that measurement = k^2 * kernel * solution * voxelVolume)

    figure(1); imagesc(real(reshape(seti.incField(:,1),[seti.nROI seti.nROI]))); axis xy; colorbar; % incident field in ROI of 1st transmitter
    % set(gca,'FontSize',20); axis square; print(1,'-depsc','guideExpSetup1.eps');

    figure(2); imagesc(real(reshape(seti.measKer(5,:),[seti.nROI seti.nROI]))); axis xy; colorbar; % measurement kernel on ROI of 5th receiver
    % set(gca,'FontSize',20); axis square; print(2,'-depsc','guideExpSetup2.eps');