function[err,f]=IAE(g,sig)
%[err,f]=IAE(g,sig)
%-------------------------------------------------------|%
%|(C)2005 Bui Trung Hieu				|%
%|Website: www.khvt.com					|%
%|Email: matlab@khvt.com				|%
%-------------------------------------------------------|%
g=minreal(g*sig);
syms x t real
    f=tf2sym(g);        
    sig=zpk(sig);
    if((length(sig.z)+length(sig.p))~=0)
        sigt=tf2sym(sig);
        sigt=ilaplace(sigt,x,t);
    else
        sigt=0;
    end
    
f=ilaplace(f,x,t);
f=abs(sigt-f);
F=inline(f);
err=quad(F,eps,1e2);


%%=========================================================================
%%=========================================================================