I have a piece of code which creates a 3 x 46 array
theta = 0;phi = pi/4;r = 1 : 0.2 : 10;theta = repmat(theta,size(r));phi = repmat(phi, size(r)); %obtain position vectors for the points along a radial line where the velocities are evaluated%
[x,y,z] = sph2cart(theta,phi,r);p = [x; y; z];
I then have a loop where you evaluate a function of the distance from each point represented in this array to a set of other points forming a 3 x 100 array:
for II=1:46 for IJ=1:100 normal=(p(1:3,II)/(norm(p(1:3,II))))'; r=(p(1:3,II)-sing(1:3,IJ))'; A((1+(II-1)*3):(3+(II-1)*3),(1+(IJ-1)*3):(3+(IJ-1)*3))=gradletSlip(r,in)'; endend
However, something seems to be wrong as when I loop through IJ the number changes as it is a different distance each time, but when II changes from 1 to 2, it is supposed to be at another point in the 3 x 46 array and yet the normal is not changing, but this is just the radial vector, so it seems like it is not moving to the new point in the 3 x 46 array as it does the loop. Any idea how to fix this?
Best Answer