clearclc


load?stockpaths;
S=stockpaths;??%各列累乘

N=3;
M=8;
K=1.1;
ex_time=N*ones（M1）;
valueF=zeros（size（S））;
valueF（:end）=max（K-S（:end）0）;
%時間間隔為一時的折現因子為0.94176，所以單個rt的值為rt=log（0.94176）=-0.06;
valueF（:1）=exp（-0.06*3）*valueF（:end）;

for?i=N:-1:2
????index=?find（S（:i）????X=S（indexi）;
????Y=valueF（indexi+1）*0.94176;
????A=[ones（size（X））?X?X.^2];
????beta=A\Y;
????conExp=A*beta;
????if?i