3-28.m

来自「《MATLAB 7.0编程基础》第3章 （源码实例）主要讲解MATLAB的数值计」· M 代码 · 共 14 行

f = @(x,y)y*sin(x)+x*cos(y);
% 编写匿名函数，将句柄赋给f
S = dblquad(f,pi,2*pi,0,pi)
% 计算二重积分

% 符号计算能得到高精度的解，语法请查阅下章
syms x y;
% 定义符号变量
S1 = int(int(y*sin(x)+x*cos(y),x,pi,2*pi),y,0,pi)
% 计算解析解

S2 = vpa(S1)
% 得高精度解