Bub in module CACSD macro freson
Reported by Lucien Povy
BUG DESCRIPTION:Bug in cacsd module, macro "freson.sci", for some continuous system freson returns [].
----------------
ERROR LOG:
----------
HOW TO REPRODUCE THE BUG: Do :
s=%s, h=syslin("c",4+4.6188*s,s+1.333*s*s+0.3333*s^3)
fre=freson(h)
This version of freson returns [] Scilab-6..
The good value is fre=0.4825 Hz in scilab-5... and by black and datatips
-------------------------
OTHER INFORMATION:
The problem is :
At lines 47 and 54 we do hh=h*horner(h,-%s) or hh=h*horner(h,1/%z),
the result must be even rationals, but in the exemple it is not :
One solution is : replace hh=clean(h*horner(h,-%s)) or hh=clean(h*horner(h,1/%z))
------------------