отрим N x 1 вектор x и N x N матрицу C. Я хотел бы оценить

 s = x'*C*x;

в Matlab для многих выборок вектора х, например, рассмотрим M выборок x как N x M матрицы X; это можно сделать с помощью

S = diag(X'*C*X);

но это плохое решение, так как в процессе выделяется матрица M x M, и это нарушается при M> 1e5. Есть ли какая-нибудь функциональность matlab, которая может предложить альтернативу?