Finite difference method for second order ode

2021 5 16
1 个回答
回答已采纳
更新时间：2021 5 16
2 次查看（30 天）

0 个投票

Hi. I am trying to solve nonlinear equation using fsolve and repeat graph bet four times with different value. Here is my code:
Function main n = 10 r0 = ones(n, 1) For bet = 0.1:0.1:0.4 Sol = solve(@(r)fun(r, n), r0) norm(fun(sol, n)) end r = ((1:n)-1)/(n-1) Plot(r, sol(0.4,:),'g-*',sol(0.3,:),'m-*',sol(0.2,:),'b-*',sol(0.1,:),'r-*') end

3 个评论
显示 1更早的评论 隐藏 1更早的评论

Jan
Jan 2021-5-16
What is your question? The code does not run. Please format it correctly to improve the readability.
Oluwaseun Adisa
Oluwaseun Adisa 2021-5-16
using fsolve. Here is my code: Function main n = 10 r0 = ones(n, 1) Sol = solve(@(r)fun(r, n), r0) norm(fun(sol, n)) r = ((1:n)-1)/(n-1) Plot(r, sol(1:n)) end
Outcome >> main Solver stopped prematurely Fsolve stopped because it exceeded the function evaluation limit, options.MaxFunEvals = 1000(the default value).
Jan
Jan 2021-5-16
I've mentioned it already: Your code would not run at all, but stop due to the uppercase "Function" immediately.
Please post code in a readable format with one command per line.

请先登录，再进行评论。

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

 采纳的回答

Jan
Jan 2021-5-16
编辑：Jan 2021-5-16

0 个投票

Matlab is case sensitive. There are no "Function", "For" and "Plot" commands. You should ge corresponding error messages.
I try t guess, how your code should be formatted:
function main
n = 10;
r0 = ones(n, 1);
for bet = 0.1:0.1:0.4
sol = solve(@(r) fun(r, n), r0);
norm(fun(sol, n))
end
r = ((1:n)-1) / (n-1);
plot(r, sol(0.4,:), 'g-*', sol(0.3,:), ...
'm-*', sol(0.2,:), 'b-*',sol(0.1,:), 'r-*');
end
end
There is an unmatched end. The loop over "bet" is useless, because the results are overwritten. It is not clear to me, what you want to achieve, so I cannot fix further errors.

更多回答（0 个）

类别

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

产品

标签

Community Treasure Hunt

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

Start Hunting!

Translated by