On 6/25/09, David Ronis <ronis at ronispc.chem.mcgill.ca> wrote:
> kill(all);
> load("eigen");
> k:columnvector([k[1],k[2]]);
> r:columnvector([r[1],r[2]]);
> v:columnvector([v[1],v[2]]);
>
> f:transpose(k).v + transpose(r).k;
>
> ev(f,v[1]=0, v[2]=0);
OK, on looking at this again I see that the problem is that you
have v[1] = [v[1]] because of v : columnvector([v[1], v[2]]).
Try vv : columnvector([v[1], v[2]) and I think then the ev stuff works.
Hope this helps,
Robert Dodier