File: newton.mc

package info (click to toggle)
maxima 5.6-17
  • links: PTS
  • area: main
  • in suites: woody
  • size: 30,572 kB
  • ctags: 47,715
  • sloc: ansic: 154,079; lisp: 147,553; asm: 45,843; tcl: 16,744; sh: 11,057; makefile: 7,198; perl: 1,842; sed: 334; fortran: 24; awk: 5
file content (9 lines) | stat: -rw-r--r-- 196 bytes parent folder | download | duplicates (2)
1
2
3
4
5
6
7
8
9

NEWTON(EXP,VAR,X0,EPS):=
	BLOCK([XN,S,NUMER],
	NUMER:TRUE,
	S:DIFF(EXP,VAR),
	XN:X0,
  LOOP, IF ABS(SUBST(XN,VAR,EXP))<EPS THEN RETURN(XN),
	XN:XN-SUBST(XN,VAR,EXP)/SUBST(XN,VAR,S),
	GO(LOOP) )$