Using bvp4c to solve a fourth order nonlinear differential equation with 4 boundary conditions and plot a graph between V vs w

5 次查看（过去 30 天）

PRATEEK KUMAR SINGH 2021-4-15

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/803281-using-bvp4c-to-solve-a-fourth-order-nonlinear-differential-equation-with-4-boundary-conditions-and-p

pks.pdf

here i have used bvp4c solver to make plot between x vs w taking V=10. But i have to plot between V vs w taking x=1.please see attachment for more clarity of question.

function mat4bvp
options = bvpset('RelTol', 1e-6, 'AbsTol', 1e-6, 'NMax', 10000);
solinit = bvpinit(linspace(0,1,10),[0 0 0 0]);
sol = bvp4c(@mat4ode,@mat4bc,solinit,options);
x = linspace(0,1);
y = deval(sol,x);
plot(x,y(1,:));
legend('BVP4C')
xlabel('x')
ylabel('w')
function dydx = mat4ode(~,y)
b=80e-6;L=200e-6;th=4.5e-6;d=3e-6;
I=b*th^3/12;
E=166e9/(1-0.06^2);V=10;
ep=8.8541878e-12;
alfa1 = (b*ep*V.^2*L^4)/(2*d.^3*E*I);
alfa2 = (2*d)/(pi*b);
dydx=[y(2)
      y(3)
      y(4)
     (alfa1+alfa1*alfa2*(1-y(1)))/((1-y(1)).^2)]; %4th order equation
function res = mat4bc(ya,yb)
res=[ya(1)
     ya(2)
     yb(3)
     yb(4)];