复化梯形公式求二重积分matlab源码

复化梯形公式求二重积分matlab源码

这段代码具有很好的交互性和通用性，将代码复制到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);    
    end
end 
 
%%%结果输出；
format long;
fprintf('各节点函数值；')
F
fprintf('复化梯形公式计算结果；')
Tnm=(b-a)*(d-c)/(4*m*n)*sum(sum(TT.*F))

之后还有辛普森公式求二重积分，希望喜欢的小伙伴持续关注篱落～～成殇～～