Метод непосредственной развёртки - C++ - Ответ 4677422

double f(double x)
{
return (2.1-x)*(2.6-x)*(3.1-x)+(2*1.1*1.1)-(3.1-x)-1.1*1.1*(2.1-x)-1.1*1.1*(2.6-x); // ф-я, корни которой мы ищем
}

// a, b - пределы хорды, eps - необходимая погрешность
double hord(double a,double b,double eps) // нахождение корней методом хорд
{
while(fabs(b - a) > eps)
{
a = b - (b - a) * f(b)/(f(b) - f(a));
b = a - (a - b) * f(a)/(f(a) - f(b));
}

// a - i-1, b - i-тый члены

return b;
}

void sis(double r1,double r2,double r3,double l) // метод итераций решение системы
{
double r0=15;
while (fabs(r0-r1) > 0.01)
{
r0=r1;
r1=(r2+1.1*r3)/(-2.1+l);
r2=(r1+1.1*r3)/(-2.6+l);
r3=(1.1*r1+1.1*r2)/(-3.1+l);
}
//return r1," ",r2," "r3;
}



double razvert()
{
double l[3],r1,r2,r3;
l[0]=hord(1,1.5,0.001);
l[1]=hord(1.5,2,0.001);
l[2]=hord(4.5,5,0.001);
for(int i=0;i<3;i++) cout<<l[i]<<endl;// характеристические числа находит правильно, не могу понять запутался с функцией sis которая решает уравнение и находит собственные вектора
r1=-0.5;r2=1.5;r3=1;
sis(r1,r2,r3,l[0]);
cout<<r1<<" "<<r2<<" "<<r3<<endl;



}

Добавлено через 1 минуту
пожалуйста помогите исправить программу, если я правильно понимаю то ошибка моя с return' ом но чё то никак не могу догнать как сделать правилно

Вернуться к обсуждению:
Метод непосредственной развёртки C++