data=sort(ft);
rc=10
tc=data(rc)
tt=10e3
rt=10
SE=sum(data)
SC=sum(data(1:rc))+(length(data)-rc)*tc
ST=sum(data(1:rt))+(length(data)-rt)*tt

x1=min(data)
xn=max(data)
xc=(x1+xn)/2

W=xn-x1
m=round(sqrt(length(data)))
w=W/(m-1)

u(1)=x1-w/2;
u(m+1)=xn+w/2;

for i=2:length(data)
    u(i)=u(i-1)+w;
    if u(i)>u(m+1)
        break
    end
end
u(end)=[];

for i=1:length(u)-1
    xcv(i)=(u(i+1)+u(i))/2;
end

dr=zeros(1,length(u)-1);
sumdr=zeros(1,length(u)-1);

for i=1:length(u)-1
    for j=1:length(data)
        if data(j)>=u(i) && data(j)<u(i+1)
            dr(i)=dr(i)+1;
        end
    end
end

sumdr(1)=dr(1);

for i=2:length(u)-1
    sumdr(i)=dr(i)+sumdr(i-1);
end

for i=2:length(u)-1
    ni(i)=dr(i-1);
end

R=ni./length(data);

dt=w;

for i=1:length(u)-1
    lambdaCI(i)=dr(i)/ni(i)/w;
end

MTBFE=SE/length(data)
MTBFC=SC/rc
MTBFT=ST/(rt+1)

s=sqrt(sum(data-MTBFE).^2/length(data)-1)

lambdaC=1/MTBFC
lambdaT=1/MTBFT