求解Aa^6+Ba^4+Ca^2=D振动方程中a与omega的幅频特性曲线，c=1;m=5;k2=1;K=2;n=0.4;k1=2;F=2;其中A=(c/(8m))^2+(3k2/(8momega))^2;
B=-c/(4m)(c/(2m)+Komega^(n-1)sin(npi/2)/(2m))+3k2/(4momega)((k1^2-momega^2)/(2momega)+Komega^(n-1)cos(npi/2)/(2m));
C=(c/(2m)+Komega^(n-1)sin(npi/2)/(2m))^2+((k1^2-momega^2)/(2momega)+Komega^(n-1)cos(npi/2)/(2m))^2;
D=(F/(2m*omega))^2;
程序编写后计算时间长一直未出结果，望解答，不胜感激
程序：
c=1;m=5;k2=1;K=2;n=0.4;k1=2;F=2;
syms omega;
A=(c/(8m))^2+(3k2/(8momega))^2;
B=-c/(4m)(c/(2m)+Komega^(n-1)sin(npi/2)/(2m))+3k2/(4momega)((k1^2-momega^2)/(2momega)+Komega^(n-1)cos(npi/2)/(2m));
C=(c/(2m)+Komega^(n-1)sin(npi/2)/(2m))^2+((k1^2-momega^2)/(2momega)+Komega^(n-1)cos(npi/2)/(2m))^2;
D=(F/(2momega))^2;
h=ezplot(@(a,omega) (Aa.^6+Ba.^4+C*a.^2-D),[-.5,3,0,3])
str1=get(get(gca,'Title'),'String')
set(h,'Color',[1 0 0],'LineStyle',':')