Incorrect derivatives when one of the variables is a function of the differential variable

[mgreminger]

See this example sheet that reproduces the issue.

[mgreminger]

Reopening as unresolved, seems to only happen if the differential variable is combined with another parameter. Works fine for d/dx(y)=20
https://engineeringpaper.xyz/rcyUDavPWn6qDuB9zQqcRt

[mgreminger]

Seems to be a side effect of Sympy's sub behavior with derivatives:

In [1]: import sympy as sp

In [2]: sp.__version__
Out[2]: '1.11.1'

In [3]: x,y,z = sp.symbols('x,y,z')

In [4]: sp.Derivative(x*y, x)
Out[4]: Derivative(x*y, x)

In [5]: test = sp.Derivative(x*y, x)

In [6]: test.subs(y, 20*x)
Out[6]: Subs(Derivative(x*y, x), y, 20*x)

In [7]: test2 = sp.Derivative(y, x)

In [8]: test2.subs(y, 20*x)
Out[8]: Derivative(20*x, x)