Merci à tous, j'ai réussi à installer la dernière version de Maxima sur
Windows XP.
Maintenant, j'ai une autre question technique.
Le programme suivant qui fait afficher le graphe d'une certaine fonction N
calculée avec romberg ainsi que celui de sa fonction réciproque marche sous
Maxima 5.14.0 :
kill(all)$
f(x):=x$
N(A):=romberg(f(q),q,0,A)$
wxplot2d([N(x)], [x,0,1])$
phi(u,v):=N(u)-v$
psi(v):=find_root(phi(u,v),u,0,4)$
wxplot2d([psi(x)], [x,0,1])$
Le programme suivant qui affiche le graphe d'une certaine fonction également
notée N et calculée en utilisant romberg marche tout aussi bien sous Maxima
5.14.0 :
kill(all)$
load(romberg)$
r:-0.5$
a:0.1$
t:10$
x2_0:0.03$
x1_0:0.05$
Q(x):=sqrt(x)-sqrt(x1_0)$
argch(x):=log(x+sqrt(x^2-1))$
G(x):=(sqrt(x^2+1-r^2)-r*x)/(1-r^2)$
H(x):=argch(G(x))$
rmin(q):=H(((a*q)/x2_0)+r)$
F(x):=sqrt(x^2+1-r^2)$
sgn(x):= if x<0 then (-1) else 1$
PSI(x):=(cosh(x)+r.sinh(x)).sqrt(F(r.cosh(x)+sinh(x)))/(F(r.cosh
(x)+sinh(x))^2)$
f(x):= float(PSI(sqrt(a^2*t)*x)*exp(-x^2/2)/sqrt(2*%pi))$
B(x):=sgn(x-x1_0)*rmin(Q(x))/sqrt(a^2*t)$
N(x):=romberg(f(u),u,-10,B(x))$
wxplot2d( [N(x)], [x, 0.03,0.08] )$
Le problème est que je n'arrive pas à afficher le graphe de sa fonction
réciproque en faisant comme précédemment.
Si je rentre le code suivant :
phi(u,v):=N(u)-v$
psi(v):=find_root(phi(u,v),u,0,0.1)$
wxplot2d([psi(x)], [x,0.3,0.8])$
Je n'obtiens rien. Maxima semble calculer en vain. Au bout d'une minute
environ, j'obtiens un rectangle rayé avec écrit "erreur" au milieu...
Quelqu'un pourrait-il m'expliquer pourquoi et me dire ce qu'il faut faire
pour voir le graphe de la fonction inverse de N comme je le voudrais ?
Cordialement,
Cyril
-------------------------------------------------------------------------
This SF.net email is sponsored by: Microsoft
Defy all challenges. Microsoft(R) Visual Studio 2005.
http://clk.atdmt.com/MRT/go/vse0120000070mrt/direct/01/
_______________________________________________
Maxima-lang-fr mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/maxima-lang-fr
