[Débutant] Utiliser une variable au début du code sachant que le but est de trouver la valeur de celle-ci!

Bonjour à tous et à toutes ,
Je sollicite vos compétences pour m'aider à résoudre un problème que j'ai eu depuis quelques jours !! en effet je devrais calculer la température d'un disque tournant dans une boite (température T1) rempli d'un mélange (air,hélium) à l'aide d'un code Matlab(j'ai utilisé un schéma de résistance thermique), le problème c'est que la densité du mélange dépend de la température dans l'enceinte . pourriez-vous m'aider ?
Merci d'avance!
Voilà le code:

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
  clear all
  close all
  T_air=288;                                        % température de l'air à 15°
  p=101325;                                         % pression en pascal
  rho=0.169;                                        % la densité de l'air
  mu=1.9388e-5 ;                                     % la viscosité dynamique de l'air
  % M=[0.029;0.004];                                  % Masse molaire respectivement de l'air et de l'hélium
  %mu=[1.802e-5;1.9388e-5];                          % viscosité dynamique respectivement de l'air et de l'hélium
  Ma=0.029;                                                     %Masse molaire de l'air
  Mhe=0.004;                                              % Masse molaire de l'hélium
  R=0.25;                                          % rayon total du disque en m
  N=10000;                                          % nombre de tours (en rpm)
  omega=2*pi*N/60;
  e= 0.010 ;                                        % epaisseur du disque  en m
  %delta=[0.01:0.005:0.03];                             %vecteur air gap
  delta=0.01;                                       % air gap
  L=0.01+2*R+0.01;                                  % Longueur de la boîte
  l=e+2*delta;                                      % largeur de la boîte
  E=2*R+0.02;                                           % Epaisseur de l'enceinte
  G=delta/R;                                        % air gap ratio
  k=0.02476;                                        % Conductivité thermique de l'air à 15°
  P_bearing=0.5*omega*0.05*(295/2)*0.015;           %Pertes produites par les roulements
  Q=P_bearing+900;                                  % Chaleur dégagée par le moteur+roulements
 
 
 
  %Introduire la densité du mélange
Xa = input('  \n entrer la fraction de l''air : \n')
 Xhe = input('\n entrer la fraction de l''hélium : \n ')
 
if Xa+Xhe==1 
    rho_m=(Xa*Ma+Xhe*Mhe)*p/(8.314*T1)     % densité du mélange(air+hélium)
else 
    input('erreur')
 
end
X=[Xa;Xhe];
%Re_theta=rho_m*omega*R^2/mu;
g1= 0.2091;
g2=1.2139;
%%%% Construction de la matrice %%%% (Appelée C au lieu de G!!!!!)
 for i=1 : 2
    for j=1 :2
        C=zeros;
        C(1,1)=(1/g1)+(1/g2);
        C(2,2)=-C(1,1);
    end
end
G1=inv(C);
T1=Q/C(1,1)+T_air
T3= T1-Ra*Q