function y=g(u,mu)
% function defining fixed points: y=0
y=u^4-u+mu^2-1;