clc; clear all; 
t=[0:0.01:1];    
m=sinc(2*pi*t); 
subplot(2,1,1) 
hold on    
 plot(m,'black') 
 title('sinc pulse') 
xlabel('time') 
ylabel('amplitude') 
delta=1*pi/100 ;
for n=1:1:100     
if n==1        
e(n)=m(n) ;        
eq(n)=delta*sign(e(n))  ;       
mq(n)=eq(n)     
else  
ix=n
e(n)=m(n)-mq(n-1);        
eq(n)=delta*sign(e(n))         
mq(n)=mq(n-1)+eq(n) ;    
end 
end
stairs(mq,'red') ;
hleg=legend('original signal','stair case approximated signal') ;
hold off   
subplot(2,1,2) 
hold on 
m1=sin(2*pi*t);
 plot(m1,'red') 
title('sin wave')
 xlabel('time') 
ylabel('amplitude') 
 for n=1:1:100     
if n==1         
e1(n)=m1(n);         
eq1(n)=delta*sign(e1(n));         
mq1(n)=eq1(n);     
else         
e1(n)=m1(n)-mq1(n-1);         
eq1(n)=delta*sign(e1(n));         
mq1(n)=mq1(n-1)+eq1(n);     
end 
end 
stairs(mq1,'blue') 
hleg=legend('original signal','stair case approximated signal') 
hold off
