K=95;
T=0.25;
S0=100;
r=0.1;
N=400;
sigma=0.5;
dt= T/N;
u=exp(sigma * sqrt(dt));
d=1/u;
q=(exp(r*dt)-d)/(u-d);
lat = zeros(N+1,N+1);
for j=0:N
lat(N+1,j+1)=max(0 ,S0*(u^j)*(d^(N-j))-K);
end
for i=N-1:-1:0
for j=0:i
lat(i+1,j+1) = exp(-r*dt) *(q * lat(i+2,j+2) + (1-q) * lat(i+2,j+1));
end
end
lat;
amcall= lat(1,1)