function binomial_distk = 0:40;pmf = binom_dist(40,0.5,k);plot(k,pmf,'b.');figure(2);cummDist = cumsum(pmf);plot(k,cummDist,'r.');end

function pmf = binom_dist(N,p,k) nValues = numel(k); pmf = zeros(1,nValues); for i = 1:nValues pmf(i) = nchoosek(N,k(i))*p^k(i)*(1-p)^(N-k(i)); endend

%define script and anonymous function in same file% data = 1:10; %# A vector of data% f = @(x) x.^2; %# An anonymous function% squaredData = f(data); %# Invoke the anonymous function

%%%%%%%%%OR%%%%%%%%%%%%

% define script and function in same file% function my_file% %script here% end% % function out = f(in)% %function here% end