function [h,ros,yy]=eqb(y,n,fai0,c,rosg,rom,lay)
% 沉積物負(fù)荷校正
        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