calculate deflection and shear force diagram and bending moment diagram for beam fixed by two support in matlab

main.m
% Define the beam properties
L = 5; % Length of the beam (in meters)
P = 10; % Load applied at a distance 'a' from support 1 (in kN)
a = 2; % Distance of the point load from support 1 (in meters)
E = 200 * 10^9; % Young's Modulus of the material (in Pa)
I = 4 * 10^-6; % Area moment of inertia of the beam (in meters^4)

% Create a vector representing the length of the beam
x = 0:0.01:L;

% Calculate the reactions at the supports
Ra = (P*a) / L;
Rb = P - Ra;

% Calculate the shear force (V) and bending moment (M) along the beam length
V = zeros(size(x));
M = zeros(size(x));

for i = 1:length(x)
    if x(i) < a
        V(i) = Ra;
        M(i) = Ra * x(i);
    else
        V(i) = Ra - P;
        M(i) = Ra * x(i) - P * (x(i) - a);
    end
end

% Plot the Shear Force diagram
figure;
plot(x,V,'r-','LineWidth',1.5);
grid on;
xlabel('Beam Length (m)');
ylabel('Shear Force (kN)');
title('Shear Force Diagram');

% Plot the Bending Moment diagram
figure;
plot(x,M,'b-','LineWidth',1.5);
grid on;
xlabel('Beam Length (m)');
ylabel('Bending Moment (kN.m)');
title('Bending Moment Diagram');
1091 chars
44 lines