3# ywh19860616 sigma=0.039974; a=0.49994; u=0.023485; lemda=0.020000; gamma=((a+lemda)^2+2*sigma^2)^(1/2); r1=0.0234; t=[1,2,3,4,5]; P=(2*gamma*exp((a+lemda+gamma)*t/2)/((a+lemda+gamma)*(exp(gamma*t)-1)... +2*gamma))^(2*a*u/(sigma^2))*exp(-(2*(exp(gamma*t)-1)/((a+lemda+gamma)... *(exp(gamma*t)-1)+2*gamma))*r1); P1=(2*gamma*exp((a+lemda+gamma)*(t+1)/2)/((a+lemda+gamma)*... (exp(gamma*(t+1))-1)+2*gamma))^(2*a*u/(sigma^2))*... exp(-(2*(exp(gamma*(t+1))-1)/((a+lemda+gamma)*(exp(gamma*(t+1))-1)+2*gamma))*r1); R=-1/t*log(P); fr=log(P)-log(P1); ??? Error using ==> mrdivide Matrix dimensions must agree.