Литература по Компьютерному практикуму / Maxima.Ильина

¡ë ì ®¡ ¨¬ . •®í ®¬ ¯®¯ë ª

- ¯¨ ì

w:i t grate(f(x),x a b $

to solve for t

changevar(w,x=g(t),t,x);

¢ë§®¢¥

-¨¥ ®¡

訡ª¥ ("Unableg 1( )" ¨ "g 1

(b)") -¥ ¬®£ ¡ë ì ¢ë ¨ «¥-ë).

Ž¤- ª®, ¢¯®«®®¡é-¥ ¤®¯ ¨¬®

[ ]

changevar(w,x=t+c,t,x);

b-c

R

f(t+c) dt

¨

a-c

changevar(w,x^2-t, t, x);

2

b

f(-sqrt(t))

dt

” -ª ¨ï

byparts

-

R

a2

2

-ï¥

- ¥£ ¨ ®¢ -¨¥ ¯® ï¬. •¥ ¥¤ ¯¥ ¢ë¬ ¢ë§®¢®¬ -ª ¨¨ -¥®¡-

¢ë¯®«¤¨¬® § £ §¨ ì ¯ ª

"bypart", ¢ ª® ® ®¬ ®-

®¯ ¥¤¥«¥- :

¢ -¨ï, "u " ¨

- ë ¤

-ë ¨¬¥ ì

¢¨¤: - ¥£ «, ¯¥

-- ï ¨-

¤

load(bypart)$

-¥ ® «¨ ¥ ®¯ ¥¤¥¬«¥--ë© ¨-

«

-¥®¯ ¥¤¥-

0

-

¥ «¨®«¦§ ¢ - | ®-

" v

" (¨ ¯® ì ¥ ï ® ¬ «

"

R

uv

0

= uv

R

u

0

v "). ” -ª ¨ï -¥

«¨èª®¬

€« --£ ®¬£¥® ¨ ¢ ¥£¤

¢®§¢ é ¥ -¥®¯ ¥¤¥«¥--ë©, ¯®í ®¬ -

¡®¥£ ë¨ ®¯®¥

® ®¢

¨

w:integrate(byparts(w,x,

*,

,x)$

,x,a,); b)$

byparts(w,x,f(x),cos(x));

¯ ¨¢¥¤ ª ®¤-®© ¨ ®© ¦¥ ¢ë¤ ¥

R

dx

(f(x))) dx

” -ª ¨ï "integrate"

¥

f(x) sin(x) -

sin(x) (

®®¡ §¨ ì, ® ¨- ¥£ « ® ¯ ®¨§¢®¤-®© ¯ ®¨§-

¢®«ì-®© -ª ¨¨ ¥ ì ¬®¦

-ª ¨ï

90

” -ª ¨ï taylor

¥¯-ë¥ ¤ ®¡¨

’¥©« .

•¥§ «ì ¢ë§®¢

-

ï¥

ª« ¤ë¢ ¥ -ª ¨î ¢

ï ® ®¡ë¬ ¢ë ¦ -¨¥¬ |

" ï¤

¬",

í ® ¢ë ¦¥-¨¥ - ¡¦

¥¦ -ï

¬¥ ª © "/T/"

¯® «¥

"%o". € £ ¬¥- ë -

ª ¢ë: ¢ë

¡-

¤ § §«®¦¥-®; ¯¥ ¥¬¥-- ï, ¯® ª

© ¨¤¥ §« ¦¥-¨¥; ® ª

¨¥¢ª, ¨¨® ® ®©¥ ¬ë

¥ ª« ¤ë¢ ¥¬¬;¥ ª¨¯® 冷ª, ¤® ª® ® ®£®

¨¤ª¥ ¨¨

§«®¦¥-¨¥:

(%i7) ta1:taylor(sin(x),x,0,3)

3

+ . . .

(%i8)

ta2:taylor(cos(x),x,0,6)%o7) /T/ x -

x6

+ . . .

(%o8) /T/ 1 -

x2

+ x

-

x6

•ï¤ë ¬®¦-

2

24

720

¤ £ , ¯ ¨

¤ë¢ ì, ¢ë ¨ ì, ¬-®¦ ì

