order of convergence of a numerical method

I have computed the approximate solution of a differential equation applying trapezoidal rule and also have the exact solution. Now I need to find the order of convergence. My question is how to code it in Matlab to find order of convergence of a numerical method?