function [B, A] = read_pezfile(filename)
% Return the filter coefficients stored in the given .mat file output by PeZ

structure = importdata(filename)
poles = getfield(structure, 'p_values')
zeros = getfield(structure, 'z_values')

% B = poly(zeros);
% A = poly(poles);

B = zeros;
A = poles;
