最近拿了个作业，其中一题是要写用FORTRAN写一个判断是否素数的program,但是我还真不知道怎样写，连公式都不知道，，，，，求高手解惑。

还有一题是按照公式写程序。。写了个出来但不知道那里错了。。。运行的时候不停显示错误。。。请各位大侠帮我检查下，或者教教小弟。这是原题。这是我写的，最近才刚学所以

PROGRAM Trigonometric
IMPLICIT NONE
REAL::x, arcsin_x
INTEGER::n, i, a, b
WRITE(*, *)"please enter the x"
READ(*, *)x
WRITE(*, *)"please enter the n"
READ(*, *)n

IF(x<=1.and.x>=1-2)THEN
WRITE(*, *)"x in the limit"
ELSE
WRITE(*, *)"x not in the limit"
END IF
b=0
a=0
x=0
DO i=1, 2*n
a=a*i
END DO
DO i=1, n
b=b*i
arcsin_x=x+(a/((4**n)*(b**2)*(2*n+1)))*(x**(2*n+1))
END DO

WRITE (*, *)"the answer is", arcsin_x
END PROGRAM Trigonometric