Naime, treba numerički rešiti vremenski zavisnu jednačinu protoka snage, koja ima sledeći oblik:

gde je: A-faktor multiplikacije drugog reda u izrazu za koeficijent gubitka snage
D-koeficijent sprezanja
z-dužina vlakna
θ-ugao prostiranja svetlosti u odnosu na osu vlakna
Granični uslovi su:
Nakon sređivanja imamo:

nakon primene Furijeove transformacije, i prelaska u frekventni domen jednačina postaje:

Kako je


i

Primenom eksplicitnog metoda konačnih razlika, koristeći:
1) šemu centralne razlike za izvode



2) šemu prednje razlike za izvod


Na taj način dobijamo dve uparene jednačine:

i

pri čemu granični uslovi postaju:

Da bismo izbegli problem singularnosti u tačkama θ=0 koristimo:

Ovo se dalje rešava numerički, pri čemu je upadni snop svetlosti, u vremenskom domenu, u obliku ravanskog talasa dat pomoću Dirakovog impulsa u vremenu, pod nekim uglom

Mene, kao prvo zanima kako će ulazni impuls izgledati u ferkventnom domenu, i kako ga implementirati na početku vlakna prilikom programiranja.
Drugo što me zanima je kako da dobijem frekventni odziv, nakon što sam dobio

