Ошибка в алгоритме - программа выводит n, вводимое с клавиатуры, а не из формул

#include <cstdlib>
#include <iostream>
#include <cmath>
float f (float x)
{
      return sin(x);
      }
      float Metod (float a, float b, int n)
      {
         float S, h, x1;
         int i;    
            S=f(a)+f(b);
            h=(b-a)/(2*n);
            x1=a+h;
            for (i=1;i<n;i++)
            {
            S=S+4*f(x1);
            x1=x1+2*h;
            }
            x1=a+2*h;
            for(int i=1;i<n;i++)
            {
            S=S+2*f(x1);
            x1=x1+2*h;
            }
            S=((b-a)/(6*n))*S;
            return S;
            }
            float Integral (float a, float b, int n)
            {
            float I1, d, I2, eps=0.001;
            I1=Metod(a,b,n);
            I2=Metod(a,b,2*n);
            d=fabs(I1-I2)/13;
            while (d>=eps)
            {
                  n=2*n;
                  I1=I2;
                  I2=Metod(a,b,2*n);
                  d=fabs(I1-I2)/13;
                  }
                  return I1;
                  }
using namespace std;
int main(int argc, char *argv[])
{
    system("chcp 1251");
    system("cls");
    float a, b, I; 
    int n;
    printf("Введите n: ");
    scanf("%d",&n);
    printf("Введите интервал a,b: ");
    scanf("%f",&a);
    scanf("%f",&b);
    I=Integral(a,b,n);
    printf("Интеграл равен =%f\n", I);
    printf ("Количество разбиений =%d\n",n);
    system("PAUSE");
    return EXIT_SUCCESS;
}

float Integral (float a, float b, int& n)