double f(double x)
{
    return x - x*x + x*x*x - x*x*x*x + sin(x*13.)/13.;
}

void example_04()
{
    TRandom3 rnd;
    
    TH1F *h1 = new TH1F("h1","A Random Histogram",120,-0.1,1.1);
    
    for(int i=0;i<100000;i++) {
        double x,y;
        do {
            x = rnd.Rndm();
            y = rnd.Rndm()*0.45;
        }while (y>f(x));
        
        h1->Fill(x);
    }
    
    h1->SetFillColor(50);
    h1->Draw();
}