фіз.-мат. наук А.Н.Полюдова і В.А.Панченко.- Москва «Мир». 1990 - 350 с.
. Д.В.Хеерман. Методи комп ютерного експерименту в теоретичній фізиці. Переклад з англ. В.Н.Задкова - Москва «Наука». Головна редакція фізико-математичної літератури. 1990 - 174 с.
. К.Жаблон, Ж.-К. Симонов. Застосування ЕОМ для чисельного моделювання у фізиці. / Переклад з франц. А.В.Арсентьевой - Москва «Наука». Головна редакція фізико-математичної літератури. 1983 - 234 с.
Додаток
# include
# include
# include
# include
# include
# include
# include metod (double x [], double y [], double vx [], double vy [], double ax [], double ay [], int N, int i); generator (int max)
{
return random (max);
} setka (double mx, double my, double k)
{
char str [4];
int i;
setcolor (11);
for (i=0; i <= 30; i +=5)
{
line (347, i * my +5 + k, 353, i * my +5 + k);
}
for (i=0; i <= 300; i +=50)
{
itoa (i, str, 10);
outtextxy (i * mx +347,147 + k, str);
}
setcolor (1);
line (320,140 + k, 640,140 + k);
line (350,0 + k, 350,160 + k);
} metod (double x [], double y [], double vx [], double vy [], double ax [], double ay [], int N, int i)
{
double axn [20], ayn [20], r [20], vxn [20], vyn [20], xn [20], yn [20];
double V, t=0.01;
axn [i]=0;
ayn [i]=0;
for (unsigned j=0; j
{
if (i!=j)
{
r [i]=sqrt (pow ((x [i]-x [j]), 2) + pow ((y [i]-y [j]), 2));
V=24 * ((1/pow (r [i], 7)) - (2/pow (r [i], 13)));
axn [i] +=(-V) * (x [i]-x [j]) / r [i];
ayn [i] +=(-V) * (y [i]-y [j]) / r [i];
}
}[i]=0.5*(axn[i]+ax[i])*t;[i]=0.5*(ayn[i]+ay[i])*t;[i]+=vx[i]*t+0.5*ax[i]*t*t;[i]+=vy[i]*t+0.5*ay[i]*t*t;[i]=axn[i];[i]=ayn[i];
} main ()
{
clrscr ();
int driver=DETECT, mode, e;
initgraph (& driver, & mode, «c: BorlandC BGI»);
e=graphresult ();
if (e!=grOk)
{
printf ("% s», grapherrormsg (e));
}
else
{
double vx [20], y [20], x [20], vy [20], ax [20], ay [20], T=0, t=0.01;
int i=0, N=...