function [sample] = epanech(parameter,method) cdm1 = 1/(pi^2); %cdm1 = 1/((4/3)*pi); sz = size(parameter); x = sqrt(parameter*parameter'); if(method==0) sample = 1/2*cdm1*(sz(2)+2)*(1-x); end if(method==1) %Only dealing with two dimensions here so Epanechnikov reduces to: if(x < 1) sample = 1/2*cdm1*(sz(2)+2)*(1-x); end if(x >= 1) sample = 0; end end