roessler_analytical_fp_rsp Subroutine

    subroutine roessler_analytical_fp_rsp(fp1, fp2)
        class(state_vector_rsp), intent(out) :: fp1, fp2

        real(sp) :: d

        d = sqrt(c_sp**2 - 4*a_sp*b_sp)

        fp1%x = ( c_sp - d)/ 2
        fp1%y = (-c_sp + d)/(2*a_sp)
        fp1%z = ( c_sp - d)/(2*a_sp)

        fp2%x = ( c_sp + d)/ 2
        fp2%y = (-c_sp - d)/(2*a_sp)
        fp2%z = ( c_sp + d)/(2*a_sp)
    end subroutine roessler_analytical_fp_rsp