Writing a matlab program to calculate Diffusion coefficients using the Sigmund correlation

I am currently trying to write a matlab program to solve the generalised sigmund correlation in order to find binary diffusion coefficients. Here is an exert from the paper that I am using to guide me:

The equation for D_ij (diffusion coefficient) itself seems relatively straightforward,

however the parameter rho_mc is causing me trouble as I am not sure how to do the sum of a series with multiple input variables such as this one. Below is what i have got so far:

%Calculating diffusion coefficients using Sigmunds Correlation%

%Use the prompt function to specify the variables%
prompt1 = ('Value of Rho_0m');

rho_0m = input(prompt1);

prompt2 = ('Value of Rho_m');

rho_m = input(prompt2);

prompt3 = ('Value of D_0ij');

D_0ij = input(prompt3);

promptn = ('n value?');

%For a specified vale of n, run file rho_mc.m to calculate rho_pr%

for n= input(promptn);
    run rho_mc
    rho_pr = (rho_m/rho_mc);
end

%Specifying the correlation to calculate the diffusion coefficient using
%all the variables specified and calculated
D_ij = ((rho_0m*D_0ij)/rho_m)*(0.99586+(0.096016*rho_pr)+(0.22035*(rho_pr^2))+(0.032874*(rho_pr^3)))

and also the incomplete sub-program for calculating rho_mc:

syms (Z,V);
X= Z*V;

S1=symsum(V^(2/3),i,1,n);
S2=symsum(V^(5/3),i,1,n);

rho_mc =(S1/S2);

any direction that anyone could give me would be massively appreciated!

Thanks again

Answer 1

[Mina]

program for coffecient of diffusion of boron against pressure

[Walter Roberson]

I do not understand how writing such a program is a solution to the questions that were asked here?