26.07.2013 Views

PDF-format

PDF-format

PDF-format

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

20<br />

løs := {x = 1 −b +<br />

2<br />

√ b2 − 4 a c<br />

}, {x =<br />

a<br />

1 −b −<br />

2<br />

√ b2 − 4 a c<br />

}<br />

a<br />

For at få en enkelt løsning ud af en sekvens af løsninger må vi indeksere i sekvensen:<br />

> løs[1];<br />

> løs[2];<br />

{x = 1<br />

2<br />

−b + √ b2 − 4 a c<br />

}<br />

a<br />

{x = 1 −b −<br />

2<br />

√ b2 − 4 a c<br />

}<br />

a<br />

Med vores generelle løsning på andengradsligingen kan vi nu bruge eval til at få løsninger til forskellige andengradsligninger:<br />

> eval(løs,{a=2,b=-3,c=1});<br />

> eval(løs,{a=2,b=4,c=2});<br />

{x = 1}, {x = 1<br />

2 }<br />

{x = −1}, {x = −1}<br />

Vi kunne selvfølgerlig også får løsningerne til denne andengradsligning ved at bruge solve direkte:<br />

> solve(2*x^2+4*x+2=0,{x});<br />

{x = −1}, {x = −1}<br />

For ligninger af mere end fjerde grad findes der ikke nogen generel løsning og Maple vil derfor ofte give løsninger<br />

i form af RootOf funktionen. Det betyder “løsningerne du leder efter er rødder i polynomiet inde i RootOf”. Maple<br />

indfører selv en hjælpevariabel i polynomiet, den vil som regel hedde _Z, _Z1 eller lignende (alle variable som Maple<br />

indfører starter med et understregningstegn “_”).<br />

Følgende sjettegradsligning har dels en løsning som Maple finder algebraisk og dels fem løsninger som er rødder<br />

i et femtegradspolynomium:<br />

> løs:=solve({4*x^6-12*x^5+2*x^2-9*x+9=0},{x});<br />

løs := {x = 3}, {x = RootOf(4 _Z 5 + 2 _Z − 3, index = 1)},<br />

{x = RootOf(4 _Z 5 + 2 _Z − 3, index = 2)},<br />

{x = RootOf(4 _Z 5 + 2 _Z − 3, index = 3)},<br />

{x = RootOf(4 _Z 5 + 2 _Z − 3, index = 4)},<br />

{x = RootOf(4 _Z 5 + 2 _Z − 3, index = 5)}<br />

For at få numeriske løsninger ud skal vi så bruge evalf:<br />

> evalf(løs);<br />

{x = 3.}, {x = 0.8086092632}, {x = 0.3993814688 + 0.8387514991 I},<br />

{x = −0.8036861004 + 0.6548567008 I}, {x = −0.8036861004 − 0.6548567008 I},<br />

{x = 0.3993814688 − 0.8387514991 I}<br />

Her fik vi atså fire imaginære rødder og en enkelt reel rod ud over roden x=3 som solve fandt algebraisk.<br />

Man bør altid checke sine løsninger ved at substituere tilbage ind i den oprindelige ligning. Dette kan man gøre<br />

med eval.<br />

> lign:=3*x-4=0;<br />

> løs:=solve(lign,{x});<br />

> eval(lign,løs);<br />

lign := 3 x − 4 = 0<br />

løs := {x = 4<br />

3 }<br />

0 = 0<br />

Med solve kan man også løse ligningssystemer. Man angiver et ligningssystem som en mængde af ligninger:

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!