(* analiticen rezultat za izbiro chi1=chi2=chi/2 *)

fun[chi_,fi_]:=(1-Cos[chi]+Cos[fi]-Cos[fi]Cos[chi])/(1+Cos[fi/2]^4-2 Cos[chi]Cos[fi/2]^2)

(* numericen izracun zacensi z vmesnimi analiticnimi rezultati *)

a[chi1_,chi2_,fi1_,fi2_]:=1/4 (Exp[2I chi1]+Exp[I (chi1+chi2)]Exp[-I (fi2-fi1)])
b[chi1_,chi2_,fi1_,fi2_]:=1/4 (Exp[2I chi2]+Exp[I (chi1+chi2)]Exp[+I (fi2-fi1)])

tgmintd[chi1_,chi2_,fi1_,fi2_]:= (Exp[I chi1]Exp[I fi1]-Exp[I chi2] Exp[I fi2])/(1-a[chi1,chi2,fi1,fi2]-b[chi1,chi2,fi1,fi2])

tg[chi1_,chi2_,fi1_,fi2_]:=Exp[I chi1] Exp[I fi1]+a[chi1,chi2,fi1,fi2] tgmintd[chi1,chi2,fi1,fi2]

td[chi1_,chi2_,fi1_,fi2_]:=Exp[I chi2] Exp[I fi2]-b[chi1,chi2,fi1,fi2] tgmintd[chi1,chi2,fi1,fi2]

t[chi1_,chi2_,fi1_,fi2_]:=(tg[chi1,chi2,fi1,fi2]+td[chi1,chi2,fi1,fi2])/2

(* Spodaj se vidi, da sta za chi1=chi2=chi/2 oba rezultata enaka*)

Plot[{Abs[t[1.3,1.3,0,fi Pi]]^2,fun[2.6,fi Pi]},{fi,0,4}]

tt[chi1_,chi2_,fi1_,fi2_]:= Abs[t[chi1,chi2,fi1,fi2]]^2

tt[chi1_,chi2_,fi2_]:= Abs[t[chi1,chi2,0,fi2]]^2

Plot[{tt[1.3,1.3,Pi fi],tt[1.6,1.0,Pi fi],tt[1.9,0.7,Pi fi],tt[2.2,0.4,Pi fi],tt[2.5,0.1,Pi fi],tt[2.8,-0.2,Pi fi],tt[3.1,-0.5,Pi fi]},{fi,0,4}]

Plot[{tt[1.3,1.3+1,Pi fi],tt[1.6,1.0+1,Pi fi],tt[1.9,0.7+1,Pi fi],tt[2.2,0.4+1,Pi fi],tt[2.5,0.1+1,Pi fi]},{fi,0,4}]

povprecje[fi_,chi2_,n_]:=(1/n)Sum[Abs[t[2 Pi i/n, chi2 ,0,fi]]^2,{i,1,n}]

povprecje2[fi_,n_]:=(1/n) Sum[povprecje[fi,2 Pi i/n,n],{i,1,n}]

Plot[povprecje2[fi Pi, 4],{fi,0,4}]