一題第一問首先在x,y的取值范圍內取若干個點，用linspace命令，如對于x則是x1=linspace(0,1,200)括號內前兩個數是x的取值范圍，最后一個數是分成的份數用meshgride命令將（x1,y1）從矢量轉化為x,y矩陣，具體的函數是x,y=meshgride(x1,y1)設一個函數z簡寫將要積分的函數，即z=cos（x）+sin(y)以便以后畫圖時方便，沒有這一步也是可以的最后一步就是利用mesh函數畫出三維圖形，具體命令為mesh(x1,y1,z)前面兩個是我們自己設的變量，順序可以換，最后一個是被積函數，如果前面沒有將它簡寫為z的話，可以直接寫原函數總結：這一問具體的操作為x1=linspace(0,1,200);y1=linspace(1,2,200); x,y=meshgrid(x1,y1); z=cos(x)+sin(y); mesh(x1,y1,z)或者x1=linspace(0,1,200);y1=linspace(1,2,200);x,y=meshgrid(x1,y1); mesh(x1,y1,cos(x)+sin(y)作出圖形如下meshz 的作用和mesh完全一樣，只是在下方加上了一些線而構成了矩形的方框，meshz截圖如下waterfall 和mesh功能完全一樣，只是會產生動態瀑布的效果 waterfall截圖 surf用法也和mesh類似，只是它做出的圖是黑色的，所以取點時不要太密，一般要少用100個，surf所做三維圖形截圖Ps:以上函數最后一行均可簡寫為函數名(z)不影響結果，如mesh(z)等價于mesh(x1,y1,z)plot用來畫三維曲線圖而不是曲面圖，它的功能就是將所有的x,y,z坐標表示的點連起來，plot3 做的三維曲面截圖 要想畫曲面圖應該用以上幾個一題第二問dblquad 用于進行數值積分，他的用法是：設一個函數名q ，令 q=dblquad（x,y）被積函數，x下限，x上限，y下限，y上限）即可算出 q=dblquad(x,y)(cos(x)+sin(y),0,1,1,2)q =1.7979fun=(x,y)cos(x)+sin(y)fun = (x,y)cos(x)+sin(y)q=dblquad(fun,0,1,1,2)q =1.7979f=cos(x)+sin(y) f = cos(x) + sin(y) g=int(f,x,0,1) g = sin(1) + sin(y) h=int(g,y,1,2) h = cos(1) - cos(2) + sin(1)fun=inline(cos(x)+sin(y),x,y)fun = Inline function: fun(x,y) = cos(x)+sin(y) jf=dblquad(fun,0,1,1,2) jf=dblquad(inline(cos(x)+sin(y),x,y),0,1,1,2)其中老師講的和不知道是不是版本不兼容問題，總是做不出來一題第三問 n=1000;x=unifrnd(0,1,n,1);y=unifrnd(1,2,n,1);z=unifrnd(0,2,n,n);frq=0;for i=1:nfor j=1:nzt=cos(x(i)+sin(y(j);if z(i,j) p=frq/n2p = 0.8969 jf=2*pjf = 1.7938第二題n=10000;frq=0;for i=1:nx1=randperm(6);x2=randperm(6);x3=randperm(6);x=x1(1)+x2(1)+x3(1);if x10frq=frq+1;endend p=frq/np =0.5036改變n的數字n=100000;frq=0;for i=1:nx1=randperm(6);x2=randperm(6);x3=randperm(6);x=x1(1)+x2(1)+x3(1);if x10frq=frq+1;endend p=frq/np = 0.5014第三題答案x=unifrnd(0,2*pi,1000,1);y=x.2.*sin(x).*cos(x);max(y)ans = 7.9543 min(y)ans = -15.3593對三題中所描述的曲線進行繪圖如下x2=linspace(0,2*pi,1000);y2=x2.2.*sin(x2).*cos(x2);plot(x2,y2)set(gca,ytick,-15,0,8)grid on圖形為第四題在已知股票價格和行權價格，無風險利率，標準差，到期天數后，到期日現金流可以簡化為Max0,50exp(2+4)-52n=10000