Interpolation de NEWTON-TCHEBYCHEV
Bonjour,
j'ai ecrit ce programme sous sage, en il me dit que 'p' n'est pas définie.
je ne suis pas trés douée en prgrammation et ne comprend pas pourquoi j'ai de type d'erreur. Merci d'avance pour votre aide
"def NEWTONTCHEBYCHEV(f,t,a,b,Ndigits,N):
corps = RealField(Ndigits); R.<X> = corps[]
T = TCHEBYCHEV(R,N+1)
ZEROS = T.roots(corps)
zeros = [ZEROS[k][0] for k in range(N+1)]
aa = corps(a); bb = corps(b)
Lx = [aa + ((bb-aa)/2)*(zeros[k]+1) for k in range(N+1)]
Ly = [f.subs({t:Lx[k]}) for k in range(N+1)]
for j in range(N):
p = NEWTONAITKEN(R,Lx,Ly)
return p"
j'ai ecrit ce programme sous sage, en il me dit que 'p' n'est pas définie.
je ne suis pas trés douée en prgrammation et ne comprend pas pourquoi j'ai de type d'erreur. Merci d'avance pour votre aide
"def NEWTONTCHEBYCHEV(f,t,a,b,Ndigits,N):
corps = RealField(Ndigits); R.<X> = corps[]
T = TCHEBYCHEV(R,N+1)
ZEROS = T.roots(corps)
zeros = [ZEROS[k][0] for k in range(N+1)]
aa = corps(a); bb = corps(b)
Lx = [aa + ((bb-aa)/2)*(zeros[k]+1) for k in range(N+1)]
Ly = [f.subs({t:Lx[k]}) for k in range(N+1)]
for j in range(N):
p = NEWTONAITKEN(R,Lx,Ly)
return p"
1 réponse
-
ContributeurSalut
p semble défini dans la boucle for uniquement.
Tu devrais formater ton code pour que qqun t'aider correctement.-
-
Contributeur@mameC'est le bouton "Code" dans l'éditeur de message du forum.
-