differentiation.m

clear all
syms a b w l r k Q x;
f1=a^2+4*(w*l+r*k-r*k*x)*(b-a)/((1-x)*Q/x);
f=symfun(0.5*(a+...
    (a^2+4*(w*l+r*k-r*k*x)*(b-a)...
    /((1-x)*Q/x))^0.5),[a b w l r k Q x]);
g=diff(f,x);

%{ 
result 
g(a, b, w, l, r, k, Q, x) =
 
-((x*(a - b)*(4*k*r + 4*l*w - 4*k*r*x))/(Q*(x - 1)^2) - 
((a - b)*(4*k*r + 4*l*w - 4*k*r*x))/(Q*(x - 1)) + 
(4*k*r*x*(a - b))/(Q*(x - 1)))/(4*(a^2 + (x*(a - b)*(4*k*r + 4*l*w - 4*k*r*x))/(Q*(x - 1)))^(1/2))
%}