function yp=Integrador(t,y,A)
%%
% Prepara para integrar o sistema yp=A*y+b
% Note que as dimensões da matriz A devem ser compatíveis com a do
% vetor-coluna y
N=length(y);
f0=1.2;
w=1.5;
b=zeros(N,1); %modificar a linha correspondente, caso queira incluir algum termo não nulo
b(2)=f0*sin(w*t);
yp=A*y+b;