¤¥« ì ¤ £ -

í ®¬ ® -® ì ª« §«®¦¥-¨ï ¨ ë¢ ¥ ï

¢ ®¬ ¨¥ª¨, - ¯ ¨¬¥

(%i9) ta1/ta2;

(%o9) /T/ x +

x3

+ . . .

ï¤ ).

.¥. §«®¦¥-¨¥ ¨¤ ¤®

¯® 浪

( ® -3

• ¨ í ®¬

ª ¨ ¥ ª¨ ¥ì¥¨¤ì¥£®® 拉 ‹® - ,

® .¥.ì ¤¯¥¯ ¢®ª£®¥ï

taylor(sin(x)/x^3,x,0,5);

6

+

120

-

5040

+ . . .

2

-

1

1

x

x

x

•®«¥¥ ®£®, é¥ ¢ ¥ íª§® ¨ ¥ ª ï ¢®§¬®¦-® ì

1

p]);1

taylor(exp(1/x),[x,0,3,'asym92 2 x2

+

6 x3

+ . . .

+

+

¨§«® ¢ª«î¥ , ¯®-ë¨ï¤ª®í¥«ìª -®, ¤¨¥®¡-®¢-n, 勞®¨®-’«ì®ª¥©«®-§ ©

¬¥¤-ª ¥¨-«ï¥©.m®¢¯….¥•¤§£ 쬥- ë |ª®«¨ï¤¥’¥¢®©

¢®

ª

¨ «¨®¢¥«ï¢ ¤ ®¡-®-ï¤

¨®

«ì- ©

-ª ¨¨ ¬¨¬-¥¥ ¤¨ï, - (¯® ª®«ì-

ª

¨ ¨§¥- ¬¥-

¥«ì ®¯ ¥¤¥«¥-ë

®«¦-® ìî

¤® ®¡é¥£®¡é¨¬® -®¦¨ ¥«ï).

ˆ-묨¨ «¨®í¥¢«ì¬¨, N

=

n + 1) + (m + 1)

1.

ta1:taylor(log(1+x),x,0,6);

+ x3

-

x4

+

x5 -

x6

+ . . .

pade(ta1,3,3)

x -

x2

112 x3

+360 x24

+ 605x

6

pade(ta1,4,2)

[3 x3

+ 36 x2

+ 90 x + 60]

[-

x4

- 12 x3

+

- 150 x2 - 180 x

]

” -ª ¨ï cf

72 x2

240 x + 180

‘®§¤

¥¯- î ¤ ®¡ì,

¯¯ ®ª

î ¤ --®¥ ¢ë ¦ -¨¥. ‚ë ¦¥-¨¥

¤®«¦-® ® ®ï ì ¨§ ¥«ë ¨ ¥«, ª¢ ¤ -ë ® -¥© ¥«ë ¨ ¥« ¨ §- ª®¢

¨ ¬¥ ¨ ¥ ª¨ ®¯¥ ¨©. ‚ë¤ ¨¬¨¨¬î饥 ¢¨¤ ¯¨ ª .

•¥-¨¥¥

¬1.¥

-- ï cflength

®¯ ¥¤¥«ï¥ ª®«¨ ¥ ¢® ¯¥ ¨®¤®¢ ¥¯-®© ¤ ®¡¨. ˆ§- «ì-® -®¢«¥-® §- -

-ª ¨ï cfdisrep

¯ ¥®¡ § ¥ ¯¨ ®ª (ª ª ¯ ¢¨«® ¢ë¤ -ª ¨¨ "cf") ¢ ®¡ ¢¥--® ¥¯- î

”¤

®¡ì.

(sq t(3));

[1,1,2]

cfdisrep(%);

float(%-sqrt(3.0));

1+

1+

-0.065384140902211

93

"•— ¡®ì -ª oat¨©,- ª®¨ «® ¬¨"ë¥

.¯®§¢®«ïî• ® -ª ¥¨¨,è ª®ì ® ë¢-¥¥-¢ë¤¨ï, ¯î¥ ¥® ¨¢¥« -¢ ¢¥ ¬¨§¤-¥«¥

float- ¨ ¥« (allroots, find_root, newton, mnewton).

” -

realroots

-ë¥ ª® -¨ ¯®«¨ ®¬¨ «ì-®£®

- -¨ï

-ª ¨ï,ª® ®

ï ¢ë¤ ¥ ¤

¤

-묨

ª

¨ ¨¥-

¬¨. • ¨

¬ -ª ¨ï ¡®

¥ ¢®¢ ¥

-¥

float-

¬,

-

¢ë¤ ¥ ® ¢¥ ¢ ¥ ¬¨-

®- «ì-ë

. … «¨

¢¥ ¥©¨  參ì

¨®- «ì-ë¥. ’® -® ì í ®©í ®¯¥ ¨¨

§ ¤ ¥ ï ¯¨ ¥«¬¥ ®¬

"ratepsilon"

( ¬.¨ §¤¥« "•¥© ¡®¢¨ ¥«ìoat- ¨ « ¬¨").

ª®í

¨ ¨¥- ë ¢

® -®í¬

- -¨¨ ïî ï float- ¨ « ¬¨, ® ®-¨ ª --

§ ¤

-- ï

rootsepsilon

ª -¥© ¤«ï -ª ¨¨ "realroots".

§ ª § -- î

-

® - ì

¯®

•®

¬®«

¨î § ¤ -® ¤

®¥ §- ¥-

1.0e-7. ˆ¬¥© ¥ ¢

¥ ï

¥ ¯¥ ¥ ¤

ª

¨®--«ì-ë¬ ª®í ¨ï,

¨¥-®£®¬. ’ ª ®

® -® ì

¥è¥-¨ï ¨

¯®¤-«®£®

®¢®«ì-¥©-¨ï float-ª®í ¨ ¨¥-

¬¨ ¡ ¤¥ § ¢¨ ¥¢¨¤ì ¥,é¥ ¨

•í ¥® ¥®¬-¥® ì

¨ ª

ª -

-¥ ¨

¨®¤ª-®£®

¢-¨¥-

¢-¥-¨ï, ª® ®

®¥

¯®«®

´

"ratepsilon".

-- ï

multiplicities

•¯®¥«¥¬¢ë¯®«¥

-¥-¨ï -ª ¨¨ "realroots" ®¤¥ ¦¨ ¯¨ ®ª ª -® ¥© ª® -¥©.

realroots(x^6+1=x^2+x^4);

” -

multiplicities;

x = -1, x = 1]

