Vielen dank, Jack

double Trapez(double l, int w)
{
	double x = 0,I,h;
	int s=0;
	h = (b-a) / l;
	for (int i=1; i <= n-1; i++)
	{
		x = a + i * h;
		s = s + f(x);
	}
	I = (h/2) * ( f(a) + 2*s + f(b));
	return I;
}

void Romberg()

{
	double T[4] [4];
	for (int i = 0; i <= m; i++)
	{
		T[i][0] = Trapez( (i+n), int(pow(2.,n+i)) );
	}

	cout << endl;
	for (int k = 1; k <= m; k++)
	{
		for (int i = k; i <= m; i++)
		{
			T[i][k] = (( T[i][k-1]-T[i-1][k-1] ) /(pow(2.,2.*k)-1));
		}
	}
	//Ausgabe der Iterationen
	for (int i = 0; i <= m; i++) 
	{
		cout <<"T["<<i<<"]["<<i<<"]="<<T[i][i];
			cout <<"\n";
	}
	return;
}

PS: Ich habe nie behauptet, dass ich nur halbwegs programmieren kann

Zweitens: Das Array T könnte möglicherweise etwas zu klein gewählt sein.

int m=3;
double a=-10, b=10, n=3; //Intervall -10 bis 10, Parameter n=m=3

double f(double x)
{
	double pi=3.141592654;
    return (1/sqrt(2*pi)*exp((-x*x)/2));
}

double Trapez(double h, int w)
{
	double temp = 0.;
	for (int v = 0.; v < w; v++)
	{
		temp += h * ( f(a+v*h) + f(a+(v+1)*h)) / 2;
	}
	return temp;
}

void Romberg()

{
	double T[10] [10];
	for (int i = 0; i <= m; i++)
	{
		T[i][0] = Trapez( (b-a)/(pow(2.,i)), int(pow(2.,i))  );
	}

	cout << endl;
	for (int k = 1; k <= m; k++)
	{
		for (int i = k; i <= m; i++)
		{
			T[i][k] = (( T[i][k-1]-T[i-1][k-1] ) /(pow(2.,2.*k)-1));
		}
	}
	//Ausgabe der Iterationen
	for (int i = 0; i <= m; i++) 
	{
		cout <<"T["<<i<<"]["<<i<<"]="<<T[i][i];
			cout <<"\n";
	}
	return;
}

Ich weiss nur, dass das korrekte Ergebnis 1.0096446526939107 lauten müsste.
Wäre schön, wenn du den Fehler finden würdest.