% unit step function
% clc;
% clear all;
% close all;
% x=-10:0.01:10;
% for i=1: length(x)
% if (x(i))>=0;
% y(i)=1;
% else 
% y(i)=0;
% end 
% end
% subplot(221)
% plot(x,y)
% title("unit function")
% xlabel('time');
% ylabel('amplitude');