Hello all I attached a new function "ffreson.sci" with my last mail; Change at line 89 of function "ffreson.sci" [quo,rest] = pdiv(der,z*z-1) by [rest,quo] = pdiv(der,z*z-1). Excuse me Yours sincerely lucien.povy at free.fr