clear; close all; l = 5; a =4.5 disp('荷重作用点'); b = l-a; x = sqrt(1/3*(l^2-b^2)) disp('最大たわみ点'); dx =linspace(0,l,20)'; dy =zeros(length(dx),1); plot(dx,dy,'r-'); hold on; plot(x,0,'go'); plot(a,0,'k*'); x1 =linspace(0,a,20)'; x2…