Теми рефератів
> Реферати > Курсові роботи > Звіти з практики > Курсові проекти > Питання та відповіді > Ессе > Доклади > Учбові матеріали > Контрольні роботи > Методички > Лекції > Твори > Підручники > Статті Контакти
Реферати, твори, дипломи, практика » Новые рефераты » Метод молекулярної динаміки

Реферат Метод молекулярної динаміки





фіз.-мат. наук А.Н.Полюдова і В.А.Панченко.- Москва «Мир». 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=...


Назад | сторінка 5 з 7 | Наступна сторінка





Схожі реферати:

  • Реферат на тему: Double Standards in Modern Politics
  • Реферат на тему: Double Entry Types of Balance Sheet
  • Реферат на тему: Double electric layer. Mechanism of formation and theory of structure
  • Реферат на тему: Застосування систем комп'ютерного моделювання для дослідження математич ...
  • Реферат на тему: Переклад англійської Наукової и технічної літератури