This line of code: asynch = C1(:)*N(:)*C(:)'/(delta-2) throws an error: Error using * Inner matrix dimensions must agree. That would seem to be right C1 is a 1X100 double and N is a 100x100 double. So I changed the code to T = transpose(C1); asynch = C1(:)*N(:)*T(:)'/(delta-2). This makes T a 100X1 double and the code should run but I get the same error. What's wrong. Thanks
The paper from which I got the code uses * and that works in this line: synch = C1(:)*C1(:)'/(delta-2). Not only do I not get an error but I get what I think is the right answer. For asynch = C1(:)*N(:)*C(:)'/(delta-2), I think the problem is that N is 2D and C1 is one D but I got that code from the paper too.