Matematyka.pl

Kod: Zaznacz cały

Ferrari := proc(P,x)
   s:=-1/4*coeff(P,x,3)/coeff(P,x,4);
   Q := subs(x=y+s,1/coeff(P,x,4)*P);
   p := coeff(Q,y,2);
   q := coeff(Q,y,1);
   r := coeff(Q,y,0);
   R := solve(8*z^3+8*p*z^2+(2*p^2-8*r)*z-q^2=0,z);
   y0 := R[1];
   if y0 = 0 then
      y0 := R[2];
   fi;
   if y0 = 0 then
      y0 := R[3];
   fi;      
   A := x^2-sqrt(2*y0)*x+p/2+y0+q*sqrt(2*y0)/(4*y0);
   B := x^2+sqrt(2*y0)*x+p/2+y0-q*sqrt(2*y0)/(4*y0);
   DA := coeff(A,x,1)^2-4*coeff(A,x,2)*coeff(A,x,0);
   DB := coeff(B,x,1)^2-4*coeff(B,x,2)*coeff(B,x,0);
   Y := [(-coeff(A,x,1)-sqrt(DA))/(2*coeff(A,x,2)),(-coeff(A,x,1)+sqrt(DA))/(2*coeff(A,x,2)),(-coeff(B,x,1)-sqrt(DB))/(2*coeff(B,x,2)),(-coeff(B,x,1)+sqrt(DB))/(2*coeff(B,x,2))];
   X := [s+Y[1],s+Y[2],s+Y[3],s+Y[4]];
   seq(k,k=X);
 end;