モンテカルロ法による積分

###実現したいこと
モンテカルロ法を利用して積分の計算がしたいです。
範囲は1,0と仮定してよいです。
積分する関数は、f(x)=1-4(x-1/2)^2です。
###問題点
自分なりに作ってみたのですが、普通に手計算でした積分の結果とだいぶ誤差があるのでモンテカルロ法がなってない可能性があります。プログラムは間違っている部分があると思うのですがどこが違うでしょうか。
実行結果は以下なのですが、f(x)=1-4(x-1/2)^2を普通に計算したら**0.6666...**あたりになると思います。。。

乱数の数:10000
cnt=3389
面積の近似値は0.338900です.
-------------------------------
乱数の数:10000
cnt=3303
面積の近似値は0.330300です.

###コード

C
1//乱数による積分
2#include<stdio.h>
3#include<stdlib.h>
4#include<time.h>
5#include<math.h>
6
7//関数
8double f(double x)
9{
10    return (1-4*pow((x-1/2),2));
11}
12
13int main(void)
14{
15    int n; //乱数の数
16    int i,cnt=0;
17    double x,y,S;
18    printf("乱数の数:");
19    scanf("%d",&n);
20
21    //乱数をn個発生させる
22    srand((unsigned)time(NULL));
23    for(i=0;i<n;i++){
24        x=((double)rand()+1.0)/((double)RAND_MAX+1.0);
25        y=((double)rand()+1.0)/((double)RAND_MAX+1.0);
26        //printf("x=%f y=%f",x,y);
27        if(y<=f(x)){
28            cnt++;
29        }
30    }
31    printf("cnt=%d\n",cnt);
32
33    //積分=面積
34    S=(double)cnt/n;
35    printf("面積の近似値は%fです.\n",S);
36
37    return 0;
38}