matlab求根的原程序,MATLAB求根程序求帮助

我有一个函数g=m*(m1/m)*(m2/m)*(r^2)*rm*Wc1*cos(2*u),我需要一个程序求这个函数在在(-2*pi,2*pi)上所有的根。

其中u是自变量，其他所有参数都是已知。已知参数如下：

kx=10037;ky=10037;kP=1500;kD=1500;kE=1500;

f1=0;f2=0;fx=125;fy=125;fP=15;fD=15;fE=15;

m=100;m1=4;m2=4;M=m+m1+m2;

r=0.05;

JP=10;JD=10;JE=10;

rm=0.04;

lz=0.3;ly=0.3;

Omiga=157;

kt=0;

Wc1=abs(((Omiga^4)*(M^2))/((((Omiga^2)*M)^2)+(Omiga^2)*(fx^2)))...

*cos(atan(fx*Omiga/((Omiga^2)*M-kx)))+...

abs(((Omiga^4)*(M^2))/((((Omiga^2)*M)^2)+(Omiga^2)*(fy^2)))*cos(atan(fy*Omiga/((Omiga^2)*M-ky)))-...

abs(((Omiga^4)*(JP^2))/((((Omiga^2)*JP)^2)+(Omiga^2)*(fP^2)))*((lz^2)*M/JP)*cos(atan(fP*Omiga/((Omiga^2)*JP-kP)))-...

abs(((Omiga^4)*(JD^2))/((((Omiga^2)*JD)^2)+(Omiga^2)*(fD^2)))*((lz^2)*M/JD)*cos(atan(fD*Omiga/((Omiga^2)*JD-kD)))+...

abs(((Omiga^4)*(JE^2))/((((Omiga^2)*JE)^2)+(Omiga^2)*(fE^2)))*((ly^2)*M/JE)*cos(atan(fE*Omiga/((Omiga^2)*JE-kE)));

求大神帮我写一个程序求这个函数在在(-2*pi,2*pi)区间内所有的根