用matlab写一段程序，两个地理AB坐标Lbh之间受地球曲率的影响，计算在地球曲率影响下从A点看到B点最远能看多少公里

在MATLAB中，我们可以使用地球的半径和两点之间的球面距离来计算从A点看到B点最远能看多少公里。以下是一个简单的MATLAB程序示例：

function max_distance = calculate_max_distance(lat1, lon1, lat2, lon2)
    % 输入参数:
    % lat1, lon1 - A点的纬度和经度
    % lat2, lon2 - B点的纬度和经度
    % 输出参数:
    % max_distance - 从A点看到B点的最大可视距离（公里）

    % 地球半径（单位：公里）
    R = 6371;

    % 将经纬度转换为弧度
    lat1 = deg2rad(lat1);
    lon1 = deg2rad(lon1);
    lat2 = deg2rad(lat2);
    lon2 = deg2rad(lon2);

    % 计算两点间的球面距离
    dlat = lat2 - lat1;
    dlon = lon2 - lon1;
    a = sin(dlat/2)^2 + cos(lat1) * cos(lat2) * sin(dlon/2)^2;
    c = 2 * atan2(sqrt(a), sqrt(1-a));
    distance = R * c;

    % 假设视线与地球表面垂直，则最大可视距离等于球面距离
    max_distance = distance;
end

这个函数`calculate_max_distance`接受两个地理坐标（纬度和经度），并返回从A点到B点的最大可视距离。这里我们假设视线是垂直于地球表面的，因此最大可视距离等于两点间的球面距离。

你可以调用这个函数并传入具体的经纬度值来获取结果。例如：

lat1 = 34.0522; % A点的纬度
lon1 = -118.2437; % A点的经度
lat2 = 40.7128; % B点的纬度
lon2 = -74.0060; % B点的经度

max_dist = calculate_max_distance(lat1, lon1, lat2, lon2);
fprintf('从A点看到B点的最大可视距离为：%.2f公里\n', max_dist);

这段代码会计算并打印出从洛杉矶到纽约的最大可视距离。