svfun_uppgc

function dy = svfun(t,y)

global FJADER FREKVENS DAMPNING

d=DAMPNING; 

k=fjader(y(1));

k=fjader(t);

%berakningsrutin for ODE-funktionen

dy=zeros(2,1);

dy(1)=y(2);

dy(2)=driv(t)-k*y(1)-d*y(2);

function s=driv(t)

global FJADER FREKVENS DAMPNING

a=FREKVENS;

s=sin(a*t);

function k=fjader(t)

%skriv in lamplig berakningsrutin nedan.

k=1+cos(t);

%snurra

function k=fjader(y)

n=0;

while n<length(y)

    n=n+1;

    if y(1,n)>=0

        k=10;

    else

        k=0.01;

    end

end