нада определить интеграл

[olegraperos]

Задали функцию. Написать программу для нахождения определенного интеграла с точностью Е для заданных границ [a, b] с помощью трех методов: 1) прямоугольников, 2) трапеций, 3) Симпсона.
a=0,b=ln2 ; функция виглядит так sinh(x)^4(в 4 степени)......
Вот мой код:

#include<iostream>
#include<math.h>
#include<conio.h>
using namespace std;
double func(double);
double pryam(double,double,double,int);
double trap(double,double,double,int);
double simp(double,double,double,int);
  void main(){
	double y,n,a=0, b=log(2.0),d,f,m;
cout<<"vvedit tochnist=";
	cin>>d;
cout<<"Vvedit promishok=";
	cin>>n;
		y=pryam(a,b,d,n);
		m=trap(a,b,d,n);
		f=simp(a,b,d,n);
	cout<<"Metodom pryamokytnuka="<<y<<"\nMetodom trapecii="<<m<<"\nMetodom simpsona="<<f;
	_getch();}
double pryam(double a,double b,double d,int n){
double x=a, h, S=0 ,e,I2,I;
h=(a-b)/n;
for(int i=0; i<=n-1; i++){
    S+=func(x+h/2); x+=h;}
  return h*S;}
double trap(double a,double b,double d,int n){
 double x=a, h=(b-a)/n, S=0;
  for(int i=0; i<=n-1; i++){
    S+=func(x)+func(x+h); x+=h; }
  return h/2 * S; }
double simp(double a,double b,double d,int n){
  double x=a, h=(b-a)/n, S=0;
  int m=n/2;
  for(int i=0; i<=m-1; i++){
    S+=func(x)+4*func(x+h)+func(x+2*h); x+=2*h;}
  return h/3 * S;}
double func(double x){
	return pow(sinh(x),4.0);}

Сделал нахождение площи но как задать с точностю не пойму..