已帮你解决了。。。。。。。

function [price sig] = mcmc(s0, K, sigma, r, T, N, M)
x=0;
for NN=N
rx=randn(M, NN-1);
dt=T/(NN-1);
S=[s0*ones(M,1) zeros(M,NN-1)];
for i=1:(NN-1)
S(:,i+1)=S(:,i)+S(:,i)*r*dt+sigma*dt^0.5*(S(:,i).*rx(:,i));
end;
x=x+1;
p1=max(mean(S(:,1:end),2)-K,0);
price(x)=exp(-r*T)*mean(p1);
sig(x)=std(p1)/M^0.5;
end
end

抱歉没及时回复，N是可以计算出一个区间的结果了，但是能不能让让所有变量可以像N这样呢？
请问应该怎么做呀？

我按照你的代码如法炮制了下，把NN换成了其它变量，比如做K的时候
function [price sig] = mcmc(s0, K, sigma, r, T, N, M)
x=0;
for KK=K
rx=randn(M, N-1);
dt=T/(N-1);
S=[s0*ones(M,1) zeros(M,N-1)];
for i=1:(N-1)
S(:,i+1)=S(:,i)+S(:,i)*r*dt+sigma*dt^0.5*(S(:,i).*rx(:,i));
end;
x=x+1;
p1=max(mean(S(:,1:end),2)-KK,0);
price(x)=exp(-r*T)*mean(p1);
sig(x)=std(p1)/M^0.5;
end
end
这样好像可以出结果，请问对吗？
如果对的话，我是不是能够直接用plot对函数输出结果作灵敏度分析的折线图呀

做MATLAB公式时，都需要进行标注与解释，对任何一个变量都需要详细说明才行。一个数与一个向量是不一样的。如果按您标注成KK=K，那么K最好定义成一个向量，且这个向量中每个值压根是正的整数比如K=1：30，或K=80：100.当然，你直接写成K=1也不会错。

是可以做图的。。。。。。