clear all
clc

t=0:pi/100:30;
y=sin(t)./abs(sin(t));

figure();
hold on;
plot(t,y);
syms s j Y;
f=square(j);
F=laplace(f,j,s)
Y1=s*Y;
Y2=s*Y1 - 1;
sol=solve(Y2 + Y1 - F, Y)
sol=ilaplace(sol,j,t)

y2= (exp(-t*pi)-1)./(t.*(exp(-t*pi)+1));
%plot(t, y2);
axis([-1 10*pi+1 -2 2]);