Error ode45 / initial condition vector
bonjour à tous
je suis sur une boucle qui contient deux équa diff assez simples à résoudre, et à chaque fois la condition initiale de l'une est la valeure finale de l'autre
Jusque là tout va bien sauf qu'avant même de faire le premier tour il me sort
Citation:
??? Error using ==> ode45
Solving eq2 requires an initial condition vector of length 0
je vois pas ce que je peux faire pour sortir de là (un prob avec les global? un clear all ne fait pas avancer le schmilblik...)
le code de la boucle:
Code:
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
| t=1;
while t<10
[Lcpp,P]=ode45('eq1',[0 Lcp],[Pv]);
figure(1)
plot(Lcpp,P)
Pa=P(end);
[Lcmm,P]=ode45('eq2',[0 Lcm],[Pa]);
figure(2)
plot(Lcmm,P)
Pv=P(end);
end |
Auriez vous des idées pour me sortir d'affaire? :)
