% Função tangente hiperbólica

x=(-4:0.01:4)';
T=(exp(x)-exp(-x))./(exp(x)+exp(-x));

f=1/16;
wo=2*pi*f;

A=[sin(wo.*x) sin(3*wo.*x)/6 sin(5*wo.*x)/120 sin(7*wo.*x)/5040 sin(9*wo.*x)/362880];

a=inv(A'*A)*A'*T; % Parametros da Decomposição de Fourier

Tf=A*a;

E=mean((T-Tf).^2)
figure(1);
plot(x,T,'b',x,Tf,'r');