[Scilab-users] Vectorization issue

[Carrico, Paul]

Dear
I failed in using vectorization in the example immediately bellow; but I do not remember if it's possible : any advice ?
Thanks
Paul
######################################################################################

mode(0)



n = 20

Km = 0;



b = 100;

d = 20;

nu = 0.3;



// using loops

for i = 0 : n

    for j = 1 : n

        Km = Km + (144*(b/d)^4)/(  %pi^6*(2*i+1)^2*j^2*((2*i+1)^2*(b/2*d)^2+j^2)  );

    end

end



Km = 1.2 + (nu/(1+nu))*Km

k = 1/Km



// using vectorization

a = [0:n]';

b = ones(n,1);

i = a.*.b;

j = b.*.a;



Km2 = zeros(n*(n+1),1);

//toto = i. *j

Km2 = (  %pi.^6*((2*i+1).^2. * j.^2)*((2*. i+1).^2*(b/2*d).^2 + j^2)  )

Km2 = (144*(b/d).^4)/ Km2

Km2 = 1.2 + (nu/(1+nu))*Km2

k2 = 1/Km2




EXPORT CONTROL :
Cet email ne contient pas de données techniques
This email does not contain technical data

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20160526/4cf3cfdf/attachment.htm>