Méthode numérique

je sais qu'il est encore temps etque c'est vachement à l'avance, mais je me demandais:

qu'est ce que monsieur Bastin demande dans son oral? car son cours est assez ...pas grand chose en fait :/

Il demande de lui débaler tout un chapitre, choisi au hasard... Ni plus, ni moins...

Si tu lui fais quinze dessins pour monte-carlo en rajoutant chaque fois un point, je suis sur qu'il en rira. =p

et il pose des questions de compréhension de l'algorithme ou un truc du genre? parce que ses chapitres c'est plutot ca, y'a peut être 3-4 démos, pas bien longues en plus et le reste c'est l'algorithme :/

Petite question : apparemment, pour l'oral, on a droit à certains slides. Quelqun saurait m'indiquer lesquels?

EDIT : C'est bon, j'ai trouvé. Ils sont sur le FTP.

si par ailleurs, certains ont des pistes de résolution pour les examens de l'année passée ou celui de cette année des géométrologues et climatologues que nous avons reçu par mail ,ce serait sympa...parce que waw quand même

Petite note dans le mail de jérémie: les questions d'avant 2009 ne sont pas de lui

Il faut bien voir qu'il y a des questions d'interpolation polynomiale dans les années précédentes il me semble... et ça c'est pas vu. Il faut tenir compte du fait que l'assistant a changé et si je ne m'abuse ils sont rares ceux qui ont raté. Donc ne nous énervons pas avant l'heure. Faut se concentrer sur optique.

Ils ont toujours été rares, ceux qui rataient l' "écrit", je pense. Sans mettre de pression...

Lindsay, mes souvenirs sont trop vague pour pouvoir en dire plus que ce que j'ai dit. Ca fait deux ans déjà... =/

l'interpolation polynomiale, je pense que c'était dans ses questions de l'an passé justement!

Maintenant il faut voir comment jérémie cote aussi, la je viens d'en essayer 2 de cette année que les géologues ont eu l'an passé, certains trucs ca va, d'autre limite >-< on verra, tu as raison gaêl, d'abord optique ^^

oui mais l'interpolation polynomiale on l'a pas vue donc pas de stress pour l'histoire de l'interpolation avec un polynôme d'ordre 12

sinon, tu tapes "numerical recipes" dans google, y'a un lien d'un livre en ligne qui te propose tous les algorithmes des tps! grmlgrmlgmrl

C'est déjà pour combien de % que les TP de l'année comptent?

L'évaluation est pondérée comme suit:

- 50% pour l'examen oral (partie théorique du cours)
- 35% pour l'examen pratique (un programme à réaliser)
- 15% pour les travaux rendus pendant l'année

Si la cote de l'examen pratique est en-dessous de 10, elle fait office de cote finale.

A votre service ! ^^

ok merci

quelqu'un pourrait m'aider? j'ai regardé ce programme à faire

Une matrice carree A(x) de dimension N X N et dependant du parametre x est definie
par :
Aij(x) = 10 x^((j-1)/i)
avec 1 <= i; j <= N.
Realisez un graphique de la valeur du determinant de A(x) pour x compris entre 0 et 1
pour les dimensions N = 3 puis N = 4. Utilisez pour ce faire la fonction det() que vous
avez dedefinie a l'aide de l'algorithme de la decomposition LU, importez donc votre chier
LUdecomp.cpp.
Pour le graphique, ecrivez dans un fichier determinants.txt 200 valeurs de paires fx; det(A(x))g
pour x compris entre 0 et 1.


et je fais ca pour remplir ma matrice et il me dit "stack around the variable a was corrupted"
sur internet ils me disent que je vais trop loin dans le tableau mais je vois pas où ni comment faire :s
et lerreur ne peut être que là car le reste c'est un tp déjà fait

bool sing;
double d;
double x=0.;
int const N=4;
double a[N*N];

for (int i=1;i<=N;i++){
for(int j=1;j<=N;j++) {
a[i*N+j]=10*(pow(x,((j-1)/i))) ;
}}

d = det(a, N, sing);

ofstream det1;
det1.open("fichier test1.txt");

for(double x=0.; x<=1.; x+=(1./200.)){
det1<< x << '\t' << d << '\n';
}
det1.close();
}

Ton tableau comporte 4*4 slots. Or, dans ta boucle, quand tu appelles l'élément pour j=i=N, tu remplis la case a[4*4] qui n'existe pas vu que, pour un tableau de 16 slots, le dernier élément est a[15]

le problème c'est que quand je me met i<N il me met des valeurs de déterminant avec des exposants +164 xD j'en déduit que ça pue

Ton problème est pas là. Tu effectues le bon nombre d'opération dans ta boucle mais tu remplis mal ton tableau : dis toi que le premier élément doit être le a[0] et le dernier a[15] et tu te rendras vite compte de ton erreur.

je peux pas commencer à 0 pour i ou j sinon dans mon exposant je divise par zéro >_< et je viens d'afficher la matrice ainsi formée...elle est singulière :/

Et comme ça?

for (int i=1;i<=N;i++){
for(int j=1;j<=N;j++) {
a[(i-1)*N+(j-1)]=10*(pow(x,((j-1)/i))) ;
}}

Ca marche?

c'est bon j'ai trouvé, je ne faisait pas varier mes indices partout de la même façon! soit une faute débile, et en exam, grâce à l'effet examen, je l'aurais jamais vue!

maintenant j'vais m'attaquer à celui sur les équa diff, il pue je le sens :p

En fin de compte, tu es bien partie comme je te l'ai expliqué?

non pas vraiment, j'ai gardé mes indices mais je fais la même chose en fin de compte merci!

De rien et bon courage alors

quand quelqu'un aura le temps (là le temps est plutôt à optique ^^)

est ce qu'il pourra m'expliquer comment traduire ça:

d²x(t)/dt²=-(GM/r(t)³)*x(t)
d²y(t)/dt²=-(GM/r(t)³)*y(t)
avec r= V(x²+y²)
.
Resolvez a l'aide de l'algorithme de Runge-Kutta ces deux equations, en posant toujours
G = M = 1;
x(0) = y(0) = 10;
sur un intervalle de temps de 400 secondes avec une precision de 0.01.
Examinez les trois situations dierentes :
vx(0) = 0.188, vy(0) = 0.188
vx(0) = 0.188, vy(0) = 0
vx(0) = 0.3, vy(0) = 0.2

en code avec le switch sur i pour les équa diff car je comprend pas du tout comment faire et ca me soule franchement,...je sens que 10/20 à cet exam c'est impossible pour moi, la moule de la progra ^^