How to find the product of each element of a matrix by itself and by the other elemets at the same row?

input:

x = [a,b,c;
     d,e,f]

output:

y = [a^2,a*b,a*c,b^2,b*c,c^2;
     d^2,d*e,d*f,e^2,e*f,f^2]

Answer 1

[Andrei Bobrov]

EDIT

[n1,n2] = ndgrid(1:size(x));
y = x(:,nonzeros(tril(n1))).*x(:,nonzeros(tril(n2)));

[Matt Talebi]

Many thanks Andrei! You've made my day

[Matt Talebi]

Hi Andrei! Sorry it seems that I should have mentioned the dimension of my actual matrix, which is 85 * 15. I'm given this error: "Error using .* Matrix dimensions must agree." I appreciate your help again!

[Andrei Bobrov]

Hi Matt! Fixed.