clear; clc; close all;

%% PASSO 1: ajuste os limites de eixo real e imaginario que vc deseja calcular. 
% Obs.:
%       1) certifique-se que raizes e polos estejam contidos
%       2) tente fazer com que o plano real=0 (ou seja, comecar ou terminar
%       em zero) seja visível para vermos a transformada de Fourier

a = -6:0.05:0;%EDITAR AQUI (ir ate zero ou comecar de zero)
w = -6:0.05:0;%EDITAR AQUI
qtdeTermosA = length(a);
qtdeTermosW = length(w);
%% PASSO 2: digite a equacao H (resposta em frequencia) do seu sistema, considerando s 
H = []; %vai guardar o resultado
for linhas=1: qtdeTermosA
   for colunas=1:qtdeTermosW
       s = a(linhas)+j*w(colunas);
       H(colunas,linhas)= abs((3*s+2)/(s^2+6*s+8)); %EDITAR AQUI
   end 
end
%% PASSO 3: apenas plota. Não tem necessidade de editar nada
mesh(a, w, H); 
xlabel('real'); ylabel('imag'); zlabel('|H|');