base2.m

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

function H=base(fai0,c,h1,y1,y2)
h2=4000;
H=2000;
eps=0.1;
n=1;
while abs(h2-H)>eps   
      H=h2;
      h2=h1+y2-y1-fai0/c*((exp(-c*y1))-exp(-c*y2)-exp(-c*h1)+exp(-c*H)); % 目标函数
      %k=-1+fai0*exp(-c*h2);   % k是数值逼近法中的斜率
      n=n+1;
      if n>10
          a=100000000
          break;
      end
end