Can someone confirm my differentiation code is correct?

I need to differentiate P(a,u,EA1,EA2) with respect to u. I get the code to run, but im not sure the answer is correct since it gives me a long answer.

syms L10(a) L1(a,u) epsilon1(a,u) L20(a) N1(a,u,EA1) N2(a,u,EA2) P(a,u,EA1,EA2)
L10(a) = 2*a;
L20(a) = sqrt(5)*a;
L1(a,u) = sqrt((2*a)^2+u^2);
L2(a,u) = sqrt((2*a)^2+(a-u)^2);
epsilon1(a,u) = (L1(a,u)-L10(a))/L10(a);
epsilon2(a,u) = (L2(a,u)-L20(a))/L20(a);
N1(a,u,EA1) = EA1*epsilon1(a,u);
N2(a,u,EA2) = EA2*epsilon2(a,u);
P(a,u,EA1,EA2) = (u/(L1(a,u)))*(N1(a,u,EA1))-((a*u)/(L2(a,u)))*(N2(a,u,EA2));
df = diff(P(a,u,EA1,EA2),u)

[darova]

Why are you surprised?

Function you want to differentiate

You didn't attach source formula or equations. How can we verify if it's correct?

[Ruslan Khoguev]

I meant does the code make sense? Am I differentiating P(a,u,EA1,EA2) with respect to u correctly?

[Walter Roberson]

Yes, that is correct syntax for differentiating the expression returned by P(a,u,EA1,EA2) with respect to u.

Answer 1

[Devineni Aslesha]

The code written to differentiate P (a, u, EA1, EA2) with respect to u is correct. It will give correct result. The reason for the answer to be long might be because of the display of all the terms in the differentiated output.