function [a,b,sig2,iflag] = mywe(x,n,ip,iq,epsilon) 
      
lag=ip+iq+1; 
r = correlat(n,lag,0,x,x); 
[a iflag] = toepeqn(r,ip,iq,epsilon);
if iflag == -1
	return 
end
y = czeroflt(x,n,ip,a);
l=round(n/5); 
[b sig2] = durbin(y,n,iq,l);