eqb.m

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

function [h,ros,yy]=eqb(y,n,fai0,c,rosg,rom,lay)
% 沉积物负荷校正
        ros=0.0;
        row=1;
        if n==1
           h=0;
           yy=0;
           return
        end
        for j=2:n
               if (y(j)-y(j-1))==0
                  fai=0;
               else
                  fai=fai0/c*(exp(-c*y(j-1))-exp(-c*y(j)))/(y(j)-y(j-1));
               end
               rob=fai*row+(1-fai)*rosg(lay-n+j);
               ros=ros+rob*(y(j)-y(j-1));
        end
        ros=ros/y(n);
        h=y(n)*(rom-ros)/(rom-row)  ;    %h是校正后的基底深度
        yy=y-(y(n)-h);
end