Re: Eigensystem

*To*: mathgroup at smc.vnet.net*Subject*: [mg39664] Re: Eigensystem*From*: Sseziwa Mukasa <mukasa at jeol.com>*Date*: Sat, 1 Mar 2003 02:47:15 -0500 (EST)*Sender*: owner-wri-mathgroup at wolfram.com

On Friday, February 28, 2003, at 04:46 AM, Angela Birk wrote: > Dear MathGroup Members, > > I would like to generate a 12x12 matrix for that I can caluclate the > Eigensystem by transposing a 12x12 matrix. This 12x12 matrix is the > result of a model that includes a bunch of functions and variables. > (Note the attached file doesn't work, because the model is missing. But > one can see that the mistake occurs.) > [contact the author to get the attachment - moderator] > > I try to derive the 12x12 matrix by using different commands: > > First, I define the first row of the matrix (see attached file): > > f[a_,b_] := f[a,b]=Table[(1-CDF[h[a,bValue,W]])* > SETP[bValue,b[[i]]]*ATP[a, j] + CDF[h[a,bValue,W]]* > SETP[bValue,b[[i]]]*ATP[a, j], {j, 1, 4}, {i, 1,3}] > > where CDF[h[a,bValue,W]] is a probability that depends on the function > h > and SETP[bValue,b[[i]]] as well as ATP[a, j] are given matrices having > the dimensions 3x3 and 4x4, respectively. > > Second, if I execute the command > > f[1, b[[1]]] // MatrixForm > > this should give a 1x12 vector. But actually I get a 4x3 matrix (see > attached file). This 1x12 vector is the first row of the 12x12 matrix, > and the sum of row is one. (Each row has to sum to one.) > > Could someone suggest how the above command(s) should be changed to get > the 1x12 vector that sums to one? > > f[a_,b_]:=f[a,b]=Flatten[Table[(1-CDF[h[a,bValue,W]])* SETP[bValue,b[[i]]]*ATP[a, j] + CDF[h[a,bValue,W]]* SETP[bValue,b[[i]]]*ATP[a, j], {j, 1, 4}, {i, 1,3}]] Regards, Sseziwa