%衍生品作業(yè)1
s=100;
k=95;
r=0.1;
sigma=0.5;
t=0.5;
n=2000;
q=0;
dt=t/n?;
u=exp（sigma*sqrt（dt））?;
d=1/u;??%選擇ud滿足的關系
pu=（exp（r*dt）-d）/（u-d）?;
pd=1-pu?;
for?b=1?:?n+1
??for?a=1?:b?
?????mx（ab）=s*power（ub-a）*power（da-1）;
?????efx（ab）=max（mx（ab）-k0）;?%給出歐式期權價格
????end;
end?;?%用for循環(huán)求出各個結點處的歐式看漲期權的價值
for?c=1:n?
???b=n+1-c;
???for?a=?1:b
???efx（ab）=（efx（ab+1）*pu+efx（a+1b+1）*pd）*exp（-r*dt）;
??end;
end;
?efx（11）?%以上用倒推的方法并考慮折現(xiàn)率來求出歐式看漲期權的精確值，所得出的矩陣EFX即為所求1.
?for?b=1:n+1
????for?a=1:b
??????afx（ab）=max（k-mx（ab）0）;%美式期權的內(nèi)在價值包括最后一列
????end;
end;
for?a=1:n+1
????ax（a