I have the function f(x,y,y′)=√1+y′(x)2 from which I want to calculate ∂f∂y as well as ∂f∂y′ and ddx(∂f∂y′).
I tried using the SymbolicRing of Sagemath with
x = SR.var('x')
y = function('y')
f = sqrt(1+y(x).diff(x,1)^2)
which only allows me to calculate dfdx. So I tried with the following sample
x,y = SR.var('x,y')
y = function('y')(x)
y1 = y.diff(x,1)
f(x,y,y1) = sqrt(1+y1^2)
fy1 = f.diff(y1)
fy = f.diff(y)
fxy1 = fy1.diff(x)
which works for calculating ∂f∂y and ∂f∂y′, but not for ddx(∂f∂y′). Because fxy1 returns 0 when it should return y″.
I'd really appreciate any idea. Thanks in advance!