nroots

ª

[2,

2]

¤¥©

-ë ª

-¥© ¯ «¨-®¬¨ «ì-®-

-ª ¨ï,ª® ®

ï ¢ë¤

£® ¢-¥- ï

¤

¥ ¥«ì®«¨-묨¥ ª¢®í ¨ ¢¨¨¥-¥«ì¬¨, ª® ® ë¥ «®ª «¨§®¢ -ë ¢

ª § --®¬

¨- ¥¢¥©« ¢¨

minf,inf);

nroots(x^6+1=x^2+x^4,0,inf);

2

96

¯¨¥

¢¥-¬ë--¨ï.ë€, ®¤£ -¬®©¥¢ë¤--¥ë¨§¢

¥-|ª-í¨¨©®.

"Š¯¨algsys®¬®ª¥ "¥®è|£®,¥-í¨©¤ ®.¯ ¯¨ª

î®®«ìªï -¢¥ª¤®®¯-¥¦¤®-¨©¥¥¤¨¥¥«è--¥ 먪¥¥

¥ ì ¯¨ ®ª §-

¥¥-¨© ª ¦¤®© ¨§

©

--ë , ®

•® ª ¥è¥-¨© | í ® ¤¢®©¯¨-®©

„«ï ¢-¥-¨©

- «¥¢®© ¯

ìî í -¯¨« ®

¯ ¢ î ì ¬®¦-®

¢«®¦¥¬¥--ë© ¯¨ ®ª.

-

§ ¥ ¢«¥-- ï,¥¢§î¢¨ ¨¬® ¨ ® ¢¨¤

¯ ª ì. -ª ¨¨ "algsys" ¤®

¥ -®©

¢-¥- ©

¢®«ì- ¯¬®¦¥ ¥¬¥ ¢ë§ë¢ ì -

