Question: Simplifier une exprission assez grande

Bonjour,

L'expression

L'exprission kw ci-dessus est une procédure en polynôme Pol de variable "p". Est-il possible de simplifier la valeur k:=kw(Pol)? Sur Mon portable, à certain moment de calcul, je vois que maple est presque bloquer.

#=============================================================:

> restart:
> with(LinearAlgebra):
> kw:=proc(Pol)    
> local a0,a1,a2,a3,R,alpha0,alpha1,alpha2,k:
> a3:=coeff(Pol,p,0):a2:=coeff(Pol,p,1):a1:=coeff(Pol,p,2):a0:=coeff(Pol,p,3):
> R:=Determinant(Matrix([[a0,a1,a2,a3,0],[0,a0,a1,a2,a3],[3*a0,2*a1,a2,0,0],[0,3*a0,2*a1,a2,0],[0,0,3*a0,2*a1,a2]]));
> alpha0:=[diff(a0,y),diff(a0,x)+diff(a1,y),diff(a1,x)+diff(a2,y),diff(a2,x)+diff(a3,y),diff(a3,x)];
> alpha1:=Determinant(Matrix([alpha0,[a0,a1,a2,a3,0],[-a0,0,a2,2*a3,0],[0,-2*a0,-a1,0,a3],[0,0,-3*a0,-2*a1,-a2]]));
> alpha2:=Determinant(Matrix([alpha0,[0,a0,a1,a2,a3],[-a0,0,a2,2*a3,0],[0,-2*a0,-a1,0,a3],[0,0,-3*a0,-2*a1,-a2]]));
> k:=factor((simplify((diff(alpha2/R,x)+diff(alpha1/R,y))))):
> #k:=collect((diff(alpha2,x)+diff(alpha1,y))*R-(alpha2*diff(R,x)+alpha1*diff(R,y))):
> end proc:
> Pol:=(-b1*x^2+c2*x*y+c1*y*x^2+a0*x^4+x*a3+c0*y*x^3+a2*x^2+a1*x^3-b0*x^3-b2*x+y*c3-b3)*p^3+(3*b0*y*x^2-2*c1*x*y^2-alpha0*x^3-2*a1*x^2*y-3*a0*y*x^3-a2*x*y-3*x^2*c0*y^2-alpha1*x^2+2*b1*x*y+beta0*x^2-c2*y^2+b2*y+beta1*x)*p^2+(-b1*y^2-3*x*b0*y^2+3*x*c0*y^3-1+c1*y^3-beta1*y+3*x^2*a0*y^2+a1*x*y^2+alpha1*x*y+2*x^2*alpha0*y-2*x*beta0*y)*p-y*(c0*y^3-b0*y^2+x*a0*y^2+x*alpha0*y-beta0*y+phi):
> k:=kw(Pol):

#=============================================================:

 

Merci d'avance,

Gérard.

Please Wait...