function dipoleMoment = obtainDipoleMoment(dminpu,wavefunc)
%
position = dminpu(:,1);
dm_temp = dminpu(:,2);
dipoleMoment = zeros(size(wavefunc,2),1);
for ii = 1:size(wavefunc,2)
dipoleMoment(ii) = trapz(position, ...
conj(wavefunc(:,ii)).*dm_temp.*wavefunc(:,ii));
end
end