开普勒方程的matlab实现

开普勒方程是描述行星运动的重要方程，其matlab实现如下：

function E = keplerEquation(M, e)
    % 输入参数：
    % M: 平均近点角（Mean Anomaly）
    % e: 偏心率（Eccentricity）
    % 输出参数：
    % E: 真近点角（Eccentric Anomaly）

    % 定义迭代初始值
    E0 = M;

    % 定义迭代精度    epsilon = 1e-8;

 % 开始迭代
    while true
        E = E0 (E0 - e * sin(E0) - M) / (1 e * cos(E0));
        if abs(E - E0) < epsilon
            break;
       
        E0 = E;
    end
end

相关问题：
1. 什么是开普勒方程？
2. 如何理解平均近点角和真近点角？
3 什么是偏心率？
4. 迭代精度的选择什么影响？