Unexpected error
function dFdx=funode(x,F)
global Lamda
global K1
global P
dFdx=[F(2) ; -K1*P+F(1)*Lamda^2];
function res=myfunbc(Fa,Fb)
global P
res=[Fa(1) ; Fb(1)-P];
solinit=bvpinit(linspace(0,l,2),[1 0]);
sol = bvp4c(@funode,@myfunbc,solinit);
x = linspace(0,l);
F = deval(sol,x);
After that i have an error The derivative function ODEFUN should return a column vector of length 2. why?

 채택된 답변

James Tursa
James Tursa 2017년 9월 29일
편집: James Tursa 2017년 9월 29일

0 개 추천

Make it a column vector as requested:
dFdx=[F(2) ; -K1*P+F(1)*Lamda^2];

댓글 수: 11

Nadim Mhanna
Nadim Mhanna 2017년 9월 29일
still the same error
What about the other function:
res=[Fa(1) ; Fb(1)-P];
Nadim Mhanna
Nadim Mhanna 2017년 9월 29일
편집: Nadim Mhanna 2017년 9월 29일
i did corrected both but didn't work
Steven Lord
Steven Lord 2017년 9월 29일
What is the size of each of the variables K1, P, and Lamda?
Nadim Mhanna
Nadim Mhanna 2017년 9월 29일
편집: Nadim Mhanna 2017년 9월 29일
They are constants i.e. k1=5 Lamda=13 P=252
Where is the value of P coming from in your lines
function res=myfunbc(Fa,Fb)
res=[Fa(1) Fb(1)-P];
You are not using the form of "function" that would allow P to be a shared variable.
Nadim Mhanna
Nadim Mhanna 2017년 9월 30일
I did it by entering global P
Nadim Mhanna
Nadim Mhanna 2017년 9월 30일
For all the variables I put the global value that i did assign but still not working
Walter Roberson
Walter Roberson 2017년 9월 30일
편집: Walter Roberson 2017년 9월 30일
When you use a global variable, it is not uncommon for the global variable to turn out to be empty at the time it is used.
Nadim Mhanna
Nadim Mhanna 2017년 9월 30일
okay got it Thanks all
Jan
Jan 2017년 9월 30일
@Nadim Mhanna: Please do not use flags to mark a "best answer". Flagging is thought to inform admins and editors, that a contribution needs attention, because it violates the terms of use, e.g. by rudeness.

댓글을 달려면 로그인하십시오.

추가 답변 (0개)

카테고리

도움말 센터File Exchange에서 Programming에 대해 자세히 알아보기

질문:

2017년 9월 29일

댓글:

Jan
2017년 9월 30일

Community Treasure Hunt

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

Start Hunting!

Translated by