accelerating code , scalar operations

Asked by Sara on 11 Jul 2012

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.

0 Comments

Sara

Products

No products are associated with this question.

0 Answers

Contact us