% Краевая задача, пример 3: прогиб балки
l= 10; EI= 1;
p= 0.0001;
f= @(x)0.0001*sin(x*pi/l);
rp= @(x,y) [y(2);y(3);y(4);(-p*y(3)-f(x))/EI];
% Концы балки на шарнирах
bc= @(ya,yb) [ya(1);ya(3);yb(1);yb(3)];
xinit= 0:l/60:l;
solinit = bvpinit(xinit,[0 0 0 0]);
sol= bvp4c(rp,bc,solinit)
plot(sol.x,sol.y(1,:));

