ODE Solver through Galerkin Method

[APPROX,EXAC,ERR] = ODEGALERKIN(POLY,BC,N) solves Ordinary Differential Equations (ODE) through Galerkin method, by inserting the characteristic polynomial matrix "POLY", boundary conditions "BC" and the finite quantity of approximative base functions "N". Outputs of the program are the approximative solution "APPROX", the analitic solution "EXAC" and the percentage error "ERR" (%). A plot of the approximative and analitic solution is also shown.