lambda_min=0.4 ;
lambda_max=5 ;

lambdalist=linspace(lambda_min,lambda_max,1000) ;
Glist=zeros(1,length(lambdalist));

Ialist=[1.0000    0.5000    0.1250];
falist=[0     1     5];

for j=1:length(lambdalist)
    Glist(j)=eulot(lambdalist(j),Ialist,falist) ;
end

figure
set(gca,'FontSize',20)
plot(lambdalist,Glist);
xlabel('lambda','FontSize',20)
ylabel('G(lambda)','FontSize',20)


lambdabar = fzero(@(lambda) eulot(lambda,Ialist,falist),1)