normaliser une fonction numeriquement

[lyoko]

bonsoir à tous!!!


mon problème est le suivant: j'ai résolu l'eq de schrodinger numériquement (en C ) est j'ai trouvé les fct d'ondes que je voudrai normaliser en utilisant la formule du trapezoide!
mais voilà: les fonctions normaliser que j'obtiens sont toutes nulles! je crois pourtant que j'ai bien défini ma norme!

voici la partie de mon programme qui definie ma norme: ( les valeurs des fct d'onde sont contenues dans un tableau appelé phi_sol[])
/****************************** ********************/

alfa_sol=gamma*ksi_sol;
phi_sol[0]=exp(-alfa_sol*h);
phi_sol[1]=1.;

norme=h*phi_sol[0]*phi_sol[0]/2 +h*phi_sol[1]*phi_sol[1];

for(k=1 ; k<(n+7)*k_max; k++)
	{
	   	  
alfa_sol=gamma*ksi_sol;
t[0]=-h;
t[1]=0.;
U[0]=0.;
U[1]=1.;

F_sol[0]= gamma*gamma*(ksi_sol*ksi_sol - U[0] );
F_sol[1]= gamma*gamma*(ksi_sol*ksi_sol - U[1] );

phi_sol[0]=exp(-alfa_sol*h);
phi_sol[1]=1.;

		t[k]= -h+k*h;
		U[k]=u(t[k]);

		t[k+1]=-h+ (k+1)*h;
		U[k+1]=u(t[k+1]);


F_sol[k]= gamma*gamma*(ksi_sol*ksi_sol - U[k] );
F_sol[k+1]= gamma*gamma*(ksi_sol*ksi_sol - U[k+1] );


phi_sol[k+1]= ( (2 + (5*h*h/6) *F_sol[k])*phi_sol[k] - (1- (h*h/12)*F_sol[k-1])*phi_sol[k-1] )/( 1- (h*h/12)*F_sol[k+1]) ;
norme= norme+ h*phi_sol[k+1]*phi_sol[k+1];
	}
norme= norme + h* phi_sol[(n+7)*k_max+1]*phi_sol[(n+7)*k_max+1]/2;
norme= sqrt(norme);

/****************************** ************************/
quelqu'un peut m'aider!!! merci:
-----

[cherwam07]

Salut,

Ca n'est peut-etre pas ca du tout, mais tu es bien certain de travailler en flottant pour TOUTES les variables que tu utilises (à part k, l'indice de boucle )

[lyoko]

salut cherwam!!!

merci pour ta suggestion c'est très gentille de répondre mais ce n'est pas ça ! toutes mes variables je les ai défini comme double sauf les i et les k c'est des long!!!!

alors si quelqu'un d'autre a une suggestion je suis prenante!!! merci

[Jack]

Fais tourner ton programme en pas à pas en observant la valeur des variables au fur et à mesure et compare avec les résultats auxquels tu t'attendais. Lorsque les valeurs seront différentes, l'analyse de la ligne fautive devrait te renseigner sur le type de l'erreur.

A+

[A voir en vidéo sur Futura]

[lyoko]

hello jack!
merci pour la réponse! mais je ne crois qu'il y a une erreur dans l'autre partie de mon programme, car si je dessine mes fct sans les normaliser, ces dernières ne sont pas nulles! en plus j'ai résolu ce problème ( resol de l'eq de schrod) avec une autre méthode et j'ai normaliser aussi mes fct ( elles ne sont pas nulle cette fois ci !!!) et j'ai pris la norme de ma deuxième méthode pour normaliser mes résultats obtenus avec numerov , et ça marche !!! j'ai le résultat souhaité! donc le problème est dans la façon avec laquelle j'ai défini ma norme en utilisant la règle du trapézoïde!!!!

alors d'autres propositions, je suis dans le rouge????

merci.