x=data(1:1758,2:4);
pdf=@(x,gamma1,gamma2,gamma3,segma1,segma2,segma3,alpha)((1./((gamma3.*(x(:,3) - 6135993778780923./1048576))./segma3 + 1).^(alpha./gamma3) + 1./((gamma2.*(x(:,2) - 17581304132125./8192))./segma2 + 1).^(alpha./gamma2) + 1./((gamma1.*(x(:,1) - 5055429226824063./524288))./segma1 + 1).^(alpha./gamma1)).^(1./alpha - 3).*(3.*alpha.*(1./((gamma3.*(x(:,3) - 6135993778780923./1048576))./segma3 + 1).^(alpha./gamma3) + 1./((gamma2.*(x(:,2) - 17581304132125./8192))./segma2 + 1).^(alpha./gamma2) + 1./((gamma1.*(x(:,1) - 5055429226824063./524288))./segma1 + 1).^(alpha./gamma1)).^(1./alpha) - 3.*alpha - 3.*(1./((gamma3.*(x(:,3) - 6135993778780923./1048576))./segma3 + 1).^(alpha./gamma3) + 1./((gamma2.*(x(:,2) - 17581304132125./8192))./segma2 + 1).^(alpha./gamma2) + 1./((gamma1.*(x(:,1) - 5055429226824063./524288))./segma1 + 1).^(alpha./gamma1)).^(1./alpha) + 2.*alpha.^2 + (1./((gamma3.*(x(:,3) - 6135993778780923./1048576))./segma3 + 1).^(alpha./gamma3) + 1./((gamma2.*(x(:,2) - 17581304132125./8192))./segma2 + 1).^(alpha./gamma2) + 1./((gamma1.*(x(:,1) - 5055429226824063./524288))./segma1 + 1).^(alpha./gamma1)).^(2./alpha) + 1))./(segma1.*segma2.*segma3.*exp((1./((gamma3.*(x(:,3) - 6135993778780923./1048576))./segma3 + 1).^(alpha./gamma3) + 1./((gamma2.*(x(:,2) - 17581304132125./8192))./segma2 + 1).^(alpha./gamma2) + 1./((gamma1.*(x(:,1) - 5055429226824063./524288))./segma1 + 1).^(alpha./gamma1)).^(1./alpha)).*((gamma3.*(x(:,3) - 6135993778780923./1048576))./segma3 + 1).^(alpha./gamma3 + 1).*((gamma2.*(x(:,2) - 17581304132125./8192))./segma2 + 1).^(alpha./gamma2 + 1).*((gamma1.*(x(:,1) - 5055429226824063./524288))./segma1 + 1).^(alpha./gamma1 + 1))
[phat,pci]=mle(x,'pdf',pdf,'start',[0.5 0.5 0.5 0.5 0.5 0.5 0.5],'lowerbound',[0 0 0 0 0 0 0],'upperbound',[inf inf inf inf inf inf inf])