function fp=fpend(y)    % Diffekvationen till pendeln
    global L c  
    fp=[y(2)  -9.81/L*sin(y(1))-c*y(2)*abs(y(2))];
%........................................................



