We present some numerical methods to solve control problems in the coefficients where the cost functional may depend on the gradient of the state non linearly. The main difficulty comes from the fact that the relaxed functional cost is not explicitly known. We prove some convergence results just using an upper or a lower approximation of this relaxed functional.