Bestimmen von Nullstellen
fnt:= 'font = [COURIER, 12]':
Erstes Beispiel
Mit solve erhält man hier offensichtlich (zunächst) nicht alle Lösungswerte:
> plot(f(x), x = -2*Pi..2*Pi, fnt);
Aber:
Eigenschaften der Konstante _Z1:
Originally _Z1, renamed _Z1~:
is assumed to be: integer
Dies ergibt:
> lgn:= 0, Pi*(k + 1/2);
Numerische Berechnung einiger Nullstellen von f(x):
zeros:= zeros union {subs(k = i, evalf(lgn[2]))}:
od:
op(sort(convert(zeros, list)));
Zweites Beispiel
fnt:= 'font = [COURIER, 12]':
> f:= x -> 2^x -
0.95*x^2;
plot(f(x), x = -2..5, y = -4..4, fnt);
Mit fsolve erhält man im Allgemeinen jeweils nur einen Lösungswert:
fsolve(f(x) = 0, x = 2);
fsolve(f(x) = 0, x = 4);
2.197
-0.7823
2.197
3.705
Aber:
> solve(f(x) = 0, x);
2.197, 3.705, -0.7823
Drittes Beispiel
Berechnung einiger Nullstellen von f:
zeros:= {}:
for i from -18 to 10 do
zero:= fsolve(f(x) = 0, x = i):
if not evalb(zero in zeros) then
zeros:= zeros union {zero}:
fi:
od:
liste:= sort(convert(evalf(zeros), list)):
for i from 1 to nops(liste) do
print(liste[i]);
od;
Die Situation in einer kleinen Umgebung von −4,712, bzw. von −7,854:
> plot(f(x), x = -4.714..-4.710, tickmarks = [3,3], fnt);
> plot(f(x), x = -7.855..-7.853, tickmarks = [3,4], fnt);
Viertes Beispiel
f:= x -> (sin(x^4))/(sqrt(x^2+1))-1;
Die Anwendung von solve liefert hier ein falsches Ergebnis:
f(x[0]);
plot(f(x), x = 0..2, y = -2..1, tickmarks = [4,3], fnt);