圆面积精度问题

#include <cstdio>
#include <iostream>
using namespace std;

int main() {
    double a,b=3.14159;
    scanf("%lf", &a);
    printf("A=%.4lf", a*a*b);
    //cout << "A=" << a*a*b << endl;
    return 0;
}