Системи нелінійних рівнянь (стр. 5 из 5)

int i;

for(int i=n-2;i>=1;i--)

{

k=0;

for(j=1;j<=i;j++)

{

p[m]=(p[m-i]-p[m-i-1])/(x[k+l]-x[k]);

// cout<<p[m]<<"&bsol;t";

m++;

k++;

}

// cout<<"&bsol;n";

l++;

}

j=0;

for(i=1;i<n;i++)

{

y[i]=p[j];

j=j+(n-i);

}

}

void main()

{

int const n=5,N=(n-1)*n/2;

double x[]={1,2,3,4,5};

double y[n];

for(int i=0; i<n; i++) y[i]=1+2*x[i]+3*pow(x[i],2)+4*pow(x[i],3)+5*pow(x[i],4);

double *p=new double[N];

Nuton1(x,y,n,p);

//for (i=0;i<n;i++) cout<<y[i]<<'&bsol;t';

double z=1.5;

cout<<Nuton(y,n,x,z)<<'&bsol;t'<<1+2*z+3*pow(z,2)+4*pow(z,3)+5*pow(z,4)<<'&bsol;n';

getch();

}

Вікно виконання програми:

Висновок

У даній курсовій роботі я описала методи розв’язання систем нелінійних рівнянь, а також навела приклади їх реалізації. В пакеті Mathcad для розв’язання систем нелінійних рівнянь я використала оператор solve, обчислювальний блок given – find, метод простих ітерацій, метод Ньютона та два модифіковані методи Ньютона, в яких на відміну від метода Ньютона матриця Якобі знаходилась не на кожному кроці, а тільки на початковому.

В пакеті Excel розв’язати систему нелінійних рівнянь можна за допомогою засобу «Пошук рішення», в якому замість системи рівнянь необхідно розв’язати рівносильне їй рівняння. Також моє знання мови програмування С++ дозволило мені продемонструвати приклади розв’язання систем нелінійних рівнянь методом ітерацій та методом Ньютона саме мовою програмування С++.

Список використаної літератури

1. Ю.Ю Тарасевич «Численные методы на Mathcad’e»

2. Д. Гурский, Е. Турбина «Вычисления в Mathcad»

3. Б. Демидович «Основы вычеслительной математики»

4. О. Кваша «Численные методы»

5. В.Н. Тарасов, Н.Ф. Бахарева «Численные методы»