Bonjour,
J'essaye de résoudre une EDP sur matlab. J'ai trouvé (par calcul) que U(x,t)=(exp(-t)+8/(pi)^3)*sin(pi*x)+8/(pi)^3*sin(x)*Somme(...)

J'ai calculé ma somme mais mon programme ne marche pas:
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
clear all
close all
clc
 
N=1000;
 
% Initialisation somme S:
S=1/27;
 
% Calcul somme:
for i=2:N
    S=S+1/(2*i+1)^3;
end
S
 
for t=0:0.1:1
    for x=0:0.1:1
        U=(exp(-t)+8/(pi)^3)*sin(pi*x)+8/(pi)^3*sin(x)*S;
    end 
end
 
 
 
% Affichage solution en 3D:
imagesc(x,U)
Est-ce que quelqu'un saurait d'où vient le problème ? Merci !