根据实验，得到了密度p与x、y的数据，用二元三次样条插值得到了得到了密度p与x、y的关系，程序如下：

function spline
x=-159:53:159;
y=[-132:44:132]';
p=[0.797243 0.84237 0.846359 0.868923 0.855447 0.849411 0.790338;
   0.865878 0.909008 0.913392 0.912192 0.914284 0.914031 0.871747;
   0.887829 0.911364 0.920252 0.919805 0.920273 0.920954 0.894406;
   0.895495 0.915 0.915352 0.921333 0.914686 0.917614 0.883075;
   0.877245 0.916904 0.920177 0.918321 0.920857 0.911496 0.880024;
   0.845419 0.920637 0.921176 0.915867 0.913707 0.902134 0.833174;
   0.800706 0.862084 0.872675 0.87048 0.867208 0.850913 0.721162];
xi=-159:1:159;yi=[-132:1:132]';
zspline=interp2(x,y,z,xi,yi,'spline');
mesh(xi,yi,zspline);
xlabel('x');
ylabel('y');
zlabel('p');

    根据实验得到了温度t与密度p的数据，拟合得到了温度t与密度p的关系，程序如下：

function wendu
t=[173 169 165 161 157];
m=[0.908 0.90114926 0.882457402 0.855259392 0.795114016];
pp=polyfit(t,m,2);
ti=157:0.1:173;
mi=polyval(pper,ti);
plot(t,m,'o',ti,mi,'k');

    现在想根据：温度t与密度p的关系，密度p与x、y的关系，得到温度t与x、y的关系，这样的程序如何实现？在同一个.m文件中，如何设计主程序和子程序？请高手指点！十分感谢！