accelerating code , scalar operations
Info
This question is closed. Reopen it to edit or answer.
Show older comments
Hi Guys,
if k <= 0
int = x ;
return
end
int = zeros(1,length(x));
c=(k-1)/factorial(k-1);
dth2=dt/2;
lenx=length(x);
ft1=(lenx:-1:0)*dt;
for n = 3:lenx
y1 = x(2:n) ;
t1 = ft1(lenx-n+3:lenx+1) ;
y2 = x(1:n-1) ;
t2 = t1 +dt;
int(n) = sum (t1.^c.*y1 + (t2.^c).*y2)*dth2 ;
end
Is there any way to do .^c in the line ft1=(length:-1:0)*dt which is before and out of the for loop and get the same result as above code? Any suggestion would be appreciated.
Answers (0)
This question is closed.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!