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

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

 {Теперь подставляем функцию, задаем
 пределы интегрирования (a,b) и число
 отрезков n. Затем вычисляем:}
 [func y=x^3+sin(x)]
 a=0 b=1 n=1000 integral I=?

 [func y=x]
 a=0 b=1 n=1 integral I=?

 Результат:

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

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

  Теперь подставляем функцию, задаем
 пределы интегрирования (a,b) и число
 отрезков n. Затем вычисляем: 
 [func y=x^3+sin(x)]
 a=0 b=1 n=1000 integral I=0.709698

 [func y=x]
 a=0 b=1 n=1 integral I=0.500000