Sekantenverfahren Maple < Maple < Mathe-Software < Mathe < Vorhilfe
|
| Status: |
(Frage) überfällig  | | Datum: | 20:22 Do 19.01.2012 | | Autor: | hubbel |
| Aufgabe | | http://www.myimg.de/?img=Unbenannt72041.jpg |
Hallo Leute,
und zwar soll ich das Sekantenverfahren mit Maple programmieren. Habe etwas gegooglet und bin auf folgenden Code gestoßen, der ähnlich dem Sekantenverfahren ist, sich aber auf Regula falsi bezieht. Habe versucht ihn abzuändern, aber irgendwie bekomme ich nicht das passende heraus. Habe auch schon eigene Versuche gestartet, aber es ist alles gescheitert, bin ziemlich verzweifel, bräuchte dabei Hilfe.
nullstelle := proc (g, ai, bi)
local f,fa,fb, a, b, x, i;
f := g;
a := ai;
b := bi;
fa := f(a);
fb := f(b);
if fa*fb>0 then ERROR(`Gleiches Vorzeichen der Funktionswerte`)
end if;
x := a-fa*(b-a)/(fb-fa);
for i to 20 do
printf("a=%f < x=%f < [mm] b=%f,\t [/mm] f(a)=%f, [mm] f(b)=%f\n", [/mm] a,x,b,fa,fb);
if 0 < f(b)*f(x) then b := x else a := x
end if;
fa := f(a);
fb := f(b);
x := a-fa*(b-a)/(fb-fa);
end do;
convert([a, b, x], list);
end proc;
Danke schonmal...
|
|
| |
|
| Status: |
(Mitteilung) Reaktion unnötig  | | Datum: | 17:20 Fr 20.01.2012 | | Autor: | matux |
$MATUXTEXT(ueberfaellige_frage)
|
|
|
|
