pro calc_incoherent,ki,kf,Qana,Qmon


    kix = ki[0]
    kiy = ki[1]
    kfx = kf[0]
    kfy = kf[1]

    kixincoherent = sqrt(kf[0]^2 + kf[1]^2)*ki[0]/sqrt(ki[0]^2 + ki[1]^2)
    kiyincoherent = sqrt(kf[0]^2 + kf[1]^2)*ki[1]/sqrt(ki[0]^2 + ki[1]^2)

    kfxincoherent = sqrt(ki[0]^2 + ki[1]^2)*kf[0]/sqrt(kf[0]^2 + kf[1]^2)
    kfyincoherent = sqrt(ki[0]^2 + ki[1]^2)*kf[1]/sqrt(kf[0]^2 + kf[1]^2)

    Qxanaincoherent = kix - kfxincoherent
    Qyanaincoherent = kiy - kfyincoherent
    Qxmonoincoherent = kixincoherent - kfx
    Qymonoincoherent = kiyincoherent - kfy

    Qana = [Qxanaincoherent,Qyanaincoherent]
    Qmon = [Qxmonoincoherent,Qymonoincoherent]


end;calc_incoherent
