2017-04-05 61 views
0

假设我有以下表述:与价值的替代变量,但不评价

(%i1) (8*x)*(log(x)/log(10)); 

(%i2) X^2; 

现在,因为我想找出恒定值,我可以选择作出声明%i1 is O(%i2)真的,我评估他们在像这样一个循环:

for a:1 thru 10 do print(%i1, "=", ev(%i1, x=a), %i2, "=", ev(%i2, numer, x=a)); 

输出是:

8 x log(x)   2 
---------- = 0.0 , x = 1 
log(10) 

8 x log(x)      2 
---------- = 4.816479930623698 , x = 4 
log(10) 

8 x log(x)      2 
---------- = 11.45091011327189 , x = 9 
log(10) 

8 x log(x)      2 
---------- = 19.26591972249479 , x = 16 
log(10) 

8 x log(x)      2 
---------- = 27.95880017344075 , x = 25 
log(10) 

8 x log(x)      2 
---------- = 37.35126001841489 , x = 36 
log(10) 

8 x log(x)      2 
---------- = 47.32549024079837 , x = 49 
log(10) 

8 x log(x)      2 
---------- = 57.79775916748438 , x = 64 
log(10) 

8 x log(x)      2 
---------- = 68.70546067963139 , x = 81 
log(10) 

8 x log(x)   2 
---------- = 80.0 , x = 100 
log(10) 

余万吨至使输出更容易眼球,是这样的:

8 1 log(1)   2 
---------- = 0.0 , 1 = 1 
log(10) 

8 2 log(2)      2 
---------- = 4.816479930623698 , 2 = 4 
log(10) 

8 3 log(3)      2 
---------- = 11.45091011327189 , 3 = 9 
log(10) 

[snip] 

8 10 log(10)   2 
---------- = 80.0 , 10 = 100 
log(10) 

我怎么能告诉Maxima在循环的每次迭代替代的a的价值x不计算表达式?

我查过手册,但没有发现任何看似相关的东西。

回答

1

在千里马很多操作是由一个名为“简化”的过程,这意味着将身份做一个“简单”的表达进行。例如。 1 + 1简化为2,sin(0)简化为0等。

为了得到您想要的效果,我们必须一般禁用简化,以便评估表达式而不是简化表达式。但要获得数值,我们需要为这些结果启用简化。

这是要做的事情。

(%i16) simp : false $ 
(%i17) for x in [1,2,3,4,5] 
     do print (ev(%i1) = ev(%i1, simp, numer), ev(%i2) = ev(%i2, simp)); 
     log(1)   2 
(8 1) (-------) = 0.0 1 = 1 
     log(10) 
     log(2)      2 
(8 2) (-------) = 4.816479930623698 2 = 4 
     log(10) 
     log(3)      2 
(8 3) (-------) = 11.4509101132719 3 = 9 
     log(10) 
     log(4)      2 
(8 4) (-------) = 19.26591972249479 4 = 16 
     log(10) 
     log(5)      2 
(8 5) (-------) = 27.95880017344075 5 = 25 
     log(10) 
(%o17)        done 

注意,我写for x in [1, 2, 3, 4, 5] ...,而不是for x:1 thru 5 ...。这是因为后者使用算术,这需要简化。尝试两种方式,我想你会看到不同之处,我相信它非常有启发性。

nota bene我已经使用了与%i1%i2相同的值。

1

使用“空”的功能:

(%i1) display2d: false $ 
(%i2) prefix("") $ 
(%i3) almost_subst(a, x, e):= subst(""(a), x, e) $ 
(%i4) almost_subst(10, x, 8*x*log(x)/log(10)); 
(%o4) (8* 10*log(10))/log(10)