gusucode.com > 润滑理论的求解程序包括弹流和刚流润滑不成功源码程序 > filmrupture.m

    function [Hr,Hi]=filmrupture(P)  %本函数的功能是油膜破裂位置，都是以轴向对称面处的数据为准

[m,n]=size(P);    
for i=1:m        %此for循环是求解油膜破裂区位置
    for j=1:n
      if(P(i,j)==0)
        Hr=i;
        Hi=j;
      break;
end
end
end