以下是mlab22.m的代码，用于在球面上绘制新旧航线：

2011040602.m

Beijing=[40 116];
Shanghai=[31 122];
Tokyo=[36 140];
SanFrancisco=[37 -123];
NewYork=[41 -76];
city=[Beijing;Shanghai;Tokyo;SanFrancisco;NewYork];
R=6400+10;
theta=city(:,1)*pi/180; %第一列，即纬度
fai=city(:,2)*pi/180; %第二列，即经度
x=R*cos(theta).*cos(fai);
y=R*cos(theta).*sin(fai);
z=R*sin(theta);
op=[x,y,z];
format bank
Dmatrix=R*acos(op*op'/R.^2)