円周率を求める

pi.sas

data _null_;
    do i = 1 to 100000000;
        x = rand("uniform");
        y = rand("uniform");
        if x ** 2 + y ** 2 <= 1 then cnt + 1;
    end;
    pi = cnt * 4 / 100000000;
    put pi = ;
run;