%scattering of the photon
function [photon_]=photon_scatter(photon_)
photon_.theta=acos(1-rand()*(1-cos(pi)));
photon_.fai=2*pi*rand();
multimatrix_=reference_rotate(photon_.theta,photon_.fai);
photon_.matrix=multimatrix_*photon_.matrix;
photon_.position=photon_.matrix*photon_.positioninoriginal;
end
