%données de l'exercice
L = 1.1;
A = 2.7;

%description du batonnet:
xmin = -L/2;
xmax = L/2;
%Nombre de points utilisés pour décrire le batonnet
N = 5100;


%Definition du pas d'espace
dx = (xmax-xmin)/N;

%Initialisation de la charge totale
Qtot = 0;

%Incrémentation de la charge le long du batonnet
for x = xmin:dx:xmax
    lambda = A * (1 + 24*(x/L).^2);%lambda est la densité linéique
    Qtot = Qtot + lambda * dx;
end
 
%Affichage de la valeur numérique obtenue
Qtot
%Affichage de la valeur théorique
Qtheo = A * L* 3

%Tracé de la densité
x = linspace(xmin,xmax,N);
lambda = A * (1 + 24*(x/L).^2);%lambda
plot(x,lambda,'+-')