error is using bvpset in bvp4c

Question

Syed 2024-9-24

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2155000-error-is-using-bvpset-in-bvp4c

评论： Syed 2024-9-25

Hi, I am Syed,

I like to have an output with 6 or more than 6 decimal place (for example: 3.033334), without bvpset i am getting the output upto 4 decimal places, but i like to have upto to 6 or 8 decimal places, please guide me to get the desired output

The error i am getting is

Error using <=

Not enough input arguments.

My program is given below

close all

clc

phi=0.1;

phi1=2;

phi2=2;

Betaf=0.00021;

Betas=0.000058;

sigmaf=0.005;

sigmas=580000;

GR=2;

EM2=20;

Rof=997.1;

Ro=5;

c=0;

Ros=4420;

h=0.5;

Hc=0.25;

K=1;

R=1;

e1=(1./((1-phi).^2.5));

e2=(1-phi)+(phi.*(Ros./Rof));

e3=(1-phi)+phi.*((Ros.*Betas)./(Rof.*Betaf));

sigma=(sigmas./sigmaf);

e4=1-((3*(1-sigma).*phi)./((2+sigma)+(1-sigma).*phi));

gama1=e2./e1;

gama2=1/e1;

gama3=e3./e1;

s1=(e4./(phi1.*(1-i.*Hc)));

s2=(e4./(phi2.*(1-i.*Hc)));

dydx=@(x,y)[y(3);

y(4);

-EM2.*gama2.*y(4)+2.*i.*Ro.*gama1.*y(1)+gama3.*GR.*x+R.*gama1;

(-e4./(1-i.*Hc)).*y(3)];

BC = @(ya,yb)[ya(1);yb(1)-c;ya(4)-s1.*ya(2);yb(4)+s2.*yb(2)];

yinit = [0.01;0.01;0.01;0.01];

options = bvpset('stats','on','RelTol',le-3);

solint = bvpinit(linspace(0,1,50),yinit,options);

U1 = bvp4c(dydx,BC,solint);

dudy1=U1.y(3,1)

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Shashi Kiran 2024-9-24

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2155000-error-is-using-bvpset-in-bvp4c#answer_1521610

在 MATLAB Online 中打开

Hi @Syed,

I understand that you are aiming to obtain an output with 6 or more decimal places but are encountering errors.

Here are my observations and corrections for your code:

The error was because le-3 is not a valid expression. It should be 1e-3, which is the proper MATLAB syntax for .

options = bvpset('stats','on','RelTol',1e-3);

The bvpinit function is used to initialize the solution guess and mesh, and it doesn't accept the options structure as an argument.

solint = bvpinit(linspace(0,1,50),yinit);

The correct place to pass the options structure is in the bvp4c function call. This allows the solver to use the relative tolerance (RelTol) and print statistics ('stats', 'on').

U1 = bvp4c(dydx,BC,solint,options);

Here is the output following these corrections:

Refer to the following documentations for more details about the functions:

Hope this solves your query.

4 个评论
显示 2更早的评论隐藏 2更早的评论

Shashi Kiran 2024-9-25

Additionally,

To set the default formatting for long in MATLAB, follow these steps:

- On the MATLAB homepage, go to Home.

- Navigate to Preferences.

- Under Command Window, select Text Display, then choose Numeric Format and set it to long.

Syed 2024-9-25

I got it now, thanks a lot

请先登录，再进行评论。

error is using bvpset in bvp4c

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

4 个评论
显示 2更早的评论隐藏 2更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

error is using bvpset in bvp4c

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

4 个评论 显示 2更早的评论隐藏 2更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

4 个评论
显示 2更早的评论隐藏 2更早的评论