function [x]=coren;
a=[1 3 3 1];
[x]=roots(a);
end;