利用卷绕数判断一点是否含在连续可微的区域D内(Fortran 语言版本