[X, y, theta] = ftheta( adj, indsE);
