C-2. (a) 编写程序：震源在地表，选择震源一定范围的入射角，追踪穿过一个地球模型的射线。图形的输出显示震源、射线路径、地球表面、核幔边界以及内核-外核边界。

(b) 用PREM模型追踪震源地表的地震射线，表示出上地幔间断面和地核对射线的影响。

(c) 该程序生成一个走时曲线图，能根据这个图给出上地幔间断面吗？

读入prem模型的速度和深度数据，计算出一组数据，运用matlab绘图。

for jiao=0:0.01:0.1

T=deta+T;

for j=1:phi

结果分析，此图在入射角为0.04的情况下得出。因为速度是连续变化的，所以绘制出图线为曲线。因为上地幔间断面为速度间断面，所以使射线弯折。入射角越小，地震波穿越地球越深。当入射角为0度的时候，地震波射线直接穿过地核，到达对称位置，入射角越大，地震波穿越深度越浅，深度很小时就到达最低点后返回。上地幔间断面和核幔边界因为速度突然减小，所以图线会下凹。核幔边界则会偏向地核。在内外核边界速度增大，本身图线应该下凹，但是下凹程度减小，近似为直线。走时曲线也为曲线，在弯折处为速度突变面。所以第一个弯折出为上地幔间断面。