въпрос за 1D_FDTD

[learnerlcy]

Здравейте ~ пускам на 1D_FDTD код на книгата, наречена ЕЛЕКТРОМАГНИТНИ СИМУЛАЦИЯ UING НА FDTD НАЧИН.Но в резултат на изчисление не е наред ~
[URL = http://images.elektroda.net/13_1258255028.jpg] http://images.elektroda.net/13_1258255028_thumb.jpg [/

Искам да знам защо?

Това е кода:
/ * 1D FDTD симулация в свободното пространство * /
# включват <math.h>
# включват <stdlib.h>
# включват <stdio.h>

# определят KE 200

Основните ()
(
поплавък предишен [KE], HY [KE];
Int N, K, KC, KE, NSTEPS;
поплавък T;
поплавък t0, разпространение, пулс;
FILE * РП / /, * fopen ();

/ * * Initialize /
за (К = 1; к <KE; K )
(
ех [к] = 0;
HY [к] = 0;
)

KC = KE / 2 / / център на проблема пространство
t0 = 40,0; / / Център за инцидента пулс
разпространение = 12; / / широчина на инцидента пулс
Т = 0;
NSTEPS = 1;

докато (NSTEPS> 0)
(
printf ( "NSTEPS -->"); / / NSTEPS е броя на основните линия е изпълнил
scanf ( "% D", & NSTEPS);
printf ( "% г \ N", NSTEPS);
N = 0;
за (N = 1; N <= NSTEPS; N )
(
T = T 1;
/ / Майн FDTD Loop
/ * Изчислете областта Ex * /
за (К = 1; к <KE; K )
(
ех [к] = предишен [к] 0,5 * (HY [K-1] - HY [к]);
)
/ / Сложете Gaussian импулс в средата
импулс = Годен (-5 * (Pow ((t0-T) / спред, 2,0)));
ех [KC] = пулс;
printf ( "% 5.1f% 6.2f \ N", t0 - Т, ех [KC]);
/ / Изчислете областта Хъ
за (к = 0; к <KE-1; K )
(
HY [к] = HY [к] 0,5 * (предишен [к] - ех [к 1]);
)
)
/ / Край на главния линия
/ / В края на изчислението, разпечатва Ex и Хъ полета
за (К = 1; к <= KE; K )
(
printf ( "% 3d% 6f% 6.2f \ N", к, ех [к], HY [к]);
)

/ / Напиши Е поле, за да файл "ех"
КП = fopen ( "Ex.dat", "W");
за (К = 1; к <= KE; K )
(
fprintf (РП, "% D% 6.2f \ N", к, ех [к]);
)
fclose (РП);

/ / Напиши областта H, за да файл "Хъ"
КП = fopen ( "Hy.dat", "W");
за (К = 1; к <= KE; K )
(
fprintf (РП, "% D% 6.2f \ N", к, HY [к]);
)
fclose (РП);
printf ( "T =% 5.0f \ N", T);
)

)Added след 1 минути:[URL = http://images.elektroda.net/53_1258255246.jpg] [IMG] http://images.elektroda.net/53_1258255246_thumb.jpg [/