clc

clear all

x=gevrnd(0.1,3,20,[100,1]);

%x=xlsread('ozone.xls');

CDF= @(x,k,sigma,mu,lamda)((1./exp((1 - (k*(mu - x))./sigma).^(-1/k))).^lamda-1)./lamda + 1

PDF=@(x,k,sigma,mu,lamda)(1./exp(1./(1 +(k*(x-mu)./sigma)).^(1./k))).^(lamda - 1).*1./(sigma*exp(1./(1 +(k*(x-mu)./sigma)).^(1./k)).*(1 + (k*(x-mu)./sigma)).^(1./k + 1))

START=[0.1,3,20,1];

est2= mle(x,'pdf',PDF,'cdf',CDF,'start',START)

but when i give lamda =real number

the code run and estimate the another paramter