Bonjour à tous,

Je recherche les racines d'un polynômes fractionnels de la forme:
sum(lambda./(lambda-mu)) où mu est l'inconnu et lambda un vecteur de taille 100.

J'ai donc pensé à la fonction solve et je l'ai utilisée de la manière suivante:

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
syms mu real;
mu=double(solve(sum(lambda./(lambda-mu))/M-1/c==0,mu,'Real', true));
Elle fonctionne bien mais prend énormément de temps, or je dois trouver les racines de ce polynôme avec plusieurs expériences aléatoires différentes.

Connaissez une méthode similaire mais plus rapide? J'ai pensé à fzeros() ou d'autres dans le même style, mais je ne peux pas savoir le voisinage des racines.

Je vous remercie d'avance!