n = -10:10;
x=(n>=0);

Energy = sum(abs(x).^2);
stem(n,x);
grid on ;
title('Unit Step Signal');
disp(['Energy = ',num2str(Energy)]);