求助关于matlab求积分上限值

[kuangdang178]

已知其他值求积分上限，上网查了好多方法，到现在弄出一个，但是总出错，求问。
theta0=0;
thetamax=pi/3;
I = 1e8;%中心光强，单位值1e8
m=20;
n=20;
fluxtotal = 0.25*I*pi*0.75;%1/4总光通量 余弦辐射体 0-pi/3
flux0=fluxtotal/(m*n); %分割每一份的光通量
theta=ones(n,1)'; %theta角分割矩阵
theta(1)=theta0; %theta初始值
f=@(x) I*cos(x)*sin(x);
for i=2:1:n
theta(i)=solve((pi/(2*m))*integral(f(theta(i)),theta(i-1),theta(i))-flux0);
end