美式期权执行边界的matlab实现

计算美式期权价格，并画出它的执行边界。%用for循环求出各个结点处的欧式看涨期权的价值 %以上用倒推的方法并考虑折现率来求出欧式看涨期权的精确值，所得出的矩阵EFX即为所求1.%比较每个节点处提前执行和不提前执行的价值，确定美式期权的内在价值,包括最后一列%通过将未来收益折现和当期执行的比较后，得出美式看跌期权的精确值，所得矩阵AX即为所求2.%增多结点数来画出执行边界

%衍生品作业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循环求出各个结点处的欧式看涨期权的价值
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） %以上用倒推的方法并考虑折现率来求出欧式看涨期权的精确值，所得出的矩阵EFX即为所求1.
 for b=1:n+1
    for a=1:b
      afx（ab）=max（k-mx（ab）0）;%美式期权的内在价值包括最后一列
    end;
end;
for a=1:n+1
    ax（a