Timeline for Basins of attraction using Newton-Raphson method for nonlinear system
Current License: CC BY-SA 4.0
17 events
| when toggle format | what | by | license | comment | |
|---|---|---|---|---|---|
| Mar 16 at 4:02 | comment | added | xzczd♦ | …Why do you directly correct the mistake in your original code (and you've only corrected part of it! ) without any clarification? This only makes the post confusing. Please make a bit more effort in editing your question properly. The new question about efficiency in your comment is definitely on-topic, and your question can be rescued if you properly update your question with this info. (I personally recommend you asking a new question though. ) | |
| Mar 15 at 15:23 | history | became hot network question | |||
| Mar 15 at 13:39 | answer | added | Alex Trounev | timeline score: 8 | |
| Mar 15 at 11:40 | answer | added | Ulrich Neumann | timeline score: 8 | |
| Mar 15 at 10:12 | history | edited | onk | CC BY-SA 4.0 |
deleted 2 characters in body
|
| Mar 15 at 9:58 | answer | added | Daniel Huber | timeline score: 5 | |
| Mar 15 at 9:54 | comment | added | xzczd♦ |
You've made many bad desgin in your code. For example the usage of Inverse. Think carefully about the following: do you really need to use Inverse again and again?
|
|
| Mar 15 at 9:48 | comment | added | onk | The code is working slowly and the number of iterations more than 12 codes it does not work. | |
| Mar 15 at 9:36 | comment | added | onk | Thanks for your help, it's working now. I changed it : Jf = {{3*X[[1]]^2, -1}, {-1, 3*X[[2]]^2}}; | |
| Mar 15 at 9:18 | comment | added | xzczd♦ |
…This won't work, of course. X in your code is a list of numbers, it cannot be used in D. With the code Clear[x, y]; D[{x^3 - y, y^3 - x}, {{x, y}}] I'm just showing you the correct expression for the Jacobian. Once again, compare your Jf with the output of my code carefully and correct your mistake.
|
|
| Mar 15 at 9:04 | comment | added | onk | I changed as below: Jf = D[{X[[1]]^3 - X[[2]], X[[2]]^3 - X[[1]]}, {{X[[1]], X[[2]]}}]; X = X - Inverse[Jf].f, {NoIter}]; | |
| Mar 15 at 8:42 | comment | added | xzczd♦ | That should not happen. How did you change your code? | |
| Mar 15 at 8:38 | comment | added | onk | I have made these changes, but it doesn't work. | |
| Mar 15 at 8:11 | review | Close votes | |||
| Mar 23 at 3:08 | |||||
| Mar 15 at 7:51 | comment | added | xzczd♦ |
And you don't calculate the Jacobian correctly. Just try Clear[x, y]; D[{x^3 - y , y^3 - x}, {{x, y}}] and compare the result to yours.
|
|
| Mar 15 at 7:44 | comment | added | Ulrich Neumann |
Change 1 e -25to 10^-25!
|
|
| Mar 15 at 7:21 | history | asked | onk | CC BY-SA 4.0 |