{Вычисление определенного интеграла
 с помощью формулы Симпсона
 I=(b-a)/(3*n) *(y[0]/2+y[2n]/2+y[2]+y[4]+..+y[2*(n-1)+]
 +2*(y[1]+y[3]+..+y[2*n-1]))}

 [IA x=x+d func I=I+2*y x=x+d func I=I+y k=k+1][IB]  
 [simp_int d=(b-a)/(2*n) x=a func I=y/2
 x=b func I=I-y/2 x=a k=0 loop(IA,n-k,IB)
 I=I*(b-a)/(3*n)]

 {Теперь подставляем функцию, задаем
 пределы интегрирования (a,b) и число
 отрезков n. Затем вычисляем:}
 [func y=1/(1+x^2)]
 a=0 b=1 n=12 simp_int I=?10
 {Проверка:}
          z=arccos(-1)/4=

 Результат:

  Вычисление определенного интеграла
 с помощью формулы Симпсона
 I=(b-a)/(3*n) *(y[0]/2+y[2n]/2+y[2]+y[4]+..+y[2*(n-1)+]
 +2*(y[1]+y[3]+..+y[2*n-1])) 

 [IA x=x+d func I=I+2*y x=x+d func I=I+y k=k+1][IB]  
 [simp_int d=(b-a)/(2*n) x=a func I=y/2
 x=b func I=I-y/2 x=a k=0 loop(IA,n-k,IB)
 I=I*(b-a)/(3*n)]

  Теперь подставляем функцию, задаем
 пределы интегрирования (a,b) и число
 отрезков n. Затем вычисляем: 
 [func y=1/(1+x^2)]
 a=0 b=1 n=12 simp_int I=0.7853981633
  Проверка: 
          z=arccos(-1)/4=0.7853981634