From: Robert D. <rob...@gm...> - 2007-12-22 00:12:53
|
On Dec 20, 2007 9:07 AM, <cyr...@uk...> wrote: > f(x):=x$ > g(A):=romberg(f(q),q,0,A)$ > phi(u):=g(u)-0.2$ > find_root(phi,u,0,4); > f(x):=x$ > g(A):=romberg(f(q),q,0,A)$ > find_root((g(u)-0.2),u,0,4); > f(x):=x$ > g(A):=romberg(f(q),q,0,A)$ > phi(u,v):=g(u)-v$ > psi(v):=find_root(phi(u,v),u,0,4); > psi(0.2); Cyril, I apologize for writing English here. I believe these examples all yield the expected result with the most recent version of Maxima (namely, the Maxima 5.14.0 release candidate, which you can obtain from Sourceforge). The evaluation of the romberg function was modified after Maxima 5.13.0. Hope this helps, Robert Dodier PS. Here is what I see: (%i1) build_info (); Maxima version: 5.13.99rc1 Maxima build date: 21:32 12/4/2007 host type: i686-pc-mingw32 lisp-implementation-type: GNU Common Lisp (GCL) lisp-implementation-version: GCL 2.6.8 (%o1) (%i2) f(x):=x$ (%i3) g(A):=romberg(f(q),q,0,A)$ (%i4) phi(u):=g(u)-0.2$ (%i5) find_root(phi,u,0,4); (%o5) 0.63245553203368 (%i6) find_root((g(u)-0.2),u,0,4); (%o6) 0.63245553203368 (%i7) phi(u,v):=g(u)-v$ (%i8) psi(v):=find_root(phi(u,v),u,0,4); (%o8) psi(v) := find_root(phi(u, v), u, 0, 4) (%i9) psi(0.2); (%o9) 0.63245553203368 |