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.

lign:={u+v+w=1, 3*u+v=3, u-2*v-w=0};<br />

lign := {u + v + w = 1, 3 u + v = 3, u − 2 v − w = 0}<br />

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

løs := {v = 3 −2 4<br />

, w = , u =<br />

5 5 5 }<br />

Vi kan igen bruge eval til at checke vores løsning:<br />

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

{0 = 0, 3 = 3, 1 = 1}<br />

Man kan også med fordel bruge eval til at få en enkelt variabel “pakket ud” ud af løsningsmængden:<br />

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

Man kan også bruge solve til at løse uligheder:<br />

4<br />

5<br />

> solve(x^2+4>8,x);<br />

RealRange(−∞, Open(−2)), RealRange(Open(2), ∞)<br />

Ovenstående er Maples måde at skrive de to intervaller [ −∞,-2[ og ]2, ∞] på. Hvis vi beder om at få løsningen på<br />

den alternative form ved at pakke x ind i en mængde får vi løsningen som to uligheder hvilket er noget mere læseligt:<br />

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

{x < −2}, {2 < x}<br />

Maple kan også løse systemer der indeholder uligheder. Løsningerne som solve finder opfylder alle ligninger og<br />

uligheder i det system man angiver:<br />

> solve({x^2-9=0,x>0},{x});<br />

{x = 3}<br />

Hvis Maple ikke kan finde en løsning får man ikke noget resultat.<br />

> solve({x4},{x});<br />

Hvor Maple ikke kan løse ligningen algebraisk giver solve en RootOf løsning :<br />

> solve(x^3=cos(10*x),x);<br />

1<br />

10 RootOf(−1000 cos(_Z) + _Z3 , label = _L2)<br />

Vi kan igen få en konkrete løsning ud af RootOf løsningen med evalf:<br />

> evalf(%);<br />

0.1566948942<br />

Men med evalf får vi kun en enkelt løsning ud, selv om der måtte være flere (i dette tilfælde er der syv, hvilket<br />

man kan forvisse sig om ved at tegne graferne for de to funktioner med plot, se afnittet om grafer for funktioner).<br />

Maple kan løse visse ligninger med periodiske funktioner hvor perioderne er indbyrdes lineære:<br />

> solve(sin(x)=cos(3*x),x);<br />

− 1 3<br />

π,<br />

4 4<br />

> solve(sin(x)=cos(1+x),x);<br />

1<br />

5<br />

π, π, −7 π, −3 π,<br />

8 8 8 8 π<br />

cos(1)<br />

arctan(<br />

1 + sin(1) )<br />

Nogle gange må Maple imidlertid helt give op og i disse tilfælde giver solve ikke nogen løsning, selv om der skulle<br />

være en:<br />

> solve(sin(x)=cos(x^2),x);<br />

Når man ikke kan finde løsninger algebraisk med solve kan man forsøge at finde dem numerisk med fsolve:<br />

> fsolve(sin(x)=cos(x^2),x);<br />

0.8493688624<br />

21

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

Saved successfully!

Ooh no, something went wrong!