function x_new = psi_Extrap(t, x, tau, k, F)

x_block = zeros(length(F));

for row = 1:length(F)
    x_block(row,1) = multi_RK2(F(row),t,x,tau);
end

for col = 2:k+1
    for row = 1:length(F)+1-col
    r = (F(row+col-1)/F(row))^2;
    x_block(row,col) = 1/(r-1)*(r*x_block(row+1,col-1)-x_block(row,col-1));
    end
end

x_new = x_block(1,k+1);

end

Published with MATLAB® R2014a