function [chisquare] = fit(probs,datahist,stdevs) % This function takes three vector arguments: Gaussian PDF predictions for % each number of events, Poisson PDF predictions, and the data histogram % for each number of events stdevs2 = (stdevs.^2)'*ones(1,length(probs)); chisquare = sum((((probs-datahist).^2)./stdevs2)'); % chisquare = sum(((probs'-datahist).^2)./probs'); % chisquarepoiss = sum(((poissprobs'-dathist).^2)./poissprobs'); return