x=[0:0.01:2 2+eps:0.01:4,4];%一行多列矩阵
y=exp(x.^2).*(x<=2)+80./(4-sin(16*pi*x)).*(x>2);
y(end)=0;
x=[eps,x];
y=[0,y];

%fill(x,y,c) x和y是要画的点，可以同为向量或者同为矩阵。

%参数c是一个颜色向量，可用符号表示，如你给的'r',也可用[r,g,b]表示其中r,g,b的值都在0到1之间，比如说t=0:0.05:2*pi，%x=sin(t)，y=cos(t)，fill(x,y,'b')。

fill(x,y,'b')

%[email protected](x,y)exp(-x.^2/3).*sin(x.^2+2*y)
%其中@(x,y)exp(-x.^2/3).*sin(x.^2+2*y)为匿名函数，第一个括号里面是自变量，第二个括号里面是表达式，@是函数指针
%[email protected](x,y)exp(-x.^2/3).*sin(x.^2+2*y)表示将匿名函数@(x,y)exp(-x.^2/3).*sin(x.^2+2*y)赋值给f，于是f就表示该函数。
[email protected](x)exp(x.^2)+80./(4-sin(16*pi*x)).*(x>2);
l=cintegral(f,0,4)   %函数0~4的定积分