Hello!

I have this scilab code to compute MSD, and I want it to be an executable program:

`//open file`

fileID = mopen('file001.txt','r');

sample = mfscanf(-1,fileID,'%i');

//initialize variables and arrays

nData = size(sample,1);

numberOfDeltaT = floor(nData/4);

msd = zeros(numberOfDeltaT,2);

deltaCoords = cell(nData,1);

//compute MSD then store in array

for dt = 1:numberOfDeltaT

deltaCoords = sample(1+dt:$) - sample(1:$-dt);

squaredDisplacement = sum(deltaCoords.^2,2);

`msd(dt,1) = mean(squaredDisplacement);`

msd(dt,2) = length(squaredDisplacement);

end

//close file

mclose(fileID);

//plot MSD

plot(msd(:,2),msd(:,1));

//end

I want the input to be the file to be read, and will produce the MSD plot.

How do I do this? Should I start with `function = calc_MSD(inputFile)`

?

I also want the `for`

loop to be using `parfor`

. Is this possible in scilab?

Thank you!

Scilab Not Spam Edit Delete

Login to add comment