有没有matlab?
两行就能把根解出来:
>> syms x y
>> solve('x^3+x-y')
ans =
[ 1/6*(108*y+12*(12+81*y^2)^(1/2))^(1/3)-2/(108*y+12*(12+81*y^2)^(1/2))^(1/3)]
[ -1/12*(108*y+12*(12+81*y^2)^(1/2))^(1/3)+1/(108*y+12*(12+81*y^2)^(1/2))^(1/3)+1/2*i*3^(1/2)*(1/6*(108*y+12*(12+81*y^2)^(1/2))^(1/3)+2/(108*y+12*(12+81*y^2)^(1/2))^(1/3))]
[ -1/12*(108*y+12*(12+81*y^2)^(1/2))^(1/3)+1/(108*y+12*(12+81*y^2)^(1/2))^(1/3)-1/2*i*3^(1/2)*(1/6*(108*y+12*(12+81*y^2)^(1/2))^(1/3)+2/(108*y+12*(12+81*y^2)^(1/2))^(1/3))]
>> syms x y
>> solve('x^3+x-y')
ans =
[ 1/6*(108*y+12*(12+81*y^2)^(1/2))^(1/3)-2/(108*y+12*(12+81*y^2)^(1/2))^(1/3)]
[ -1/12*(108*y+12*(12+81*y^2)^(1/2))^(1/3)+1/(108*y+12*(12+81*y^2)^(1/2))^(1/3)+1/2*i*3^(1/2)*(1/6*(108*y+12*(12+81*y^2)^(1/2))^(1/3)+2/(108*y+12*(12+81*y^2)^(1/2))^(1/3))]
[ -1/12*(108*y+12*(12+81*y^2)^(1/2))^(1/3)+1/(108*y+12*(12+81*y^2)^(1/2))^(1/3)-1/2*i*3^(1/2)*(1/6*(108*y+12*(12+81*y^2)^(1/2))^(1/3)+2/(108*y+12*(12+81*y^2)^(1/2))^(1/3))]