Problem with the shooting method

2 views (last 30 days)
Melda Harlova
Melda Harlova on 10 May 2019
Commented: Jan on 10 May 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 Comments
Torsten
Torsten on 10 May 2019
Replace "e^(-t)" by "exp(-t)".
Jan
Jan on 10 May 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.

Sign in to comment.

Answers (0)

Community Treasure Hunt

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

Start Hunting!