amp+=sin(tav*k);

Biztos sokaknak szemetszurt ez a sor. Mert bizony nem igy szokas hullamokkal szamolni,

http://en.wikipedia.org/wiki/Wave

Nalam az x a tav valtozo,  Igenam, de ez az egyenlet idoben valtozo hullamot ad, mig az enyem egy 2 dimenzios megfagyott hullam statikus allapotat szamolja. Mind irtam, a terido hullam nem mozog. Most csak egy terdimenzio van, az x, amiatt az egesz egy 2 dimenzios fagyott allapotu hullamra hasonlit. Ezert eleg a tavolsaggal {x} szamolni, csak elobb a terbeli hullamhossz teridobelive kell atszamolni, ahogy a programban is teszem,

 l1*=sin(a1);//4d terido hullamhossz!
 l2*=sin(a2);

De meg lehet probalni atirni, kihagyva a fenti atszamitast, es kiszamolni az elektron frekvenciajat. A sin(kx-wt) egyenlet ugyan ezt  fogja kirajzolni,