[幸儿李]

各位大虾，利用下面的程序求出矩阵c之后，若想通过“c（2，1）=0”来求解W（也就是系统的各阶固有频率），该采用什么样的方法？新手上路，请大虾们指教一二
% function f=cdjz(n)
w=sym('w');
c=ones(2);
% f1=sym('[x;m]');
j=[46474 879.69 462.3 445.15 528.63 193.45 5.376 12.124 2.912 11.611 4.935];
k=[0 36.186 0.0297 4.1489 22.21 0.2329 19.529 0.1338 0.0117 0.0146 0.1717];
for i=2:11
c=c*[1 1/k(i);-j(i)*w^2 1-j(i)*w^2/k(i)];
end
% f=c*f1;