国土面积测量问题

如图所示一个国家的地图边界线，有关方面需确认国家的国土面积的数值。现在将地图放入平面直角坐标系中，以由西向东方向为 x 轴，由南向北方向为 y 轴，在边界上选取若干 点测出个点对应的坐标数据（?,?）和（?,?）(i=1,2,…),地图 的比例：18mm 相当于 40km,试请你利用测量数据计算该国的国土面积的大小，并与已有参考值 41288??进行比较分析。

x=[7 10 13.0 17.5 34 40.5 44.5 48 56 61 68.5 76.5 80.5 91 96 101 104 106.5 111.5 118 123.5 136.5 142 146 150 157 158]'; 
y1=[44 45 47 50 50 38 30 30 34 36 34 41 45 46 43 37 33 28 32 65 55 54 52 50 66 66 68]'; 
y2=[44 59 70 72 93 100 110 110 110 117 118 116 118 118 121 124 121 121 121 122 116 83 81 82 86 85 68];
xi=(0:0.001:158)';%插值步长 
yi1=interp1(x,y1,xi,'spline');%插值 
yi2=interp1(x,y2,xi,'spline');%插值 
axis([0 160 20 140]);%设置坐标轴 
hold on; plot(x,y1,'b',x,y2,'g')%绘图 
integ=trapz(xi,[yi1,yi2]);%用梯形法计算两个积分 
disp('此国国土面积为：'); 
disp((integ(2)-integ(1))/18/18*1600); disp('km^2')