T(M, k, Ec, Pr, y) = (A*Ec*Pr/(Ec*Pr-4*A))*( exp(-2*y*sqrt(A)) - exp(-y*sqrt(Ec*Pr)) ) + (y*sqrt(Ec*Pr)/2 + 1)*exp(-y*sqrt(Ec*Pr));
DT = gradient(T(M, k, Ec, Pr, y),Vp); Vp = [M k Ec Pr y]; Vs = [5 1 0.1 5 0.1];
% ANSWER is coming as this for the above approach: 0.14785682021063908572684904500442
% But when I ran the derivative separately and adding them answer is coming as : 0.2014966764864323281894134426085 ( correct ANSWER according to my work as % I have drawn fig which is matching with previous work)
% How to modify or ...