Shivanand

@Shaun I finally figured out why the solution was not converging. It turned out to be one of my hunches - the precision is not sufficient. Though the paper mentions that they use double precision, perhaps the method that Eigen library uses to solve a system of linear equations does not result in double precision when dealing with matrices of the double data-type. I ended up using an arbitrary precision data type and the solutions converge in a few steps!

@Shaun Thanks. What is MO?

@JasperLoy Both :)

You'll figure out when you actually do research. Shadow someone who does research :D

Howdy, I am implementing a boundary value problem solver in C++ called the initial value adjusting method. Although I have followed the algorithm correctly, the solution does not converge. I'd appreciate if someone can help me. I can share the files with you. The paper is: sciencedirect.com/science/article/pii/0022247X79901495. Asking a question on the forum may not be a good idea as it involves considerable exchange of ideas

mathematica.stackexchange.com/questions/120723/… is the question which needs attention.

@JasonB The second question uses a theta of 85 degrees and that has a non-zero volume of intersection. I will delete the first one.

@JasonB @JasonB What do you suggest?

Nobody is answering my question: mathematica.stackexchange.com/questions/120723/… :(

;3

:)

Hi People. This is my first time on chat