base.m

我自己编的计算地质埋藏史用的一维回拨程序。

function H=base(fai0,c,h1,y1,y2)
h2=4000;
H=0;
eps=0.001;
while abs(h2-H)>eps   
      h2=H;
      f=h1-h2+y2-y1-fai0/c*((exp(-c*y1))-exp(-c*y2)-exp(-c*h1)+exp(-c*h2)); % 目标函数
      k=-1+fai0*exp(-c*h2);   % k是数值逼近法中的斜率
      H=-f/k+h2;
end