Решение интеграла методом трапеции

[Тёма777]

помогите написать программы для вычесления интеграла методом трапеции :
Интеграл от sin*mx*cos*nx*dx

[_ILYA_]

сколько платишь?

[Гром]

Сделаю за деньги. Аська в профиле.

[Stilet]

Интеграл от sin*mx*cos*nx*dx

Что-то я не пойму че за интрегал такой...
Можешь нарисовать формулу как полагается математически?

[Гром]

Stilet
Благодаря

#include <telepat.h>

получается
[sin(m*x) * sin (n*x)] dx

[davinci]

Могу кодик подбросить, у меня лаба в универе такая была, правда интеграл(формулу) туда свой сам положишь

[Stilet]

#include <telepat.h>

Нету. Пропил я его.
Ждем от автора пояснений...

[davinci]

Такая программка, которую я делал для лабораторной работы

#include <iostream>
#include <conio>
using namespace std;

int main()
{
	float a,b,n,h;//a и b - границы интегрирования, n - количество итераций
	cout<<"Enter a,b and n"<<endl;
	cin>>a>>b>>n;
	h=(b-a)/n;//шаг итерации
	float I;
	I=(exp(sin(a))*sin(a)-exp(sin(b))*sin(b))/2;//Ваш интеграл от a и от b
	int i;
	for (i = 0; i<=n; i++)
		I+=exp(sin(a+i*h))*sin(a+i*h);//тоже Ваш интеграл
	cout<<I*h;
	getch();
	return 0;
}

Прошу прощения за корявость, но программа работает...
Если надо, есть еще и схожие программки для метода срединных прямоугольников и метода парабол (метод Симпсона)

[dxdy]

внесу свои 5 копеек...

#include <iostream>
#include <math.h>
#include <locale.h>
using namespace std;
const float p = 3.141592653;

float ingdv(float a, float b, float e, float (*f)(float,float,int));
float f1(float a, float b, int n);
float f2(float a, float b, int n);
float f3(float a, float b, int n);

int main()
{
	setlocale(LC_CTYPE,"Russian");
	float a,b,rez,e;
	printf("Введите A, B и точность вычисления\n");
	scanf("%f%f%f", &a,&b,&e);
	printf("\nРешение интеграла 3-мя методами:\n\n");
	rez = ingdv(a,p/b,e,f1);
	printf("Метод Трацеция %f\n",rez);
	rez = ingdv(a,p/b,e,f2);
	printf("Метод Прямоугольник %f\n",rez);
	rez = ingdv(a,p/b,e,f3);
	printf("Метод Симпсона %f\n\n",rez);
	return 0;
}
/*-----------------Метод трапеции-------------------*/
float f1(float a, float b, int n)
{
	float h,sum,x;
	int i;
	sum = 0;
	h=(b-a)/n;
	sum += sqrt(tan(a));
	for(i = 1; i < n; i++)
	{
		x = a + i*h;
		sum += 2*sqrt(tan(x));
	}
	x = a + n*h;
	sum += sqrt(tan(x));
	return(0.5*h*sum);
}
/*---------------Метод прямоугольника------------------*/
float f2(float a, float b, int n)
{
	float h,sum,x;
	int i;
	sum = 0;
	h=(b-a)/n;
	for(i = 1; i <= n; i++)
	{
		x = a + (i-0.5)*h;
		sum += sqrt(tan(x));
	}
	return h*sum;
}
/*----------------Метод Симпсона--------------------*/
float f3(float a, float b, int n)
{
	float h, sum,x;
	int i, k;
	sum = 0;
	h=(b-a)/n;
	sum += sqrt(tan(a));
	for(i = 1; i < n; i++)
	{
		x = a + i*h;
		if(i%2 == 0) {k = 2;}
		else {k = 4;}
		sum += k*sqrt(tan(x));
	}
	x = a + n*h;
	sum += sqrt(tan(x));
	return 0.333*h*sum;
}
/*-----------Общая функция-----------*/
float ingdv(float a, float b, float e, float (*f)(float,float,int))
{
	int n = 10;
	float p,k;
	k = f(a,b,n);
	do{
		p = k;
		k = f(a,b,2*n);
		n = n*2;

	}while(abs(p-k)<e);
	return k;
}

[Тёма777]

ребят, спасибо, проверю код отпишусь !!!