[TTT_IOU]

试试这个程序：
syms m x
w=[pi/2,pi,pi*3/2];
N=[pi/2-1,-2,-1.5*pi-1];
for i=1:length(w)
f=int(x^m*cos(x),x,0,w(i));
ff(i)=solve(f-N(i));
end
ff=double(ff);
望有帮助！！！