Вычислить значения интегралов по двум формулам (задача решена, проблема с рекурсией)

#include <stdio.h>
#include <stdlib.h>
#include <cmath>
#include <locale.h>
 
int func(float a, float b, float y, int* N);
 
int main(void)
{
    setlocale(LC_ALL, "");
 
    float a = 0.5, b = 2.5, y = 0;
    int N = 10; 
    int* pN = &N;
    func(a, b, y, pN);
}
 
int func(float a, float b, float y, int* N)
{
    float mas[1001] = { 0 }, I1 = 0, I2 = 0, Ib = 0;
    int i = 0;
 
//x меняется с шагом 0.2  от 0.5 до 2.5
//вычисляется y и записывается в массив
    for (float x = a; x <= b + (b - a) / *N; x += (b - a) / *N)
    {
        y = pow(x, 2) * sqrt(pow((2 * x + 3), 3));
        mas[i] = y;
        /*printf("\nmas[%d] = %f", i, mas[i]);*/
        i++;
    }
 
//нужно сравнить результаты ф. трапеции и ф.Симпсона
 
/////////////////////////////////////////вычисление формул/////////////////////////////////////////////
 
    I1 += (mas[0] / 2) + (mas[*N] / 2); //к формуле трапеций
 
    for (int j = 1; j < *N; j++)
    {
        I1 += mas[j];  // к формуле трапеций
        if ((j > 0) && (j < *N) && (j % 2 != 0)) I2 += mas[j];       // к формуле Симпсона
        if ((j > 1) && (j < *N - 1) && (j % 2 == 0)) Ib += mas[j]; // к формуле Симпсона
    }
 
 
    I1 *= (b - a) / *N; // Итог ф. трапеций
    I2 = ((b - a) / *N) / 3 * (mas[0] + mas[*N] + 4 * I2 + 2 * Ib); // Итог ф. Симпсона
 
/////////////////////////////////////////////////////////////////////////////////////////////////////////////
 
 
 
//вывод результата:
    if (I1 > I2)
        printf("При N = %p\nФормула трапеций: %f  >  Формула симпсона: %f\n\n", N, I1, I2);
    else
        printf("При N = %p\nФормула трапеций: %f  <  Формула симпсона: %f\n\n", N, I1, I2);
 
 
//рекурсия
 
    if (*N == 1000) { getchar();  return 0; }
 
    else { *N *= 10; return func(a, b, y, N); }
}