www.gusucode.com > 界面看成高一维空间中某一函数ψ源码程序 > levelsetcode/Heaviside.m
function H = Heaviside(phi,epsilon) % Heaviside(phi,epsilon) compute the smooth Heaviside function % % created on 04/26/2004 % author: Chunming Li % email: li_chunming@hotmail.com % Copyright (c) 2004-2006 by Chunming Li H = 0.5*(1+ (2/pi)*atan(phi./epsilon));