function [f]=acceptanceloop(trials);

%Acceptance.m: The input argument to this function is a function 
%This function uses acceptance-rejection to generate the distribution
%P(x)=.75*(1-x^2) where -1<=x<=1

%Largest value the probability will reach

for i=1: trials
    
while j=0
    xtry=2*rand-1;
    if (1-xtry^2)>=rand, f(i)=xtry;j=1;
    else j=0;
    
end