function fx(x, y, t : real) : real; 
begin
    fx := y
end;
function fy(x, y, t : real) : real; 
begin
    fy := (1.0 - x * x) * y - x
end;