Ошибка при выводе результата

#include <cstdlib>
#include <iostream> 
#include <string>
 
using namespace std;
 
int i, j, n;
double h;
double func(double x, double y);
double integr(double x0, double y0, double xn, int n);
int main()
{
    double x0, y0, xn;
    setlocale(0, "");
    cout << "х0=";
    cin >> x0;
    cout << endl;
    cout << "y0=";
    cin >> y0;
    cout << endl;
    cout << "xn=";
    cin >> xn;
    cout << endl;
    cout << "n=";
    cin >> n;
    cout << endl;
 
    integr(x0, y0, xn, n);
    cout << "result: " << integr << endl;
 
    system("pause");
 
}
double func(double x, double y)
{
    if ((pow(y, 3) + x + 1) != 0)
        return ((3 * pow(y, 2)) / (pow(y, 3) + x + 1));
}
 
double integr(double x0, double y0, double xn, int n)
{
    double h = (xn - x0) / n;
 
    double* x = new double[n + 1];
    double* y = new double[n + 1];
    double* y1 = new double[n + 1];
 
    for (i = 0; i <= n; i++)
    {
        x[i] = x0 + i * h;
        cout <<"h"<< i << " = " << x[i] << endl;
    }
 
    y[0] = y0;
    for (i = 1; i <= n; i++)
    {
        y1[i] = y[i - 1] + h * func(x[i - 1], y[i - 1]);
        y[i] = y[i - 1] + h * (func(x[i - 1], y[i - 1]) + func(x[i], y1[i])) / 2.0;
    }
    delete[] x;
    delete[] y;
 
    return 0;
}