计算中间数据以字母代替后，后续计算仅显示代数式结果，无数值结构问题

[hahaha]

我计算一个 应力张量的主应力 即 三阶矩阵的特征值，方法大概是，求三个应力不变量 然后以这三个不变量为系数 构造三次方程 求解
这里的应力不变量我用中间变量 I1 I2 I3表示，求解方程，结果是以 I1 I2 I3 表示的，无法以数值显示。下面是我的实验代码，
clear all
clc
s=eye(3);
%计算应力不变量
I1=trace(s);%计算第一应力不变量I1
I2=det(s(2:3,2:3))+det(s(1:2:3,1:2:3))+det(s(1:2,1:2));%计算I2
I3=det(s);%计算I3
解三次方程，求主应力。
D=solve('S^3-I1*S^2+I2*S-I3')

我的目的就是想得出具体的数值，而没有用I1 I2 I3 表示的代数式。请大虾知道，3q。

PS.我知道有专门的求特征值得函数（eig）。提这个问题就是想 解决 sym 转化为 数值的问题，我在帮助中看到了double（）函数，但是老是出错。如下:

double(D(1))
??? Error using ==> mupadmex
Error in MuPAD command: DOUBLE cannot convert the input expression into a double
array.

If the input expression contains a symbolic variable, use the VPA function instead.

Error in ==> sym.sym>sym.double at 936
Xstr = mupadmex('symobj::double', S.s, 0);

[mathjiang]

表达式
D=solve('S^3-I1*S^2+I2*S-I3')
中，原I1，I2，I3的值并没有传入方程
S^3-I1*S^2+I2*S-I3=0
因为solve是符号求解器。
破解之法是求数值解：
clear all
clc
s=eye(3);
%计算应力不变量
I1=trace(s);%计算第一应力不变量I1
I2=det(s(2:3,2:3))+det(s(1:2:3,1:2:3))+det(s(1:2,1:2));%计算I2
I3=det(s);%计算I3
%解三次方程，求主应力。
c = [1 -I1 I2 -I3];%这就是三次多项式 S^3-I1*S^2+I2*S-I3
D=roots(c);

[hahaha]

非常感谢mathjiang的解答，醍醐灌顶！对有些东西还是认识太浅，还需加强学习。