function wave
x10=...;
x20=...;
[T,Y] = ode45(@rhs,[0 12],[x10 x20]);
plot(T,Y(:,1),'-',T,Y(:,2),'-.');
function dy = rhs(t,y)
w0=...;
n=...;
x=...;
L=...;
An=...;
Ant=...;
Ft=...;
dy = zeros(2,1); % a column vector
dy(1) = y(2);
dy(2) = -w0^2*(n*x/L)^2*y(1)+(An+Ant)*Ft;
end
Best wishes
Torsten.
