Question: Fixed point recursion (25 points) Consider the Kepler's equation in celestial mechanics (see Lec. 7, p. 3) of the form r=m+esin c, where we fix

Fixed point recursion (25 points) Consider the Kepler's equation in celestial mechanics (see Lec. 7, p. 3) of the form r=m+esin c, where we fix m = = 0.5. Then, 9 : [0,1] [0,1]. 9(2) Using format long, write a MATLAB code with filename YourlastnameyourfirstnameHW4p1.m to compute the fixed point of the Kepler's equation via fixed point recursion with 20 iterations. Your code should generate a figure plotting the recursion index k (in the horizontal axis) versus Ik in the vertical axis) for 50 randomly chosen initial guesses in [0,1]. For generating random initial guesses in (0,1), you can use the MATLAB command rand. Fixed point recursion (25 points) Consider the Kepler's equation in celestial mechanics (see Lec. 7, p. 3) of the form r=m+esin c, where we fix m = = 0.5. Then, 9 : [0,1] [0,1]. 9(2) Using format long, write a MATLAB code with filename YourlastnameyourfirstnameHW4p1.m to compute the fixed point of the Kepler's equation via fixed point recursion with 20 iterations. Your code should generate a figure plotting the recursion index k (in the horizontal axis) versus Ik in the vertical axis) for 50 randomly chosen initial guesses in [0,1]. For generating random initial guesses in (0,1), you can use the MATLAB command rand
Step by Step Solution
There are 3 Steps involved in it
Get step-by-step solutions from verified subject matter experts
