function [S_n,normal_density]=S_n_star_compare_ND( n,p )
%S_N_STAR_COMPARE_ND compares S_n_star to the normal density for 


%a dummy variable to check against.
X=0:n;
X_star=(X-n*p)/sqrt(n*p*(1-p));

S_n=binopdf(X,n,p);

normal_density=(1/sqrt(2*pi))*exp(-(X_star).^2/(2));

%olots analysis.
figure(1); plot(X_star,S_n,X_star,normal_density);

end

