% y=[];
% for t=0:0.1:3;
%     if(t>=0&&t<1)
%         y=[y,0];
%     elseif(t>=1&&t<2)
%         y=[y,1];
%     elseif(t>=2&&t<=3)
%         y=[y,0];
%     end
% end
% t1=0:0.1:3;
% it=find(t1>0.5)
% nit=size(it,2);
% it(2:nit)=[];
% u=size(y,2);
% y(it:u)=y(it:u)+0.5*y(1:u-it+1)
% y;
% plot(t1,y)
% x=1:10;
% z=x'*x
% t=abs('a')
% c=abs('i')
% g=abs(' ')
% d=char([t c g c t])
z=-10:0.1:10;
p=[3 4 7 2 1];
y=polyval(p,t);
yf=fft(y);
plot(abs(yf))

s=0;
dt=0.000001;
for t=-3:dt:1-dt;
  s=s+polyval(p,t)*dt; 
end
 s
syms x
p=[3 4 7 2 1];
A = poly2sym(p,x);
S=int(A,-3,1);
double(S)