house_c.m

来自「FISMAT accommodates different arithmetic」· M 代码 · 共 11 行

function ok = house_c(ik,x,T)

% ok = house_c(ik,x,T)
% 
a=5*(5/9);
b=-a;
if (ik >= a) & (x(2) < 0.5), ok=1; end;
if (ik <= b) & (x(2) > 0.5), ok=0; end;
if (ik < a)  & (x(2) < 0.5), ok=0; end;
if (ik > b)  & (x(2) > 0.5), ok=1; end;