城市轮廓线： 输入文件：input1.dat 第一行为建筑个数 之后每行格式为三个数x1 h x2，其中x1,x2是建筑的开始和结束坐标 h是建筑的高度 输出最后轮廓-Urban contour line: the input file: input1.dat first acts after the building number format for each line number 3 x1 h x2, where x1, x2 is the beginn

function dg=gplus(x,y,z,x1,x2,y1,y2,z1,z2,p) G=6.67e-11 t222=myfun(x,y,z,x2,y2,z2) t122=myfun(x,y,z,x1,y2,z2) t212=myfun(x,y,z,x2,y1,z2) t221=myfun(x,y,z,x2,y2,z1) t211=myfun(x,y,z,x2,y1,z1) t121=myfun(x,y,z,x1,y2,z1) t112=myfu