assume(n,integer);

f1:=z*a/d;

f2:=(L-z)*a/(L-d);

B(n):=int(f1*sin(n*Pi*z/L),z=0..d)+int(f2*sin(n*Pi*z/L),z=d..L);

simplify(%);

L:=7*d; a:=1; d:=1;

S2:=(2/L)*sum(B(n)*sin(n*Pi*x), n=1..25):

plot(S2,x=0..1);

S3:=plot(S2,x=0..1,color=magenta):

with(plots):

S4:={(2/L)*(int(f1*sin(i*Pi*z/L),z=0..d)+int(f2*sin(i*Pi*z/L),
z=d..L))*sin(i*Pi*x) $ i=1..5}:
S5:=plot(S4,x=0..1):

display([S5,S3]);