Bonjour à toutes et à tous,

Voilà j'ai un projet à rendre dans 2 semaines exactement sur la résolution de l'équation de Schrödinger stationnaire à 1D.
Nous avons une méthode qui nous est donnée, sur laquelle nous devons nous baser pour programmer cela pour plusieurs types de potentiels (potentiel plat, Lennard-Jones etc ...).

Néanmoins, j'ai appris la programmation seulement depuis 2 mois et je ne sais mais alors pas du tout comment débuter mon projet ...

Je vous mets le début de mon code :

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
#include <cstdlib> 
#include <vector> 
#include "matprod.h" 
#include <fstream>
 
#define h 6.62606957*pow(10,-34) 		// Définition de la constante de Planck
#define Pi 3.1415926 				// Définition de la valeur de Pi
 
// -------------------------------------------------- //
// ------ Définition du potentiel carré simple ------ //
// -------------------------------------------------- //
 
double V (double x, double xmin, double xmax,  double V0)
{
	if(x<xmin || x>xmax)
		return 0;
	else
		return V0;
}
 
 
int main() {
 
 
cout<<"###########################################################################"<<endl;
cout<<"###########################################################################"<<endl;
cout<<"########################## Projet Informatique ############################"<<endl;
cout<<"############# Résolution de l'équation de Schroedinger à 1D ###############"<<endl;
cout<<"###########################################################################"<<endl;
cout<<"###########################################################################"<<endl;
 
return 0 ;
 
}

Et ci-joint les quelques documents qui m'ont été donné :

Le pdf : projet M1 ==> Mon sujet correspond au premier
Le pdf : Notice ==> Explication d'une méthode de résolution

projets-M1.pdf
Notice 2014.pdf

Je ne souhaite pas forcément que l'on me donne des codes déjà fait, car je veux coder plus ou moins moi-même afin de comprendre ce que je fait.
En espérant avoir un petit coup de pouce

Merci d'avance !


-----------------------------------------------------
Je viens d'ajouter j'espère de façon correcte au niveau de la syntaxe la fonction potentiel J(x) :

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
#include <iostream> 
#include <cmath> 
 
using namespace std ; 
 
#include <cstdlib> 
#include <vector> 
#include "matprod.h" 
#include <fstream>
 
#define h 6.62606957*pow(10,-34) 		// Définition de la constante de Planck
#define Pi 3.1415926 				// Définition de la valeur de Pi
 
// -------------------------------------------------- //
// ------ Définition du potentiel carré simple ------ //
// -------------------------------------------------- //
 
double V (double x, double xmin, double xmax,  double V0)
{
	if(x<xmin || x>xmax)
		return 0;
	else
		return V0;
}
 
double J (double V (double x))
{
	double J = ((2*m)/(pow(h,2)))*(E - V (double x)) ;
 
	return J ;
}
 
 
int main() {
 
 
cout<<"###########################################################################"<<endl;
cout<<"###########################################################################"<<endl;
cout<<"########################## Projet Informatique ############################"<<endl;
cout<<"############# Résolution de l'équation de Schroedinger à 1D ###############"<<endl;
cout<<"###########################################################################"<<endl;
cout<<"###########################################################################"<<endl;
 
cout << V << endl;
 
return 0 ;
 
}