%Pendulum_fun1.m

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


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

u=-K*y;

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

dy=[dy1; dy2];

