code:
 
%в эту переменную импортируется выборка
x=sigma1; 
>> % Плотность распределения Парето 
p=inline('(x^gamma)/(1+x^gamma))',... 'x','gamma'); 
% Функция правдоподобия 
L=inline('-prod(feval(p,x,theta))',... 
'theta','p','x'); 
% Наиболее правдоподобное значение параметра 
m=fminsearch(L,[1],[],p,x) 

code:

Error using ==> inlineeval at 15 
Error in inline expression ==> -prod(feval(p,x,theta)) 
 Error in inline expression ==> exp(-mu)*mu.^x/factorial(x(1)) 
 N must be a matrix of non-negative integers. 

Error in ==> inline.subsref at 27 
    INLINE_OUT_ = inlineeval(INLINE_INPUTS_, INLINE_OBJ_.inputExpr, 
    INLINE_OBJ_.expr); 

Error in ==> fminsearch at 205 
fv(:,1) = funfcn(x,varargin{:}); 
  

code:
 
% Плотность распределения Пуассона
p=inline('exp(-mu)*mu.^x/factorial(x(1))',...
'x','mu');
% Функция правдоподобия
L=inline('-prod(feval(p,x,theta))',...
'theta','p','x');
% Наиболее правдоподобное значение параметра
m=fminsearch(L,[1],[],p,x)
m = 5.6500
% График функции правдоподобия
mu2=m-2:0.1:m+2;
 
for i=1:length(mu2), L2(i)=L(mu2(i),p,x); end
L2=L2/L(m,p,x); plot(mu2,L2)