Bonjour,

Je suis débutante en matlab. ça me parait difficile. Je veux résoudre ce système d’équations

Code:
  P= exp(-r.^2.*P+R.^2.*Q).*sum((r.^2.*P +R.^2.*Q).^m./factorial(m)) // with m=(0:1:k-1); ;k=5

  Q=exp(-r.^2.*P + R.^2.*Q.*sum((r.^2*P + R.^2.*Q).^n./factorial(n)) // n=(k:1:K_avg); K_avg=10
J'ai essayé de le faire avec "solve" Résultat d'exécution ------> error eqns
avez vous une idée ?

Code:
th = -80:-3:-106;
list1 = zeros(size(th));
list2 = zeros(size(th));
r=[ 15.8 20.8 27.5 36.3 47.6 63.9 83.1 109.6];
R=[ 8.12 10.9 14.4 19.2 25.6 33.1 43.5 75.6];

for i = 1:numel(th)
 K_avg=10;
 k=5;
 m=(0:1:k-1);
n=(k:1:K_avg);
syms P Q
eqns= [exp(-r.^2.*P+R.^2.*Q).*sum((r.^2.*P +R.^2.*Q).^m./factorial(m))-P==0,exp(-r.^2.*P + R.^2.*Q.*sum((r.^2*P + R.^2.*Q).^n./factorial(n))-Q==0)];%---> error eqns
vars = [P Q];
[sol1,sol2]=solve(eqns,vars);
list1(i) = sol1;
list2(i) = sol2;
end