"allr ots

"real--

‹®® £®¨ª®© ¨ª¦¥¬ë¢¥ ì¬

§¢¥ ¢«¥-- ï, ¬®¦¥ , ¢

®ª ¥¨¨¥¤ì, ¢ë§ë¢ ì -ª

roots", "solve". (Š

¨,

- ª® ® ë ¨ ¨ï -ª ¨ï "solve",

«®£¨ª

” -ª ¨ï "algsys" -¥ ª®-¢¥ ¨ ¥ oat- ¨ « , ¢¢®î¤ï騥 ¢ ¨ ¥¬ , ¢ ¨®-

¨î "algsys").

- «ì-ë¥.

algepsilon

•¥ ¥¬¥-- ï

ï

¤«ï

-

"algsys". •® ¬®«-

- § ¤

ì

® ì ¥è¥

¨ §-¥¬ë

-¨î -

«¥- §-

-¨¥ 10 ,

¥ 8 ¢¥ -ë §- ª®¢. •®¢ëè¥-¨¥

„®«¦-® ¨ ®® ¢¥ ¢ ¥

¢¥«¨ ¥

¨î ¯®ª §

¥«ï ¥¯¥ª- ¨¨.

Š ®¦ «¥-¨î, í ® ¯ ¥ ¬¥

8

®¡® ¥ .

-

” -ª ¨ï " algsys" ¨£-® ¨ ¥

¯¥ ¥¬¥-- î "

algepsilon".

--

%rnum list

¯®ª -« ¨ î ï

¨§ ¯ ¥ ¨ª "%r" ¨ ¥«®£® ¨ « , - ¯ ¨¬¥ "%r7".

•¥ ¥¬¢ë§®¥

-ª ¨¨ "algsys"

®

¦¨ ¯¨ ®ª -¥®¯ ¥¤ « --ë ¯

¢ ¤ïé¨ ¢ ¥è¥-¨¥ ¤«ï -¥¤®®¯ ¥¤¥« -- © ¨ ¥¬ë.

ˆ¬¥- í ¨ ¯ ¬¥ ®¢,

3*y+2

=0],[x]);

1], [x = 2]]

);

x=y],[x,y][[x 2, y = 2], [x = 1, y = 1]]

algsys([x^2-y],[x,y]);

[x = %r13, y = %r132 ]]

%rnum_list;

[%r13]

97

‚-¨ «¯¨¥¢®©®ª¯¯¥®¢®©¥¬¥---ëª ì, í"¢ë¤algsys¯ ¢|",î ¯¨- ®ª ì¨ï¬®¦¥"èsolve¥-®¨©®.¯" •®ª-¨ íì.®¬ ¤«ï oat¢- -¨¨©« ,

®¤ï騥 ¢ ¨ ¥¬ , ¢ ¨®- «ì-ë¥. •®í ®¬ ¥¥ ¯®¢¥¢¤¥ -¨¨¥ ¬®¦¥ ¥ §

¢¨ ¥ ì

¢® ¯ «¨ ¬¨¥¥

"ratepsilon" ( ¬. §¤¥« "• ¡®

oat- ¨ « ¬¨").

^2-3*x+2=0],[x]);=

1, x = 2]

solve([sin(x)-1/2],[x]);

to

