?? expl.m
字號:
function values=expl(Smax,M,T,N,r,sig,X)
dS=Smax/M;
dt=T/N;
matf=zeros(M+1,N+1);
vetS=linspace(0,Smax,M+1)';
matf(:,N+1)=max(X-vetS,0);
matf(1,:)=X;
matf(M+1,:)=0;
veti=0:N;
vetj=0:M;
a=0.5*(-r*dt*vetj+sig^2*dt*(vetj.^2))/(1+r*dt);
b=(1-sig^2*dt*(vetj.^2))/(1+r*dt);
c=0.5*(r*dt*vetj+sig^2*dt*(vetj.^2))/(1+r*dt);
for i=N:-1:1
for j=2:M
matf(j,i) = a(j)*matf(j-1,i+1) + b(j)*matf(j,i+1)+ ...
c(j)*matf(j+1,i+1);
matf(j,i) = max(matf(j,i),X-(j-1)*dS);
end
end
values = flipud([vetS, matf(:,1)]);
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -