Variational formulation based on duality to solve partial differential equations: Use of B-splines and machine learning approximants