[SOLVED] what is error of this 3D FDTD matlab code?

[ghasem_008]

hi friends.
I wrote a 3D FDTD matlab code in cartesian coordinates with first order Mur ABC.
but I don't konw why my Fields become very big (in order of 1e20!!!).
can anyone find it's reason and find my error?
I don't konw how import a current source to bottom surface of my antenna.
my antenna is dielectric with high negative permitivity such as -100).
My Courant limit satisfied and mesh is regular.
I attach my code.please help me.
where I mistake?
View attachment my code.txt

 

[sara67]

hi ghasem
for a courent source you should study this thesis
it's useful