How do we get the final formula of the Bernoulli number? I was trying to understand Bernoulli numbers. When I googled, I found this link. 
It starts by saying that, The Bernoulli numbers are defined via the coefficients of the power series expansion of 
$\frac{t}{e^{t}-1}$,
then they write the expansion,for $m \geq 0$.
$\frac{t}{e^{t}-1} = \sum_{m=0}^{\infty} \frac{B_m}{m!}t^{m}$ (1)
Then, the tutorial says that we multiply both sides by $(e^{t} -1)$ and get 
$B_0=1$, 
$B_m= -\frac{1}{m+1}\sum_{k=0}^{m-1}\binom{m+1}{k}B_k$. (2)
How is this happening? 
I tried doing this with hand and I get a recursive relation, that looks like, 
$\frac{t}{e^{t}-1} = \sum_{m=0}^{\infty} \frac{B_m}{m!}t^{m}$
$\implies t = e^{t}\sum_{m=0}^{\infty} \frac{B_m}{m!}t^{m} - \sum_{m=0}^{\infty} \frac{B_m}{m!}t^{m}$
$\implies t = e^{t}\sum_{m=0}^{\infty} \frac{B_m}{m!}t^{m} - \frac{t}{e^{t}-1}$
What next?
How do I get to (2)? 
I am new to number theory, please excuse me if the doubt is very basic. Any kind of help will be appreciated. :)

With recurrences of this type I like the representation in a matrix format. Note, that the binomial-coefficients in your recursive formula occur in the manner of the lower triangular Pascal-matrix ( = "P") and that we need a modified/trimmed version (= "Q" ) of it.         
The following is the multiplication-scheme for the Bernoulli-numbers [update] corrected :
        |                               |      1  |
        |                               |   -1/2  |  vector B of
        |                               |    1/6  |  Bernoulli-numbers 
        |  -dZ * Q * B = B'             |      0  |
        |                               |  -1/30  |
        |                               |      0  |
        |                               |   1/42  |
        |                               |      0  |
     -  +  -  -   -   -   -   -   -  -  +      -  +
    -Z  |  modified Pascalmatrix Q      |         |
     -  +  -  -   -   -   -   -   -  -  +      -  +
    -1  |  .  .   .   .   .   .   .  .  |      0  |  reduced (first entry =0)
  -1/2  |  1  .   .   .   .   .   .  .  |   -1/2  |  vector B' of 
  -1/3  |  1  3   .   .   .   .   .  .  |    1/6  |  Bernoulli-numbers 
  -1/4  |  1  4   6   .   .   .   .  .  |      0  |  
  -1/5  |  1  5  10  10   .   .   .  .  |  -1/30  | 
  -1/6  |  1  6  15  20  15   .   .  .  |      0  | 
  -1/7  |  1  7  21  35  35  21   .  .  |   1/42  |
  -1/8  |  1  8  28  56  70  56  28  .  |      0  |
     -  +  -  -   -   -   -   -   -  -  +      -  +

This scheme illustrates the matrix-product $$ -\,^dZ \cdot Q \cdot B =  B' $$        where $\,^d Z = \operatorname{diagonal}([1,1/2,1/3,...]) $ .
This is an explication of your formula (2) where $m$ is the matrix-row-index beginning at $0$ .
Of course: so far it is no proof, just an illustration how to understand the ingredients of formula (2). However, manipulations at this formula can be made to related it to known, simpler and proven relations between binomial-coefficients and Bernoulli-numbers, if that is what you want.          
(P.s.: you might also be interested in this small treatize on exactly that subject / method)
