How can we use a Heaviside function for differential equation having moving load ?

I have a differential equation mentioned in function file where i have to apply moving load. I want to verify the approach that i have tried to get a moving load using heaviside function. so please help me implimenting correctly the moving load.

clc
clear all

%%

% time span
 
tspan=[0:0.1:10];

 y0=[0;0];
%% moving load

F0=9.8;      % magnitude of moving load
L=4;     % length
V=0.1;    % velocity
m=28.152;
C=500;
   [t,y]=ode45(@(t,y) Lin(t,y,C,w,F0,m,L,V),tspan,y0);

function [f]=Linear12(t,y,C,F0,m,L,V)
%% moving load

 p1= F0*sqrt(2/m*L)*sin((pi*V*t)/L).*heaviside(L/V-t);
%%
f= zeros(2,1);

%%  for beam with linear absorber

f(1)=y(2);
f(2)= -50^2*y(1)-C*y(2)+p1;


end