Аппроксимация.Линейная интерполяция

@slevin_45 · Регистрация: 17.12.2013

Здравствуйте!Помогите,пожалуйста,найти ошибку в коде.Высчитывает неверные корни.

C++

#include<iostream>
#include<math.h>
using namespace std;
int main (void)
{
    int i,j,m=12;
    double x_t,x_m[12]={0},y[22]={0},y_m[12]={0},x[22]={0},polinom,pogr,n=2,a_v[12]={0},b_v[12]={0},a=-4,b=2;
    for(i=1;i<m;i++)
    {
        x_m[i]=a+(i-1)*(b-a)/(m-1);
    }
    for(i=1;i<m;i++)
    {
        y_m[i]=0.1*pow(x_m[i],3)+pow(x_m[i],2)-10*sin(x_m[i]);
    }
    for(j=1;j<22;j++)
    {
        x[j]=a+(j-1)*(b-a)/20;
    }
    for(j=1;j<22;j++)
    {
        y[j]=0.1*pow(x[j],3)+pow(x[j],2)-10*sin(x[j]);
    }
    for(j=1;j<22;j++)
    {
        x_t=x[j];
        i=0;
stop:
        i++;
        if(x_t<=x_m[i]&&x_t>=x[i-1])
        {
            polinom=a_v[i]*x_t+b_v[i];
            a_v[i]=(y_m[i]-y_m[i-1])/(x_m[i]-x_m[i-1]);
            b_v[i]=y_m[i-1]-a_v[i]*x_m[i-1];
            cout<<endl;
            pogr=y[j]-polinom;
            cout<<pogr;
        }
        else
        {
            goto stop;
        }
    }
    system("pause");
    return(0);
}

@S_el · 13.03.2014, 10:39

Сообщение от slevin_45

Высчитывает неверные корни.

О каких корнях речь?

@slevin_45 · 13.03.2014, 18:23 **[ТС]**

Ошибка найдёна.

@Ev[G]eN · 13.03.2014, 18:29

slevin_45, будьте добры напишите в чем она заключалась и как была исправлена. Возможно, это пригодится другим пользователям. Заранее спасибо.

@S_el 2431 / 1831 / 404 Регистрация: 15.12.2013 Сообщений: 8,124
	13.03.2014, 10:39	2
	Сообщение от slevin_45 Высчитывает неверные корни. О каких корнях речь? 1

@slevin_45 1 / 1 / 0 Регистрация: 17.12.2013 Сообщений: 66
	13.03.2014, 18:23 [ТС]	3
	Ошибка найдёна. 0

@Ev[G]eN Mobile Software Engineer 5365 / 1584 / 956 Регистрация: 23.01.2011 Сообщений: 3,214
	13.03.2014, 18:29	4
	slevin_45, будьте добры напишите в чем она заключалась и как была исправлена. Возможно, это пригодится другим пользователям. Заранее спасибо. 0

Опции темы