回复武汉某大学老师的一个问题：北太天元求根和范德蒙行列式

%北太天元求行列式给出的多项式的根

% det( [1 1 1; 2 4 x; 3 9 x^2 ] ) = 0 的根

% 实际上有范德蒙行列式的知识，我们知道 x = 2 或者 3

% 但是我们假设没有学过范德蒙矩阵的知识，

A = [ 1 1 ;

         2 3 ;

         4 9 ; ]

   c = zeros(3,1);

   for j=1:3

      if j==1

         ind = [2,3];

      elseif j==2

               ind = [1, 3];

      else

                ind = [1,2];

      end

      c(j) = (-1)^(j+3) * det( A( ind,:));

   end

   %p(x) = c(1) +c(2)*x+c(3)*x^2

  % 但是使用roots的时候要求输入的系数是按降幂排列

   c = c(end:-1:1);

   [r,a] = roots(c)

   % roots 是 a 的特征值

   [v,d] = eig(a);

   xx = diag(d)