I would like to write a code in R for POWER of smooth Neyman test. I generate n observations and repeat it 10 000 times.

Function should view as:

neyman <- function(n,k){

..........

............ }