calculate shear force diagram and bending moment diagram and deflection by using integration method in matlab

main.m
% Define beam properties
L = 5; % Beam length (m)
P = 10; % Point load at a distance 'a' from left support (kN)
a = 3; % Distance of point load from left support (m)

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

% Distance along the beam
x = 0:0.1:L;

% Shear force equation
V = @(x) Ra * (x > a) - P;

% Bending moment equation
M = @(x) Ra * x .* (x <= a) - P * (x - a) .* (x > a);

% Deflection equation
E = @(x) Ra * x.^2 / (2 * 24 * 10^3 * 8) .* (x <= a) - P * (x - a).^2 / (2 * 24 * 10^3 * 8) .* (x > a);

% Calculate shear force, bending moment, and deflection
shear_force = V(x);
bending_moment = M(x);
deflection = E(x);

% Plotting the diagrams
subplot(2, 2, 1)
plot(x, shear_force, 'r')
title('Shear Force Diagram')
xlabel('Distance (m)')
ylabel('Shear Force (kN)')
grid on

subplot(2, 2, 2)
plot(x, bending_moment, 'b')
title('Bending Moment Diagram')
xlabel('Distance (m)')
ylabel('Bending Moment (kNm)')
grid on

subplot(2, 2, 3)
plot(x, deflection, 'g')
title('Deflection Curve')
xlabel('Distance (m)')
ylabel('Deflection (m)')
grid on
1075 chars
48 lines