Формула Симпсона

[Raz0r]

Подскажите пожалуйста, как по формуле Симпсона вычислить площадь фигуры, заключенной между функциями? Саму формулу я нашёл, но как её применить к данным функциям, ума не приложу.
Функции:
y=-x^2+2x+2
y=-x^2-4x-1
y=3
Формула Симпсона

[JTG]

Кто-нить помнит, как интегралы в паскале считать? Хых, сессия скоро, "а моя совсем отупела однако"

[_Dmitry]

Нужно вычислить площадь заштрихованной фигуры. Я бы сделал так: спустил параболу вниз на 3, перевернул ее, вычислил площадь желтого треугольника и умножил на 2, т.к. красная линия - ось симметрии.
Вот код с контролем точности, написан давно, сейчас заменил только функцию.
В данной задаче, контроль точности не нужен, т.к. для квадратной параболы формула Симпсона дает точное значение. Поэтому цикл repeat until можно убрать.

var
  i,n: integer;
  d,S0,S,x_min,x_max,eps: real;

function F(x: real): real;
begin
  F:=x*x+4*x+4;
end;

begin
  eps:=0.0001;
  x_min:=-2; x_max:=-0.5;
  n:=1; S:=0;
  repeat
    n:=n*2; d:=(x_max-x_min)/n; S0:=S;
    S:=F(x_min)+F(x_max);
    for i:=1 to n-1 do
      if odd(i) then S:=S+4*F(x_min+i*d)
                else S:=S+2*F(x_min+i*d);
    S:=S*(x_max-x_min)/3/n;
  until abs(S-S0) <= eps;
  S:=2*S;
  writeln(S:1:4);
  readln;
end.

[Raz0r]

Спасибо, только я совсем забыл написать, что надо под BС...

[_Dmitry]

На C будет так

#include <stdlib.h>
#include <stdio.h>
#include <conio.h>

float F(float x)
{return x*x+4*x+4;}

main()
{ int n;
  float d,S0,S,x_min,x_max,eps;
  eps=0.0001;
  x_min=-2; x_max=-0.5;
  n=1; S=0;
  do {
    n=n*2; d=(x_max-x_min)/n; S0=S;
    S=F(x_min)+F(x_max);
    for (int i=1; i<n; i++)
      if (i%2 != 0) {S=S+4*F(x_min+i*d);}
               else {S=S+2*F(x_min+i*d);}
    S=S*(x_max-x_min)/3/n;}
  while (abs(S-S0) <= eps);
  S=2*S;
  printf("%f",S);
  getch();
}

[Raz0r]

Спасибо, будем разбираться.