function [A,a,E]= myfun(T,K)logT=log(T);logK=log(K);daoT=T.^(-1);lenT=length(T);A=ones(3);x=[];for r=1:lenT-2    A(:,2)=logT(r:r+2);    A(:,3)=daoT(r:r+2);    b=logK(r:r+2);    xt=A\b;    xt=xt';    x=[x;xt];endlogA=mean(x(:,1));A=exp(logA);a=mean(x(:,2));E=mean(x(:,3));h1=stem(T,K,'bo');    set(h1,'MarkerFaceColor','green');    set(h1,'LineStyle','none');    set(get(h1,'BaseLine'),'LineStyle','none');    hold on;for t=200:10:1300    k=A*(t^a)*exp(E/t);plot(t,subfun(t,A,a,E));    function k = subfun(t,A,a,E)k=A*(t^a)*exp(E/t);