Численное интегрирование

[А_л_е_н_а]

Здравствуйте уважаемые мастера и гуру програмирования.

С языком СИ я лично познакомилfcm только сегодня, и поэтому понимаю пока не очень хорошо )

есть выражение

y=∫(снизу а сверху b)f(x)dx


нужна ваша помощь в написании программы в VisialStudio (ver. от 6 до 10)

решить этот определенный интеграл

честно говоря сама без понятия как это делать, вот сижу читаю умные книжки по СИ. Пока читаю о вводе переменных =)

если не очень заняты и можете помочь, буду вам очень признательна.

среда Visial Studio язык C консольное приложение

решается методом трапеции или Симпсона

трапеция - (f(a)+f(b))/2*(b-a)

f(x) принадлежит <math.h>
тип интерфейса <window.h>
<conio.h>
<stdio.h>


В большей степени интересует диалог пользователя с программой, т.е. вывод функций на экран, предположим это будут Sin x, Cos x, x*x. С вопросом "выберите функцию"
и подставить выбранную функцию в интеграл и решить методом трапеции или симпсона.
пределы интегралла задаем а шаг меняется программой для большей точности.(адаптивные алгоритмы численного интегрирования, когда на каждом шаге оценивается ошибка, и если ошибка слишком велика, то шаг интегрирования уменьшается (и наоборот, если ошибка достаточно мала, то шаг увеличивается). Это позволяет выполнить интегрирование за меньшее число шагов, но с такой же или даже меньшей ошибкой)


Написала ссуть задачи, как могла ))) Буду очень признательна есть поможете...Спасибо вам.

[Stilet]

Ну хорошо, допустим описать это не сложно, но ты сможешь по пунктам рассказать алгоритм математически?
Как что считается.

[А_л_е_н_а]

http://www.cleverstudents.ru/method_of_trapezoids.html вот тут написано как решать методом трапеции, т.е. математику я понимаю, а вот с кодом и синтаксисом СИ еще не совсем в ладах.
НАдеюсь на вашу помощь

[rlib]

#include <math.h>
#include <stdio.h>
double f(double x) {
   return pow(x,3)-5*cos(x);
}
int main() {
   int i;
   int np; // number of points
   double a, b; // limits
   double s; // integral
   double dx; //step
   double x; 
   a = 1.0; // from 1
   b = 2.0; // to 2
   np = 1000;
   dx = (b-a)/np; // step
   s = 0.0; //init sum
   x = a;
   for (i=0; i<np; ++i) {
	s += f(x)*dx;
	x += dx;
     
   }
   printf("Int(a=%f, b=%f)=%f", a, b, s);
   return 0;
}

[Stilet]

http://easylab.net.ua/chislennyie-me...etod-trapetsiy
http://cybern.ru/chislennoe-integrir...ciya-na-c.html
Ну и так далее поисковик рулит