求解三元一次方程时，​得出的解是一个代数式​而不是数。

deltax=200;

deltay=-300;

deltaz=50;

taoxy=500;

taoyz=0;

taozx=300;%x,y,z方向的正应力和切应力

I1=deltax+deltay+deltaz;

I2=deltax*deltay+deltay*deltaz+deltaz*deltax-taoxy^2-taoyz^2-taozx^2;

I3=deltax*deltay*deltaz-deltaz*taoxy^2-deltax*taoyz^2-deltaz*taozx^2+2*taoxy*taoyz*taozx;%I1，I2，I3为有正应力和切应力组成的系数

syms x;

x = solve(x^3-I1*x^2+I2*x-I3 == 0,'x'),

x1 = 50;

x2 = 50*161^(1/2) - 50;

x3 = - 50*161^(1/2) - 50;%x1,x2,x3为上述方程的解

a1= deltax-x1;a2= deltay-x1;a3=deltaz-x1;%a1,a2,a3上述方程的解组成的系数

[l1,m1,n1]=solve('a1*l1+taoxy*m1+taozx*n1=0','taoxy*l1+a2*m1+taoyz*n1=0','l1^2+m1^2+n1^2=1','l1','m1','n1')%求解l1,m1,n1

Answer 1

[lefow]

最后加上：

l1=subs(l1)

m1=subs(m1)

n1=subs(n1)