fricas (1) -> digits 20
Type: PositiveInteger?
fricas -- n:=x^3+a1*x^2+a2*x+a3 ::Polynomial Fraction Integer
Type: Polynomial(Fraction(Integer))
fricas R:=(9*a1*a2-27*a3-2*a1^3)/54
Type: Polynomial(Fraction(Integer))
fricas S:=(R+(Q^3+R^2)^(1/2))^(1/3)
Type: Expression(Integer)
fricas T:=(R-(Q^3+R^2)^(1/2))^(1/3)
Type: Expression(Integer)
fricas x1:=S+T-a1/3
Type: Expression(Integer)
fricas x2:=-(S+T)/2-a1/3 + %i*sqrt(3)*(S-T)/2
Type: Expression(Complex(Integer))
fricas x3:=-(S+T)/2-a1/3 - %i*sqrt(3)*(S-T)/2
Type: Expression(Complex(Integer))
fricas a5:=x^3+a1*x^2+a2*x+a3 ::Polynomial AlgebraicNumber
Type: Polynomial(AlgebraicNumber?)
fricas a6:=(x-x11) ::Polynomial AlgebraicNumber; Type: Polynomial(AlgebraicNumber?)
fricas a7:=monicDivide(a5, fricas a77:=a7.quotient; Type: Polynomial(AlgebraicNumber?)
fricas a78:=a7.remainder; Type: Polynomial(AlgebraicNumber?)
fricas qu1 :=eval(a77,
Type: Expression(Integer)
fricas rem1:=eval(a78,
Type: Expression(Integer)
fricas eval(rem1,
Type: Expression(Float)
... --wyscc, Mon, 14 Nov 2005 02:53:03 -0600 reply How about this:
fricas pkg:= SOLVEFOR(UP('x, |