A newbie question: I have: void foo(matrix<double>& Y, const matrix<double>& X) { ... Y[0] = X[0]; // the compiler chokes on this ... } What is wrong with this? Please tell me. Thank you. Roland