% Biegeschwingungen gerader Traeger, Rayleighscher Quotient, Aufgabe 32-6.
function Aufg32_6_Rayl
clear all

% Parameter:
tl = 1 ;  EI = 3000 ;  rhoA = 3 ;

% Ansatzfunktion und deren 2. Ableitung:
global vS vSDD
vS   = [2/tl^4 -5/tl^3 3/tl^2 0 0] ;  % Biegelinie bei konstanter Linienlast
vSDD = polyder (polyder (vS));        % 2. Ableitung der Ansatzfunktion

k11 = EI  *quad(@k11fun,0,tl) ;
m11 = rhoA*quad(@m11fun,0,tl) ;

Eigenfrequenz = sqrt(k11/m11) / (2*pi)

function y = m11fun(x) 
global vS vSDD
y = (polyval(vS , x)).^2 ;

function y = k11fun(x) 
global vS vSDD
y = (polyval(vSDD , x)).^2 ;