一个关于“同构数”的Fortran程序的问题

题目要求1-1000中的同构数
program ex309
implicit none
integer x,n
do x=1,1000
do n=1,4
if((x/10**n)==0)then
goto 1
endif
enddo
1 if(mod(x**2-x,10**n)==0)then
print*,x
endif
enddo
end
这是我个人编写的程序~貌似和书上的答案有所出入~
书上的答案是
program ex309_1
implicit none
integer m,l,n,j
do m=1,1000
k=m*m
l=1
n=0
do while((n<4)).and.(l==1))
n=n+1
j=10**n
if(m/j==0)l=0
enddo
if(mod(k-m,j)==0)
print*,"m=",m,"m**2=",k
enddo
end
但是书上的答案输入电脑里却出现错误
ex309_1.F95(8) : error 58 - Unpaired right bracket(s)
ex309_1.F95(13) : error 363 - Expected label/statement/THEN after IF statement, inserting THEN
ex309_1.F95(15) : error 763 - Nesting error - this END DO statement is matched with the IF construct starting on line 13
ex309_1.F

请教一个同构数的C程序编写方法? 用C＋＋编一个1～1000的同构数 一个数恰好等于它的平方数的右端，这个数称为同构数，如，5的平方等于25，求1到1000之间的所有同构数， 什么是同构数 所谓同构数 到底什么是“同构数” C语言编程100以内所有的同构数 用C++编1到1000的同构数 求1~9999之间的全部同构数 用PASCAL 求同构数它出现在它的平方数的右端。例如：5的平方是25，5就是同构数，25也是构数。