function [xp] = x_expectation_value(op,data)

xp = sum(sum( op .* abs(data).^2 )) / sum(sum( abs(data).^2 ));
