Error with Cheb function

6 次查看（过去 30 天）
Chala Kedir
Chala Kedir 2022-4-9
回答： Bhavana Ravirala 2023-2-15
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更早的评论
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
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!!

类别

在 Help Center 和 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