Frobenius method