how to calculate deflection of a beam in matlab

main.m
% Define beam parameters
L = 1;              % Length of the beam
E = 2e11;           % Young's modulus of the beam material
I = 8e-6;           % Second moment of area of the beam
w = 1000;           % Distributed load on the beam

% Define discretization of the beam
N = 100;            % Number of elements
x = linspace(0,L,N); % Position vector along the beam

% Calculate the deflection of the beam
d = zeros(1,N);
for i = 1:N
    d(i) = (w*x(i)^2)/(6*E*I)*(3*L-x(i));
end

% Plot the deflection of the beam
plot(x,d);
xlabel('Position along beam');
ylabel('Deflection');
title('Deflection of the beam');
610 chars
22 lines