clear, clc, clf

(dane…)

Ix=b*h^3/12;

k=48*E*Ix/L^3; %belka wolnopodparta lub: 3*E*Ix/L^3 wspornik

m=b*h*ro*L/2; %belka wolnopodparta lub: m=b*h*ro*L/4 wspornik

wn=sqrt(k/m);

ckr=2*m*wn;

c=ckr*ksi;

%drgania tłumione

wd=wn*sqrt(1-ksi^2);

Td=2*pi/wd;

%drgania swobodne

A=u0

B=v0/wn

un=A*cos(wn*t)+B*sin(wn*t);

unt1=un(find(t==t1)) %ważne podwójne = !!!; szukane przemieszczenie w drganiach swobodnych

%drgania wymuszone

dtkr=2/wn

dt=… %(dobieramy sporo mniejszy)

tk=…

t=[0:dt:tk];

t_1=[0:dt:td];

p1=p0*sin(pi*t_1/td); %plus ewentualnie odpowiednia liczba w nawiasie, warto narysować ten wykres i sprawdzić czy wygląda jak ten dany

t_2=[td+dt:dt:tk];

p2=0*t_2;

p=[p1 p2];

[u,v,a]=mrc(m,c,k,p,t,u0,v0); %ważne!!! u0 i v0 teraz mają być zerowe

ut2=u(find(t==t1)) %szukane przemieszczenie w drganiach wymuszonych