clc; %清屏
clear; %清除变量
close all; %关闭
syms x;
f(x) = x;
df = int(f(x),x);%求不定积分
figure;
subplot(1,2,1);
ezplot(df);
%{
int(F,var,a,b) 函数F对变量var求积分，下标为a,上标为b,
不要a，b表示不定积分
%}
%二重积分
df2 = int(int(f(x),x),x);
subplot(1,2,2);
ezplot(df2);
原文：http://www.cnblogs.com/TheoryDance/p/6135586.html