這段**具有很好的互動性和通用性，將**複製到matlab編輯器之後，按照提示操作即可。傻瓜式操作，結果一目了然～
%%%
%%%%
%%%2020.6.5%
%%%%
%%%%
%%%%
%%%%
%%復化梯形公式求二重積分%%%
%%%%
%%%%
%%%%
%%%% liu deping %%%
%%%%
%%clear all;%%
%被積函式及積分上下限匯入；
s=input
('請輸入函式表示式：f = '
,'s');
f=inline
(s);
a =input
('請輸入積分變數x左邊界a的值:');
b =input
('請輸入積分變數x右邊界b的值:');
c =input
('請輸入積分變數y左邊界c的值:');
d =input
('請輸入積分變數y右邊界d的值:');
h1 =
input
('請輸入積分變數x步長h1的值:');
h2 =
input
('請輸入積分變數y步長h2的值:');
m=round
((b-a)
/h1);n=
round
((d-c)
/h2);%
%%係數矩陣t，t[i,j]為復化梯形公式的係數；
tt=zeros
(m+1
,n+1);
tt(1,
1)=1
;tt(m+1,1
)=1;
tt(1,n+1)
=1;tt
(m+1
,n+1)=
1;for i=2:m
tt(i,1)=
2;tt(i,n+1)
=2;end
for j=2:n
tt(1
,j)=2;
tt(m+
1,j)=2
;end
tt(tt==0)
=4%%係數表 %
%%計算各二維節點函式值，並存放於矩陣f中；
f=zeros
(m+1
,n+1);
for i=0:m
for j=0:n
f(i+1
,j+1)=
f(a+i*h1,c+j*h2)
;
endend %
%%結果輸出；
format long
;fprintf
('各節點函式值；')f
fprintf
('復化梯形公式計算結果；'
)tnm=
(b-a)
*(d-c)/(
4*m*n)
*sum
(sum
(tt.*f)
)
復化辛普森公式求二重積分matlab原始碼及例題
2020.6.5 復化simpson公式求二重積分 liu deping clear all 被積函式及積分上下限匯入 s input 請輸入函式表示式 f s f inline s a input 請輸入積分變數x左邊界a的值 b input 請輸入積分變數x右邊界b的值 c input 請輸入積...
matlab二重定積分 二重積分 matlab
第六章 用matlab 計算二重積分 由於二重積分可以化成二次積分來進行計算，因此只要確定出幾分區域，就可以反覆 使用int 命令來計算二重積分。例6.4.1 計算二重積分yd ixedxdy d是由直線 x 0,y 1,y x 所圍區域 解該積分可以寫成yy idyxe dx或yy idxxe d...
用python求一重積分和二重積分
首先是對一元函式求積分，使用scipy下的integrate函式 from scipy import integrate def g x return 1 x 2 0.5 用integrate下的quad函式可以同時求出積分結果和誤差 res,err integrate.quad g,1,1 1和1...