Quadrati perfetti

program quadrati_perfetti

implicit none

integer :: i,j

j = 1
do i=1,100
        if (i == j**2) then
                print*,i,"Quadrato perfetto!"
                write(unit=10,fmt=*),i,"Quadrato perfetto!"
                j = j+1
        else
                print*,i
                write(unit=10,fmt=*),i
        end if
end do

end program