function Bogenlaenge
clear all

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

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

function y = integr(phi)
adR = 2 ;
y = sqrt(1+adR^2-2*adR*cos(phi)) ;