function x=solve(p,y)
% solves the polynom p for the value at y
%
%   INPUT VALUES:
% input: a polynom
%  
%   RETURN VALUE:
% return value(s) are the values, where the polynom is equal the
% y-value
%
 
% (c) 2003-2008, University of Cambridge, Medical Research Council 
% Maintained by Tom Walters (tcw24@cam.ac.uk), written by Stefan Bleeck (stefan@bleeck.de)
% http://www.pdn.cam.ac.uk/cnbh/aim2006
% $Date: 2008-06-10 18:00:16 +0100 (Tue, 10 Jun 2008) $
% $Revision: 585 $


p(end)=p(end)-y;
x=roots(p);
