function Nullstelle1

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

while lam < 10
    lam = lam + 0.001 ;
    f1  = ewg(lam) ;  
    if f0*f1 <= 0  
        lam1 = fzero (@ewg , lam)
        break ;
    end
    f0 = f1 ;
end

function f = ewg (lam)
f = cosh(lam)*sin(lam)-cos(lam)*sinh(lam) ;