数学实验:通过Matlab软件操作，对函数性质有进一步了解 f(x)=sin^2(1.006)x

f = @(x)(sin(1.006*x).^2);
format long;
m = [100,200,1000];
for i=1:3
    %区间分n段等长子区间
    n = m(i);
    xn = (0:n-1)/n*pi;   %xn(i)为第i段区间左端点值
    d = pi/n;       %子区间宽度
    xnm = xn+d/2;   %区间中点
    ynm = f(xnm);   %区间中点的函数值
    S(i) = sum(ynm); %求和
end
%求平均值
Smean = sum(S)/sum(m) %函数均值

%解方程
x = fsolve(@(x) sin(1.006*x)^2-Smean,0.1)

%f(x)=sin^2(m/100)x，x属于[0,pi],m=403
%(1)-----------------
x=0:0.01:pi;
m=403;
f = sin(m*x/100).^2;
plot(x,f);
%(2)---------------
n=100;
x1 = 0:pi/n:pi;
x1 = x1 + pi/2*n;
x1 = x1(1:n-1);
f1 = sin(m*x/100).^2;
n=200;
x2 = 0:pi/n:pi;
x2 = x2 + pi/2*n;
x2 = x2(1:n-1);
f2 = sin(m*x/100).^2;
n=1000;
x3 = 0:pi/n:pi;
x3 = x3 + pi/2*n;
x3 = x3(1:n-1);
f3 = sin(m*x/100).^2;
%(3)--------------------------
m1 = (99*mean(f1)+199*mean(f2)+999*mean(f3))/(99+199+999);
disp('由2产生的函数值的平均值')
disp(m1)
%(4)-----------------------
xx = fsolve(@(x) sin(m*x/100)^2-m1,0.1);
disp('对于x解方程f(x)=均值')
disp(xx)

请问有答案了吗