calculate deflection of the beam by integration method in matlab

main.m
% Given parameters
L = 5; % Length of the beam
P = 10; % Point load applied at the center of the beam
E = 200e9; % Young's modulus of the material
I = 2.5e-6; % Moment of inertia of the beam

% Define the load function
w = @(x) P/(2*L) * (L - abs(x - L/2));

% Calculate the deflection using integration method
syms x;
syms C1 C2;
theta = int(w(x)*x/(E*I), 0, L) + C1;
y = int(theta, 0, x) + C2;

% Applying boundary conditions
eqn1 = subs(y, x, 0) == 0; % Deflection at x=0 is 0
eqn2 = subs(y, x, L) == 0; % Deflection at x=L is 0

[CSolC1, CSolC2] = solve([eqn1, eqn2], [C1, C2]);

deflection = subs(y, [C1, C2], [CSolC1, CSolC2]);
disp(deflection);
652 chars
24 lines