1
$\begingroup$

I am using matlab to solve the following ODE

$$0.02K^2 F''(K)-3.45 F'(K)-3.45=0.08 F(K) $$

subject to boundary conditions

$$F(K=0)=63.09$$ and $$ F(K=700)=0.$$

I found that I had to give bvpinit an offset such as bvpinit(linspace(1e-3,700,10000),[1 0.1]) from the left boundary 0 because otherwise, I would get an error

Error using ==> bvp4c Unable to solve the collocation equations -- a singular Jacobian encountered

But this isn't really elegant, because the offset seems arbitary. This also did not extend the solution all the way to the left end point 0. Would appreciate if someone can provide a better solution to this problem. Many thanks.

  • 0
    QIang Li, I have a problem similar to yours. Could you solve it? Or with other software different from Matlab? Thank you, Damian.2011-08-11
  • 0
    @Damian: I have converted your answer to a comment. Because you do not have 50 reputation points yet, [you can only comment on your own questions and answers](http://meta.stackexchange.com/questions/19756/how-do-comments-work/19757#19757). So, you didn't do anything wrong; the "add comment" button will only appear for you once you gain 50 points. Here is an [explanation of reputation points](http://meta.stackexchange.com/questions/7237/how-does-reputation-work/7238#7238).2011-08-11

1 Answers 1