2017-04-05 66 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)