matlab积分怎么积分？

Question

i=1:0.01:pi;
x=1:0.01:pi;
N1=13;

N2=15;

for i=1:pi+1
if 0.4*pi+1<=x(i)+1 && x(i)+1<=0.6*pi+1
phi_2(i)=(N2-N1)*x(i)+1/2;
end
if 0.0*pi+1<=x(i)+1 && x(i)+1<=0.2*pi+1
phi_2(i)=((N2-N1)*x(i)+1+pi)*0.5;

end
end

syms i;
Er=1-abs(exp(1j*(phi_2(i)-((N1+N2)/2)*i)).*cos(((N2-N1)/2)*i-phi_2(i) )).^2;

s=int((Er),0.4*pi,0.6*pi);

Answer 1

f=@(x);%积分函数
y0=f(0);%定积分初值
y1=f(2);%定积分终值
s=0;%积分结果
n=1;%将定积分区间分成了n份
while abs(s-pi)>=0.0001%确定精度为万分之一
n=n+1;%精度不满足的话份数就加1
s=0;%每次判断精度满足与否进行下一次循环，设定积分结果的初始值为0
h=2/n;%分成n份后的步长
for x=h:h:2-h%积分循环
s=s+h*f(x);
end
s=2*(s+(y0+y1)*h/2);%积分结果
end
disp('梯形公式积分得到圆的面积为'),s