Problem with the shooting method

1 visualización (últimos 30 días)
Melda Harlova
Melda Harlova el 10 de Mayo de 2019
Comentada: Jan el 10 de Mayo de 2019
Helloo,
Can someone help me with that code in matlab:
function main
y0=0;
y0=fsolve(@fun1,y0);
x0=[0.6,y0];
[t,x]=ode45(@fun2,[0 4], x0);
plot(t,x(:,1),t,x(:,2),'r')
function f=fun1(y)
x0=[0.6;y];
[t,x]=ode45(@fun2,[0 4],x0);
n=length(t);
f=x(n,1)+0.1;
function dx=fun2(t,x)
dx(1,1)=x(2);
dx(2,1)= -2*x(2)-2*x(1) + e^(-t) + sin(2*t);
That code must solves this task
x'' = -2x'(t) – 2x(t) + e^(-t) + sin(2t)
t ϵ [0,4]; x(0) = 0,6; x(4) = -0,1 with shooting method
but code doesnt run in matlab because of errors.
Thanks in advance.
  2 comentarios
Torsten
Torsten el 10 de Mayo de 2019
Replace "e^(-t)" by "exp(-t)".
Jan
Jan el 10 de Mayo de 2019
@Melda: A general rule in the forum: If you mention, that there is an error, post a copy of the complete error message. It is a good idea to share this important information with the readers.

Iniciar sesión para comentar.

Respuestas (0)

Categorías

Más información sobre Ordinary Differential Equations en Help Center y File Exchange.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by