`solve' is using arc-trig f

get a solunctions.

Some solutions will be lost.

•¥ ¥¬¥-- ï

multiplicities

[x =

%pi

6 ]

®¤¥ ¦¨ ¯¨ ®ª ª -® ¥© ª® -¥©, - ©¤¥--ë -ª ¨¥© "solve"

so ve([x^4-3*x^3+2*x^

]);

multiplicities;

2=0],[x1, x = 2, x = 0]

[1,

1, 2]

„«ï ¨ ¥¬ ¢-¥-¨© ¥è¥-¨¥ | í ® ¤¢®©-®© ¢«®¦¥--ë© ¯¨ ®ª ( ¬. -ª-

¨î "algsys"):

3, y

=

1]]

solve([x+y=4,x-y=2],[x,y]);[[x

(¢ ¤ -®¬ « ¥è¥-¨¥ ®¤-®).

„«ï -¥¤®®¯ ¥¤¥« --ë ¨ ¥¬ ¢ ¥è¥-¨¥ ¢ ®¤ï -¥®¯ ¥¤¥«¥--ë¥ ¯

¢¨¤ "%r4",

¯¥ ¥¬¥-- ï "%rnum list ¯® «¥ ¢ë§®¢

-ª ¨¨ "solve"¬¥®¤ ë-

¦¨ ¨ ¯¨ ®ª ( ¬. -ª ¨î "algsys"):

]]

solve([x^2-y],[x,y]);

[x = %r5, y = %r52

%rnum_list;

[%r5]

®

” -ª ¨ï "solve" ¨¬¥¥ ¤®¢®«ì-®

§¢¥ ¢«¥-- î «®£¨ª . ‚ §

ª®-ª ¥ -®£® ¢¨¤ ¢-¥-¨ï ¨«¨ ¨ ¥¬ë

®-

¢¥¤¥ ¥¡ï ® ¥-ì ¢¨¯®-¨¬®§- ¬¨ ¨

98

‘«è¨¬¥-¥¨î,¤ª¥ ¯ ¨ª¦ «ì ¯¥© ¥¢¬¢¨¤¥¡®--ë¥, -(¨®«¯®£-®¢),ª¨¨ï¡¥"¯®«solve륧-¬"ë¥.-¯ïî„¥¢«ï«¥¥¥¯®¢®¬,磻¥®¢®«ì-¨®¥. -Š®-¡®«ìª®¦¨ï-

"solve" - ¤¥¦-¨¬®

¥ ì¡® ¥ £« ¢-ë¬ ®¡

®§ ®¬ ¤«ï ¥

¢- - © ¨«¨

ª ® ë¥

¥¢¨¤

® ìî

¨¬¥î ¥è¥-¨¥. ‚

¬ «

¥ « £¨

-¥ ¦-ë. ‚

¡ «¥¥ ®«¨- ¨¢¨¥ ¢®«ì-ë

ï

¥è¥-¨¥ ®¡ë -í ®¯®« ¨ ì -¥¢®§¬®¦-®,

- § -

— ®¡ë ¯ ®¤¥¬®- ¨ ®¢ ì, - ª

® ¯ ¨ ¤«¨¢ë¬ ¬®¦¥ ¡ë ì ¯®¢¨¥¤ -¥¨¬,¥

¢¨ ¨¬® ¨-ª ¨¨ "

", ¯ ¨¢¥¤¥¬ -¥ ª®«ìª

¯ ¨¬¥ ®¢:

® §- ¥-¨ï « £®¢ë ¯¥ ¥¬¥--ë .

8 32 %i + 0.

solve([x^5+y=7,x=y],[x,y]);1

36

1 36 164

32 %i + 0.5084694089,

6169408

- 1.

%

y =

- 1.

%i,

-5084694.9241 81109 %i -2138763.

33450],

0.

.9241

089%i -

.

.

450,

[x =

-

1

%i -

832450],

1.41081382385, y

36861648= .41081382385]]

solve(x^5+x-7);

0.92418 11

09

%i - .

[0 = x

5

+ x -

solve(x^6+x-7.0000000123456789);

`rat' replaced -7].00000001234568

[0 = x

6

by -7/1 = -7.0

+ x - 7]

” -

solve([sin(x)-y,sin(x)+y=1][,[x,y]);

¨ï

eliminate

¢- -

ª § --ë¥ ¯¥ ¥¬¥--ë¥. Ž

¢-¥-

¨ ª«î

¥ ¨§ ¨ ¥¬ë

-¨ï ¯ ¨¢®¤ï ï ª ¢¨¤ - «¥¢®¨©¯ ¢®©

ìî, ª® ®

ï ®¯ ª

. ” -ª ¨ï

"eliminate" ª®-¢¥ ¨ ¥ oat- ¨ «

¢ ®¤ï騥 ¢ ¨ ¥¬ , ¢ ¢è¨¨®¥-ï«ì-ë¥.

-y z=2,x+y-z=3],[z]);

- 2), - 2 y - 1]

y+z=1,

[2 (y +

eliminate([x-y+z=2,x+y-z=3],[x,y]);

99