function Nullstelle2

lam = 1e-6 ;
f0  = ewg (lam)  ;

while lam < 10
    lam = lam + 0.001 ;
    f1  = ewg(lam) ;  
    if f0*f1 <= 0  
        [lam1 flam1 exitflag] = fzero (@ewg , lam) ;
        if exitflag == 1
            Nullstelle_lam1 = lam1
            break ;
        end
    end
    f0 = f1 ;
end

function f = ewg (lam)
f = tan(lam)-tanh(lam) ;