beam fixed by two supports and concrete load calculate shear force diagram and bending moment diagram and deflection in matlab

main.m
% Define parameters
L = 5; % Length in meters
E = 200e9; % Young's Modulus in Pa
I = 1e-5; % Moment of inertia in m^4
w = 10; % Uniform load in kN/m
P = 20; % Point load in kN

% Create distance vector
x = linspace(0, L, 1000);

% Calculate shear force
V = zeros(size(x));
for i = 1:length(x)
    V(i) = w*x(i) - P;
end

% Calculate bending moment
M = zeros(size(x));
for i = 1:length(x)
    M(i) = w*x(i)^2/2 - P*(x(i)-2) - w*x(i)*(x(i)/2);
end

% Calculate deflection
D = zeros(size(x));
for i = 1:length(x)
    D(i) = -w*x(i)^4/(24*E*I) + P*x(i)^3/(6*E*I) + w*x(i)^2/(2*E*I) + C1*x(i) + C2;
end

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

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

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