Function 函数(函数)在极大值中的函数
我想在极大值中定义函数的函数(可能是函数)。例如,我想在代码中定义T:Function 函数(函数)在极大值中的函数,function,user-defined-functions,maxima,Function,User Defined Functions,Maxima,我想在极大值中定义函数的函数(可能是函数)。例如,我想在代码中定义T: f(r,GN,M)::=(2*GN*M^2)*(2*GN*M/r-(2*GN*M/(2*r))^2); X(r,GN,M) ::= 2*GN*M/(2*r)+2*GN*M/r*(1-2*GN*M/(4*r)); T(r,GN,M) ::= diff(f(r,GN,M),r)+X(r,GN,M); 但是我不知道代码。在maxima中,使用表达式很方便 定义两个表达式 (%i1) display2d: false $ (%i2
f(r,GN,M)::=(2*GN*M^2)*(2*GN*M/r-(2*GN*M/(2*r))^2);
X(r,GN,M) ::= 2*GN*M/(2*r)+2*GN*M/r*(1-2*GN*M/(4*r));
T(r,GN,M) ::= diff(f(r,GN,M),r)+X(r,GN,M);
但是我不知道代码。在maxima中,使用表达式很方便 定义两个表达式
(%i1) display2d: false $
(%i2) f: (2*GN*M^2)*(2*GN*M/r-(2*GN*M/(2*r))^2) $
(%i3) X: 2*GN*M/(2*r)+2*GN*M/r*(1-2*GN*M/(4*r)) $
(%i4) T(f, X) := diff(f, r) + X $
(%i8) subst([r=1, GN=2, M=3], dfpX);
(%o8) 2142
定义一个对两个表达式进行操作的函数
(%i1) display2d: false $
(%i2) f: (2*GN*M^2)*(2*GN*M/r-(2*GN*M/(2*r))^2) $
(%i3) X: 2*GN*M/(2*r)+2*GN*M/r*(1-2*GN*M/(4*r)) $
(%i4) T(f, X) := diff(f, r) + X $
(%i8) subst([r=1, GN=2, M=3], dfpX);
(%o8) 2142
叫它。结果是一个表达式
(%i5) dfpX: T(f, X);
(%o5) (2*GN*M*(1-(GN*M)/(2*r)))/r+(GN*M)/r+2*GN*M^2
*((2*GN^2*M^2)/r^3-(2*GN*M)/r^2)
您可以基于dfpX
(%i6) define(f(r, GN, M), dfpX);
(%o6) f(r,GN,M):=(2*GN*M*(1-(GN*M)/(2*r)))/r+(GN*M)/r
+2*GN*M^2*((2*GN^2*M^2)/r^3-(2*GN*M)/r^2)
叫它
(%i7) f(1, 2, 3);
(%o7) 2142
但你可以继续使用表达式
(%i1) display2d: false $
(%i2) f: (2*GN*M^2)*(2*GN*M/r-(2*GN*M/(2*r))^2) $
(%i3) X: 2*GN*M/(2*r)+2*GN*M/r*(1-2*GN*M/(4*r)) $
(%i4) T(f, X) := diff(f, r) + X $
(%i8) subst([r=1, GN=2, M=3], dfpX);
(%o8) 2142
将
T
定义为函数,而不是宏,以便它计算其参数
(%i) T(r,GN,M) := diff(f(r,GN,M),r)+X(r,GN,M)$
现在在t
上定义一个函数t
,该函数引用”
其参数(以便t
知道要区分哪个变量),然后使用ev
计算结果%%
(%i) t(r,GN,M) := (T('r,'GN,'M),ev(%%))$
(%i) t(1,2,3);
(%o) 2142