a2d0[]=1.0;
atm12d0[]=1.0;

func real defIC(real s, real t){
  real res=0;
  if (H2d0(s,t)<0.5)
    res=0;
  else{
    if (s>0.7 && t > 0.7)
      res=1;
  }
  return res;
}

b2d0=defIC(x,y);
btm12d0=defIC(x,y);
plot(b2d0,cmm="b IC",wait=true);