def riccati_jn(n,x):
if not (isscalar(n) and isscalar(x)):<EOL><INDENT>raise ValueError("<STR_LIT>")<EOL><DEDENT>if (n != floor(n)) or (n < <NUM_LIT:0>):<EOL><INDENT>raise ValueError("<STR_LIT>")<EOL><DEDENT>if (n == <NUM_LIT:0>):<EOL><INDENT>n1 = <NUM_LIT:1><EOL><DEDENT>else:<EOL><INDENT>n1 = n<EOL><DEDENT>nm,jn,jnp = specfun.rctj(n1,x)<EOL>return jn[:(n+<NUM_LIT:1>)],jnp[:(n+<NUM_LIT:1>)]<EOL>
Compute the Ricatti-Bessel function of the first kind and its
    derivative for all orders up to and including n.