calculate shear force diagram and bending moment diagram and deflection for indetermined beam in matlab

main.m
% Define the beam properties
L = 5; % Length of the beam (in meters)
E = 200e9; % Young's modulus of the beam material (in Pascal)
I = 1.3333e-5; % Moment of inertia of the beam (in meters^4)
q = 1000; % Uniformly distributed load on the beam (in Newton/meter)

% Create symbolic variables
syms x

% Define the load and boundary conditions
V = @(x) -q*x; % Shear force equation
M = @(x) -q*x.^2/2; % Bending moment equation
w = @(x) -q*x.^3/(6*E*I) - (q*L^2*x)/(2*E*I) + (q*L^3)/(6*E*I); % Deflection equation

% Plotting the shear force diagram
figure
fplot(V, [0 L], 'LineWidth', 1.5)
title('Shear Force Diagram')
xlabel('Position along the beam (m)')
ylabel('Shear Force (N)')
grid on

% Plotting the bending moment diagram
figure
fplot(M, [0 L], 'LineWidth', 1.5)
title('Bending Moment Diagram')
xlabel('Position along the beam (m)')
ylabel('Bending Moment (Nm)')
grid on

% Plotting the deflection diagram
figure
fplot(w, [0 L], 'LineWidth', 1.5)
title('Deflection Diagram')
xlabel('Position along the beam (m)')
ylabel('Deflection (m)')
grid on
1051 chars
38 lines