Hi Vikor,
I somehow don't get "depends" to work. What I try is
- Writing the matrix
(%i1) K_L : matrix([f_L,beta_L,c_Lx],[0,alpha_L*f_L,c_Ly],[0,0,1]);
(%o1) matrix([f_L,beta_L,c_Lx],[0,alpha_L*f_L,c_Ly],[0,0,1])
- For a test, I want to consider f_L, c_Lx, c_Ly as variables and the rest as constants
(%i2) depends([f_L,c_Lx,c_Ly],K_L);
(%o2)
[f_L(matrix([f_L,beta_L,c_Lx],[0,alpha_L*f_L,c_Ly],[0,0,1])),c_Lx(matrix([f_L,beta_L,c_Lx],[0,alpha_L*f_L,c_Ly],[0,0,1])),c_Ly(matrix([f_L,beta_L,c_Lx],[0,alpha_L*f_L,c_Ly],[0,0,1]))]
- Hoping that Maxima understands that it has to differentiate regarding to the three depending variables. It doesn't.
(%i3) diff(K_L);
Maxima encountered a Lisp error:
Error in PROGN [or a callee]: Bind stack overflow.
Automatically continuing.
To reenable the Lisp debugger set *debugger-hook* to nil.
- Trying to explicitly tell, what variables are to differentiate - no luck.
(%i4) diff(K_L, [f_L,c_Ly,c_Lx]);
Non-variable 2nd argument to diff:
[f_L,c_Ly,c_Lx]
-- an error. To debug this try debugmode(true);
- What works is the following:
jacobian(flatten(args(transpose(K_L))),[f_L, c_Lx, c_Ly]); /*works!*/
But how would I do that with "depends"?
Any more tips I could try?
Best regards,
weaker