C++, вычисление значения интеграла по формуле трапеций

[Кирилл335632₽]

Помогите пожалуйста написать программу для шестого задания, заранее благодарю

[digitalis]

Помогаю.
Делишь интервал на много частей - допустим, 64. Вычисляешь и суммирушь значения функции в начале и конце интервала, делишь затем на 2. В цикле вычисляешь значения в промежуточных точках a+h, a+2h и т. д. Добавляешь их в сумматор .
y = (f(a) + f(b) )/2 + f(a+h)+ .....
Умножаешь на h - и наслаждаешься результатом. Потом, слегка подумав, удваиваещь n и делаешь то же самое. Если разница предыдущего расчёта и следующего по абс. величине больше епс - удваиваешь ещё раз и так - до результата.
Получаешь заслуженную нами пятёрку и говоришь преподу: "Это ещё что! Я вот то же самое могу по и Симпсону отчекрыжить!"

[Кирилл335632₽]

Спасибо большое

[digitalis]

Спасибо скажешь, когда воплотишь это в коде и представишь нам на поругание восхищение.

[Кирилл335632₽]

{
	float a, b, h, x, k, yx, S, I, Ip;
	int i, n;
	float const eps=1e-3;
	printf("\n a=");
	scanf("%f", &a);
	printf("\n b=");
	scanf("%f", &b);
	n=10;
    I=0;
    Ip=0;
	while(abs(k)>0)
	{
		h=(b-a)/n;
		for(x=a, S=0; x<b; x+=h)
	    {
		  yx=1/(1+exp(x/2));
		  S+=yx;
        }
        I = h*((1/(1+exp(b/2)) - 1/(1+exp(a/2))) / 2 + S); 
        k = (I - Ip);
        Ip = I;
	} 
	printf ("\n I=%7.3f",I);

Почему то выдается что I=0.000
Что может быть?

[BDA]

Используйте кнопку CODE для оформления кода. eps нигде не задействовали. Перед циклом while переменная k не имеет значения (инициализируйте чем-нибудь). Скорее всего цикл while ни разу не выполнился. Ну и есть мелкие недоработки: нужно суммировать f(a) и f(b), цикл for надо начинать с a + h, а не с a.

[digitalis]

	float a, b, h, x, k, yx, S, y, yp;
	int i, n;
	float const eps=1e-3;
	printf("\n a=");
	scanf("%f", &a);
	printf("\n b=");
	scanf("%f", &b);
	n=10;
        k=1000;
   	while (abs(k)>eps)  
	{
	    h=(b-a)/n;
            S = (1/(1+exp(a/2)) + 1/(1+exp(b/2)) ) / 2 ;
            x = a+h ;
	    for(i=1; i<n; i++)
	       {
		  yx=1/(1+exp(x/2));
		  S+=yx;
                  x += h ;
               }
             y = h*S;
             if (n>10) k = y - yp;
             yp = y ;  
	     printf ("\n n=%6d y=%7.3f",n,y);
             n = n+n;
  	}

Я бы так написал... Во время первого прохода k не вычисляем - не с чем сравнивать. Большим начальным значением k гарантируем ещё один проход, как минимум.
Почему выдало I=0.000 ? Ну кто же старательно вычисленное float запихивает в int ? Там было небольшое значение, и в int оно обрезалось до 0.
---------------------------
Ну и ?

[T4UV0R]

Спасибо большое

[digitalis]

Задание от Кирилла, а мерси - от T4UV0R ? Один в двух лицах, как Ляпис-Трубецкой ?