[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: octave-speed
From: |
Paul Thomas |
Subject: |
Re: octave-speed |
Date: |
Sun, 6 Feb 2005 08:35:18 +0100 |
Rafael,
Just in case I have confused things totally, following Paul K's last
message:
octave:14> system("cat version2.m")
tic;
Ns=1e5;
m=[1:(Ns-1)];
T=(sqrt(m.*(m+1)));
n=(8^(0.5^Ns))*cumprod(T.^(0.5.^(Ns-m)));
printf("time=%g\n",toc);
printf("n=%g\n",n(Ns-1));
ans = 0
octave:15> version2.m
time=0.179999
n=99998.5
Paul T
-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.
Octave's home on the web: http://www.octave.org
How to fund new projects: http://www.octave.org/funding.html
Subscription information: http://www.octave.org/archive.html
-------------------------------------------------------------