Читать курсовая по математике: "Задача Діріхле" Страница 6


назад (Назад)скачать (Cкачать работу)

Функция "чтения" служит для ознакомления с работой. Разметка, таблицы и картинки документа могут отображаться неверно или не в полном объёме!

Nach_Prav(y);i=2:ny; for j=2:nx;(i,j)=(u(i,1)+u(i,nx+1)+u(1,j)+u(ny+1,j))/4;

//Метод Лібмана=1d>1e-5 a = u(ny, nx) for i = 2:ny for k = 2:nx(i, k) = (u(i - 1, k) + u(i + 1, k) + u(i, k - 1) +(i, k + 1)) / 4; end= abs(a - u(ny, nx)); end

Після закінчення розрахунків розв’язок - вектори у , х, матрицю u можна вивести на екран Результати наведені у табл. 1.

Таблиця 1 - Результат розв’язання прикладу

х=0

х=1

х=2

х=3

х=4

у=0

0

0

0

0

0

у=і

0

0,0405074

0,0898696

0,1331149

0,1111111

у=2

0

0,0722196

0,1858920

0,3314878

0,4444444

у=3

0

0,0625

0,25

0,5625

1

Для більш точного розв’язання можна розділити область на більш дрібну сітку. Для наочного представлення результатів доцільно користуватися побудовою трьохвимірного графіка за допомогою комп’ютерних засобів.

2. ТЕКСТ ПРОГРАМИ # include < stdio.h >

# include < fstream.h >

# include < conio.h >

# include < iostream.h >

# include < math.h >i, j, k;// Variablesh, x, y, tmp, E1;point {xx;yy;BelongsToDh_;BelongsToDh;F;F_;

}, arrayP [13] [33];arrayX [13];arrayY [33];diff [500];CreateNet (void);// Procedure PrototypesIsLineFit (float Param);CrMtrD (void);RegArrayX ();RegArrayY ();CreateDh_ ();IsFit (point Par);FillF ();CreateDh ();IsInner (int i, int j);FillF_ ();CountDif ();MakeFile ();main (void)// MAIN

{();.xx = 3;.yy = 5;= 0.2;.BelongsToDh_ = 1;.BelongsToDh = 1;();();();();_ ();();();_ ();();(E1> = 0.005) {(i = 0; i



Интересная статья: Основы написания курсовой работы