Someone please optimize this code or write is a little different, so that it doesn't take that much time to run. I really need it, I'd appreciate it alot!

close all
format long
caught = false;
w = 1.61803;
while ~caught
w = w + 0.0000010000000000;
N = 100000000;
x = zeros(N, 1);
y = zeros(N, 1);
t = 0;
dt = 1 / (N - 1);
for k = 2:N
xp = t - x(k-1);
yp = 1 - y(k-1);
C = w / sqrt(xp * xp + yp * yp);
xp = C * xp;
yp = C * yp;
x(k) = x(k-1) + dt * xp;
y(k) = y(k-1) + dt * yp;
if (y(k-1) < 1) && (y(k) >= 1)
caught = true;
w
y(k-1);
y(k);
k;
end;
t = t + dt;
end
end

Someone please optimize this code or write is a little different, so that it doesn't take that much time to run. I really need it, I'd appreciate it alot!

[Walter Roberson]

Your for k loop continues even after you have caught something. Is that what you want? I think you should probably be breaking out of the for k loop as soon as you hit caught the first time.

[Catalytic Kys]

How do I do this?

[Walter Roberson]

if (y(k-1) < 1) && (y(k) >= 1)
caught = true;
w
y(k-1);
y(k);
k;
break
end;
if caught; break; end

[Catalytic Kys]

clear;
close all;
clc;

format long
    caught = false;
    omega = 0;
    while ~caught
        omega = omega + 0.000001;
        T = 1;
        steps = 100;
        x = zeros(steps, 1);
        y = zeros(steps, 1);
        t = 0;
        dt = T / (steps - 1);
        for i = 2:steps
            xSecond = t - x(i - 1);
            ySecond = 1 - y(i - 1);
            C = omega / sqrt(xSecond * xSecond + ySecond * ySecond);
            xSecond = C * xSecond;
            ySecond = C * ySecond;
            x(i) = x(i-1) + dt * xSecond;
            y(i) = y(i-1) + dt * ySecond;
            if (y(i-1) < 1) && (y(i) >= 1)
                caught = true;
                omega
                y(i-1);
                y(i);
                    i;
            break;
            end;
            if caught; break; end;
        t = t + dt;
        end;
 end;

Does this script work for you?

[Walter Roberson]

omega =
   1.609883999957745

is output within a reasonable time.

[Catalytic Kys]

What happens if you change the step to 1000 instead of 100?

Thanks!