Séance 11 - Tir au but

Il s'agit bien d'une équation linéaire,non? vous obtenez quelle solution pour C? Pcq j'arrive jms au bon but...(mon C vaut 6.200 et mon y(end)=2.098 au lieu de -1).
Et, en fait, c'est quoi le z(end) ? (ou u2(end) selon les préférence) Il est pas sensé être égal à C, vu que c'est la dérivée de y, et donc la pente?
Merci bcp si qqn pouvait me dire quoi

j'ai ça comme résultat :

C =

-4.200169112285897


yfin =

-1.000000000000004


mon code ne manipule aucun u2(end) mais des u1(end). Le but est de rapproché ce u1(end) de la limite. (dans mon cas je juge qu'il est suffisamment proche de la valeur limite, si tu veux encore être plus proche augmente ton N)

je joins mon code :

Code:

%% TIR AU BUT (PROBLEME LINEAIRE)

clear;clc;clf;
format long
hold on
grid on
a=0;
b=pi;
N=1000;
h=(b-a)/N;
f=inline('x-1-(sin (x))*y');
x=a:h:b;
lim=-1;
C1=1;
C2=-1;

u1(1)=1;
u2(1)=C1;
for i=1:N
    u1(i+1)=u1(i)+h*u2(i);
    u2(i+1)=u2(i)+h*f(x(i),u1(i));
end
yfin1=u1(end);
plot(x,u1,'g');
u2(1)=C2;
for i=1:N
    u1(i+1)=u1(i)+h*u2(i);
    u2(i+1)=u2(i)+h*f(x(i),u1(i));
end
yfin2=u1(end);
plot(x,u1,'r');
C=C1+(C1-C2)*(yfin1-lim)/(yfin2-yfin1)
u2(1)=C;
for i=1:N
    u1(i+1)=u1(i)+h*u2(i);
    u2(i+1)=u2(i)+h*f(x(i),u1(i));
end
yfin=u1(end)

plot(x,u1,'k');
hold off
pause
plot(x,u1,'k');
grid on