% load VMP
vmp = xff('*.vmp');

% compute conjunction map of first two maps, default options
vmp.ComputeFormula('conjvalp(#1, #2)', struct('pvalues', true));

% save VMP
vmp.Save;