単純支持・集中荷重のはりの曲げ：荷重点と最大たわみ点は異なる

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 = linspace(a,l,20)';
y1=-b/l*x1.*(a*(a+2*b)-x1.^2);
y2=-a/l*(l-x2).*(b*(2*a+b)-(l-x2).^2);
plot(x1,y1,'b-');
plot(x2,y2,'b-');
xline(x);
xline(a);
big;