%Pendulum_fun2.m

function dy=Pendulum_fun2(t,y,a,b,c,delta,Tss,K)


y1=y(1);
y2=y(2);
sigma=y(3);

u=-K*y;

dy1=y2;
dy2=-a*sin(y1+delta)-b*y2+c*(u+Tss);
dsigma=y1;

dy=[dy1; dy2; dsigma];

