As @Kalyan suggested you could use the symbolic toolbox if you have it. Otherwise you could just go to the definition of a partial derivative. Evaluate your function to quantify how much it changes when you make small perturbations to the variable you are taking the derivative with respect to. Then the partial derivative is the small change in the function divided by the small change to the input variable. You have to be careful though on the size of the perturbation, too small and the round off causes problems, too big and it isn't a very good approximation. Maybe plot your estimate of your derivative versus perturbation size to get an idea of the behavior. Hopefully has a relatively stable plateau between too small and too big.