Em 16 de abril de 2010 Allin <cottrell(a)wfu.edu&gt; escreveu: On Fri, 16 Apr 2010, Henrique Andrade wrote: Dear Allin, thanks a lot for the script, but unfortunately I think it doesn't do what I need. My aim goal is to build a matrix with this elements: R6_1, R6_2, R6_3, R6_4, R6_5, R6_6, ..., R6_12. I'm not interested in the elements R6_1, ..., R6_5, so it can be set to 0 (or any other value). The number 6 (in R6_x) indicates that we have six elements to sum. To calculate these elements I'm using this: R6_6 = r[6,1] + r[5,2] + r[4,3] + r[3,4] + r[2,5] + r[1,6] R6_7 = r[7,1] + r[6,2] + r[5,3] + r[4,4] + r[3,5] + r[2,6] R6_8 = r[8,1] + r[7,2] + r[6,3] + r[5,4] + r[4,5] + r[3,6] R6_9 = r[9,1] + r[8,2] + r[7,3] + r[6,4] + r[5,5] + r[4,6] R6_10 = r[10,1] + r[9,2] + r[8,3] + r[7,4] + r[6,5] + r[5,6] R6_11 = r[11,1] + r[10,2] + r[9,3] + r[8,4] + r[7,5] + r[6,6] R6_12 = r[12,1] + r[11,2] + r[10,3] + r[9,4] + r[8,5] + r[7,6] Where r[i,j] represents the row "i" column "j" of the "r" matrix. I'm trying to simplify these calculations using the following instructions (I'm just starting to learn how to use loops in my scripts, therefore, I would like to apologize in advance for the simplicity and the errors of my scripts). <script> matrix R = zeros(12,1) loop i=6..12 R_$i = r[$i,1] + r[$i-1,2]+ r[$i-2,3]+ r[$i-3,4]+ r[$i-4,5]+ r[$i-5,6] R[i] = R_$i end loop </script> This script does partly what I want. I need to make it more general in order to allow me to calculate these kind of expressions: R12_6 = r[12,1] + r[11,2] + ... + r[2,11] + r[1,12] R12_7 = r[13,1] + r[12,2] + ... + r[3,11] + r[2,12] I'd attached a .pdf file into this e-mail to show you a little bit better what exactly I need. Thanks for this fix! Best regards, Henrique

On Sat, 17 Apr 2010, Henrique Andrade wrote: On that basis, here's an example: <script> matrix r = mshape(seq(1,144),12,12) print r scalar m = 7 # number of values to compute scalar n = 6 # number of elements to sum matrix R_n = zeros(m,1) loop i=1..m -q loop j=1..n -q row = i + n - j R_n[i] += r[row, j] # see what's going on printf "i=%d: adding r[%d,%d] (%g)\n", i, row, j, r[row, j] endloop endloop print R_n </script> Allin