Ialist=[1 .5 .125] ;
falist=[0 1 5];

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