calculate shear force diagram and bending moment diagram and deflection in matlab

main.m
% Beam properties
L = 5; % Length of the beam (in meters)
P = 10; % Point load at a distance of 3 meters from the left end

% Calculate reactions at the supports
R1 = P*(L-3)/L;
R2 = P*(3)/L;

% Initialize arrays to store values
x = 0:0.01:L; % Positions along the beam
V = zeros(size(x)); % Shear force array
M = zeros(size(x)); % Bending moment array
D = zeros(size(x)); % Deflection array

% Calculate shear force and bending moment
for i = 1:length(x)
    if x(i) <= 3
        V(i) = R1;
        M(i) = R1*x(i);
    else
        V(i) = R1 - P;
        M(i) = R1*x(i) - P*(x(i)-3);
    end
end

% Calculate deflection using the bending moment
for i = 1:length(x)
    D(i) = (1/(6*E*I))*integral(@(x) M(i), 0, x(i));
end

% Plot the diagrams
figure;
subplot(3,1,1);
plot(x,V);
title('Shear Force Diagram');
xlabel('Position (m)');
ylabel('Shear Force (N)');

subplot(3,1,2);
plot(x,M);
title('Bending Moment Diagram');
xlabel('Position (m)');
ylabel('Bending Moment (N.m)');

subplot(3,1,3);
plot(x,D);
title('Deflection Diagram');
xlabel('Position (m)');
ylabel('Deflection (m)');
1084 chars
50 lines