Неверный вывод одной из подпрограмм, реализующей расчет формулы

#include "stdafx.h"
#include <Windows.h>
#include <math.h>
 
int main()
{
    SetConsoleCP(1251);//Первый способ
    SetConsoleOutputCP(1251);
    double formula(int a, int b,double &s);
    double R = 0, r = 0, S = 0, s=0;
    printf("Введите R,r:", R, r); //
    scanf("%f%f", &R, &r);
    S = formula(R,r,s);
    printf("Ответ: %f \n", S);
    system("pause");
    return 0;
}
 
double formula(int a, int b, double &s)
{
    double rez = 0;//Если закомментить расчет и поставить в rez = 1; все выводит. если оставить расчет формулы, то выдаёт 0
    rez = M_PI*(pow(a, 2) - pow(b, 2));
    /*printf("%3.1f^2-%5.1f^2=%3.1f\n", a, b, rez);*/
    system("pause");
    return rez;
}
 
//double rs(double R, double r)//второй способ
//{ 
//  return 3.14159265*(R*R - r*r); 
//}
//int main() 
//{
//  float S, R, r; 
//  printf("R r ? "); 
//  scanf("%f%f", &R, &r); 
//  S = rs(R, r);
//  printf("S=%f", S); 
//  system("pause");
//  return 0;
//}

double sposob1(double a, double b)
{
 return M_PI*(a*a - b*b);
}
 
int main()
{
    SetConsoleCP(1251);//?????? ??????
    SetConsoleOutputCP(1251);
    double R = 0, r = 0, S = 0;
    printf("Введите R,r: ");
    scanf("%lf %lf", &R, &r);
    S = sposob1(R,r);
    printf("%3.2lf^2-%3.2lf^2=%3.2lf\n", R, r, S);
    system("pause");
    return 0;
}