I am working on a program to find the value of the current in a coil. This value satisfies the following equation:

y'=sin(2t)-[(e

^{y}-1)/(e

^{y}+1)]

which is of the form y'=f(t,y)

I know that in order to solve this I need to use the trapezoidal method to solve a differential equation, the formula is:

y

_{n+1}=y

_{n}+.5*h(f(t

_{n},y

_{n})+f(t

_{n+1},y

_{n+1}) where h=t

_{n+1}-t

_{n}

I have found examples of the standard trapezoidal method but I do not think they will help because of the difference in the formulas. Can anybody give me some guidance here?