BVP=Boundary va2bvplue problem
4-4th order,r, the tolerance is also set as 10-6
It needs a format structure.
clear all;
init=bvpinit(linspace(1，3,10)，[0,0]) % [0,0] is the initial guess
sol=bvp4c(@bvp_rhs, @bvp_bc,init);
x=linspace(1,3,100)
BS=deval(sol,x);
plot(x,BS)
function rhs=rhs_bvp(x,y);
rhs=(y(2);-3*y(2)-6*y(1)+5);
% dveval- already know how to
http://www.doc88.com/p-8856140196042.html
https://www.youtube.com/watch?v=i7U-ra5lTN8
https://www.youtube.com/watch?v=iEep1-WnjlM
https://cn.mathworks.com/help/matlab/ref/bvp4c.html