Page 1 sur 1

csolve et conj

Publié : sam. nov. 02, 2013 8:54 pm
par alb
csolve(z^2-conj(z)=-1/4,z) oublie la solution z=1/2
mais je pense que le pb a deja ete evoque ailleurs:
z:=a+i*b;
Z:=z^2-conj(z)+1/4;
solve([re(Z),im(Z)],[a,b]) renvoie [[-1/2,1],[-1/2,-1]]

Re: csolve et conj

Publié : dim. nov. 03, 2013 8:16 am
par parisse
Il y a en effet une simplification qui ne tient pas compte d'un cas particulier.

Code : Tout sélectionner

diff solve.cc solve.cc~
5101c5101
< 	  if (is_linear_wrt(eq[i],var[j],a,b,contextptr) && is_zero(derive(eq[i],var,contextptr)) && !is_zero(simplify(a,contextptr))){
---
> 	  if (is_linear_wrt(eq[i],var[j],a,b,contextptr) && !is_zero(simplify(a,contextptr))){