No graphs, Mechanical Vibrations desperate for help

b) Let the force be accelerating such that (F sin(5t^2 /1000)). Plot the response (for c=0.2) for each mass as function of angular frequency. Explain the result

I have already solved part a) and i am currently stuck at part B)

It requires me to plot the displacemnet as a function of driving frequencies and i should observe the natural frequencies at the right frequencies.

BUT after trying, plot (w,x), all i get is an empty graph

Any idea how to go about it? I'm "close" to the right solution.

% Simulation of forced 2DOF system with ode45
clear all
clc
global A B F Z I wdr
m=1;
c=0.2;
k=100;
F0=10;
c=0.2;
M=[m 0 0 ;0 m 0 ; 0 0 m];
K=[2*k -k 0;-k 2*k -k;0 -k k];
C=[2*c -c 0;-c 2*c -c;0 -c c];
f=[0 0 F0]'; % Force vector
 t0=0;
 t1=100;
    for t=t0:t1
        wdr=(t/100)%driving frequency
    end
A=M\K;
B=M\C;
F=M\f;
Z=zeros(3); % Zero matrix
I=eye(3); % Diagonal matrix
x0=[0 0 0 0 0 0]'; % Initial conditions
t0=0;
t1=100;
[t,x]=ode45('Fun1',[t0 t1],x0);
w= 0:1:3632;
plot(wdr,x(:,1),wdr,x(:,2))

Fun1 file:

function dx=Fun1(t,x)
global A B F Z I wdr
dx=zeros(6,1);

dx=[Z I ; -A -B] *x+[[0 0 0]' ; F]*sin(5*t^2/1000);

end