hw.m

function [ res ] = ie505_hw1( n,a,b )
%İE505_HW1 Summary of this function goes here
%   Detailed explanation goes here

[~,bins]= hist([a,b],1000);

r =unique( round( n*bins));

res = arrayfun(@(x)nchoosek(n,x)*double(1/(2^n)),r);

res = sum(res)
% [~,k] = hist([0,1],length(res));
% plot(res)
% [sigma, mu] = gaussfit(k,r);