function Bogenlaenge
clear all

a = 0 ; b = pi*4  ;              % Integrationsgrenzen

disp(['Bogenlaenge: s = ', num2str(quad (@integr , a , b , 1.e-10))]) ;

function v = integr(phi)
Rdr = 1.5 ; adr = 2 ;
v = (Rdr+1)*sqrt(1+adr^2+2*adr*cos(Rdr*phi)) ;