Creo que donde esta el error (por asi decirlo) es en los:
x[i+1]=x[i]+dx;
fx[i+1]=(x[i+1]*x[i+1])+1;
nArea[i+1]=dx*(fx[i]+fx[i+1])/2;
areaTotal = areaTotal + nArea[i+1];
mas especifico los:
ya que por ejemplo le das un tamaño a los vectores de
n+1 los cuales tienen las posiciones desde
0 hasta
n, entonces mientras
i<n todo va bien, pero cuando llega a
i=n al usar
i+1 estas intentando acceder a la posicion
n+1 del vector, que no existe :S
Saludos!!