à l'aide je suis nulle en programmation

b = 3.0e-10;
e_f = 0.3;
y = 0.1;

r_m0 = 10E10 ;
r_f0 = 10E-11;
h = 0.001;
e0 = 0;
ef = 1;

[e,r_m,r_f] = syskutta('rho_prime',e0, r_m0, r_f0, ef, h);

figure(1);
plot(e,r_m,e,r_f); 

% Initialisation des vecteurs permettant la sauvegarde des valeurs des résultats C1, C2, C3 et C4
C1 = zeros(length(r_m),1);
C2 = zeros(length(r_m),1);
C3 = zeros(length(r_m),1);
C4 = zeros(length(r_m),1);

% Boucle de calculs de C1, C2, C3 et C4
for ii = 1:length(r_m):
A = [ 1  0  0  -b^2*(r_f(ii) );
1  -b^2*(r_m(ii))  -b*(sqrt(r_f(ii))) 0;
0  0  0  1;
0  r_m(ii)  (-y/b)*sqrt(r_f(ii))  0];

B=[0;0;1/e_f;0]; 

[C1(ii),C2(ii),C3(ii),C4(ii)] = coefficientsC_i(A,B); 
end;


figure(2);
plot(e,C1, e,C2, e,C3,e,C4); 

function[c_1,c_2,c_3,c_4] = coefficientsC_i(A,B)

X = B\A; 
end;