This code solves the parallel RLC circuit using GSL for differential equation solver. The system parameters are L=47mH, C=47 nF, Is=5 A, and R=20 kohm. See the MIT OpenCourseWare text for more details.
Use the following command to compile.
gcc rlc-parallel.c -o rlc-parallel -lgsl -lgslcblas -lm