Вычисление интеграла методом Ньютона-Котеса

@Илья007 · Регистрация: 17.10.2015

Author24 — интернет-сервис помощи студентам

Дан интеграл: $https://www.cyberforum.ru/cgi-bin/latex.cgi?\int_{0}^{\pi /2}\frac{x+sin(x)}{1+cos(x)}$

Методом Ньютона-Котеса вычислить интеграл (n=6)
У меня такой код, но есть ошибки, не запускается, скажите пожалуйста как их устранить? Заранее благодарен.

C++

#include<iostream>
#include<math.h>
using namespace std;
typedef double D;
D nyut_k_6(D a,D b,D f(x),eps);// прототип
int main()
{ 
setlocale(LC_ALL,"rus_rus.1251"); 
int n=6;
D a,b,eps,f(x);
D H[6]={216,27,272,27,216,41};
D H0=41;
D INT1,INT2;
D h,S,hk,xi; 
int k=1; 
int i,m; 
cout<<"допустимая точность вычисления = "<<eps<<endl; 
do 
{ 
    nyut_k_6(0,90,fun,0.0001);// Обращение к функции
k=k*2; 
INT1=INT2; 
hk=(b-a)/k; // шаг интегрирования
h=hk/n;
S=(f(a)+(b))*H0;
m=n*k;
for(i=1;i<m;i++) 
{ 
    xi=a+i*h;
    S+=H[(i-1)%n]*f(xi);
} 
INT2=hk*S;
} 
while(fabs(INT1-INT2)>=eps) ;//модуль
return INT2; 
}
 
 
 
D nyut_k_6(D a,D b,D f(x),eps)// метод Ньютона-Котеса
{
 D fun(D x)
{
    
    int k_obr;
    k_obr++;// Кол-во обращений к функции
 
    return (x+sin(x))/(1+cos(x));
 }
return 1;
}

@Igrich · 07.03.2016, 17:38

C++

D nyut_k_6(D a,D b,D f(x),eps){// метод Ньютона-Котеса
 D fun(D x){
  int k_obr;
  k_obr++;// Кол-во обращений к функции
  return (x+sin(x))/(1+cos(x));
 }
 return 1;
}

Определяешь внутри функции функцию?(D fun(D x) )

@Илья007 · 08.03.2016, 07:14 **[ТС]**

ну а как исправить?

C++

#include<iostream>
#include<math.h>
using namespace std;
typedef double D;
D nyut_k_6(D a,D b,D f(D),D eps);// прототип
 
int main()
{ 
setlocale(LC_ALL,"rus_rus.1251"); 
int n=6;
D a,b,eps,f(D);
D H[6]={216,27,272,27,216,41};
D H0=41;
D INT1,INT2;
D h,S,hk,xi; 
int k=1; 
int i,m;
cout<<"допустимая точность вычисления = "<<eps<<endl; 
do 
{ 
    nyut_k_6(0,90,fun,0.0001);// Обращение к функции
k=k*2; 
INT1=INT2; 
hk=(b-a)/k; // шаг интегрирования
h=hk/n;
S=(f(a)+(b))*H0;
m=n*k;
for(i=1;i<m;i++) 
{ 
    xi=a+i*h;
    S+=H[(i-1)%n]*f(xi);
} 
INT2=hk*S;
} 
while(fabs(INT1-INT2)>=eps) ;//модуль
return INT2; 
}
 
 
D fun(D x)
{
    
    int k_obr=0;
    k_obr++;// Кол-во обращений к функции
 
    return (x+sin(x))/(1+cos(x));
 }
D nyut_k_6(D a,D b,D f(D),D eps)// метод Ньютона-Котеса
{
 
return 1;
}

@Igrich · 08.03.2016, 12:00

C++

1
2
3

D nyut_k_6(D a,D b,D f(D),D eps){
    return 1;
}

Мне чертовски нравится эта функция)
А вообще ну уж очень много тут ошибок,причем даже с самим синтаксисом. Серьезно,почитай учебник какой-нибудь,или одну из сотни статей,что бы не писать такие функцию.
ДА просто посмотри на прототип функции:

C++

1	D nyut_k_6(D a,D b,D f(x),eps);

1)зачем сокращать double до D? Как по мне лишнее;
2)D f(x) - в аргументе ты таким образом объявляешь еще одну функцию,к тому же с аргументом типа x, тоже нельзя так делать;
3)eps - что это такое? Он у тебя даже выше не объявлен.

Прочитай про объявление,переменные,типы,функции...и тогда ты сам сможешь все легко написать. Удачи тебе в изысканиях!

@Илья007 · 08.03.2016, 15:45 **[ТС]**

Сокращал D специально, так сказано по заданию. Сделал по новой, но что-то опять не так(((

C++

#include<iostream>
#include<math.h>
using namespace std;
typedef double D;
D fun(D x);
D nyut_k_6(D a,D b,D f(D),D eps) ;
void main()
{
setlocale(LC_ALL,"rus_rus.1251"); 
D a=0;
D b=1;
D eps;
D prib_nyut;
cout<<"Введите допустимую точность вычисления:  ";
cin>>eps;
prib_nyut=nyut_k_6(a,b,fun,eps);
cout<<"Ньютон-К. для n=6, Интеграл= "<<prib_nyut<<endl;
}
D nyut_k_6(D a,D b,D f(D),D eps) // метод Ньютона-Котеса
{
    setlocale(LC_ALL,"rus_rus.1251"); 
int n=6;
D H[6]={216,27,272,27,216,41};
D H0=41;
D INT1,INT2=1;
D h,S,hk,xi; 
int k=1; 
int i,m;
do 
{ 
k=k*2; 
INT1=INT2; 
hk=(b-a)/k; // шаг интегрирования
h=hk/n;
S=(f(a)+(b))*H0;
m=n*k;
for(i=1;i<m;i++) 
{ 
    xi=a+i*h;
    S+=H[(i-1)%n]*f(xi);
} 
INT2=hk*S;
}
while(fabs(INT1-INT2)>=eps) ;//модуль
return INT2;
}
D fun(D x)
{
    return (x+sin(x))/(1+cos(x));
 }

@Juffin · 08.03.2016, 15:56

C++

1	D nyut_k_6(D a,D b,D f(D),D eps)

Передавать функцию как параметр так, как это делаешь ты, нельзя. Это делается с помощью указателя на функцию:

C++

1	D nyut_k_6(D a,D b,D (*f)(D),D eps)

А вообще читай хотя бы, какие ошибки компилятор выдает.

	08.03.2016, 15:56

Опции темы