function dy=phase_pendulum_fun(t,y)

%[m/s^2]
g=10;
%[m]
l=10;
%[kg*m^2/s]
b=5;
%[kg]
m=0.1;

dy1=y(2);
dy2=-(g/l)*sin(y(1))-(b/m/l^2)*+y(2);

dy=[dy1; dy2];