Error with Cheb function

2022 4 9
1 个回答
更新时间：2023 2 15
8 次查看（30 天）

0 个投票

Where the ODE is
with boundary condition
y(1)=sin(1)
D = [1/pi*3 ,1];
N = chebop(D);
N.lbc = [1];
N.rbc = 1 - sin(1);
N.op = @(t,x,y) diff(y,2) + (2/x)*diff(y,1) + (1/(x^4))*y;
[x, y] = N\0;
plot([x, y]);

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

Walter Roberson
Walter Roberson 2022-4-9
What difficulty are you experiencing?
Chala Kedir
Chala Kedir 2022-4-9
@Walter Roberson I couldn't plot the output, please if you could help me, it would help me a lot
Walter Roberson
Walter Roberson 2022-4-10
What error message are you getting, at what line?
Have you installed the Chebfun package? https://www.mathworks.com/matlabcentral/fileexchange/47023-chebfun-current-version
Chala Kedir
Chala Kedir 2022-4-10
编辑：Chala Kedir 2022-4-10
@Walter Roberson
this is the error message I get
Error using chebop/linearize
Failed to evaluate operator on the initial guess passed (or the one constructed
by CHEBOP). A potential cause might be division by a zero CHEBFUN. Please supply
a valid initial guess via the 'init' field of the CHEBOP.
Error in chebop/solvebvp (line 109)
[L, residual, isLinear] = linearize(N, u0, x);
Error in \ (line 56)
[varargout{1:nargout}] = solvebvp(N, rhs, varargin{:});
Error in untitled (line 6)
[x, y] = N\0;
Torsten
Torsten 2022-4-10
Are you sure about your D ?
Shouldn't it be
D = [1, 1/3*pi];
?

请先登录，再进行评论。

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

回答（1 个）

Bhavana Ravirala
Bhavana Ravirala 2023-2-15

0 个投票

Hi,
The error you received indicates that the operator fails to evaluate the initial guess. To resolve this error, we can assign the initial guess with the boundary conditions.
N.init = [0; sin(1)];
Hope this helps!!

类别

在 帮助中心 和 File Exchange 中查找有关 Programming 的更多信息

产品

版本

R2022a

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by