MATLAB: How to write a code for an iteration

codecomputeiterationiterativescript

Hi there I would like to run an iterative code.

I would like to start with a matrix f0=[a,b,c,d]

Then I would like to compute the following operation

f1=[a+b,c+d,a-b,c-d]

I would like this to be done iteratively and so that the script uses the value before.

that is f100 uses the a,b,c and value from f99.

How would I do this?

Best Answer

f=zeros(100,4);
f(1,:)=[1 5 10 15];
for i=2:100
    f(i,:)=[f(i-1,1)+f(i-1,2),f(i-1,3)+f(i-1,4),f(i-1,1)-f(i-1,2),f(i-1,3)-f(i-1,4)];
end

Related Solutions

MATLAB: How to iterate the Newton method so that I get an array of values instead of just the last number

function R = yourNewton(f, df, x0, tol)
    R = [];
    ff = f(x0);
    while abs(ff) > tol
        B = x0 - ff/df(x0);
        ff = f(B);
        x0=B;
        R=[R;B,ff];
    end
end

MATLAB: Index exceeds the number of array elements at y2 for the finite difference formulas

In this line

y2(i) = (f(i+2) - 2*f(i+1) + f(i))/(deltaX.^2);

index of f is i+2, therefore, i should not go above n. Change the range of the for loop like this

for i = 1:n-2 % change this range
    %forward difference - first and second derivative
    y(i) = (f(i+1)- f(i))/deltaX;
    y2(i) = (f(i+2) - 2*f(i+1) + f(i))/(deltaX.^2); 
end

Similarly in the line

y4(i) = (f(i) -- 2*f(i-1) + f(i-2))/(deltaX.^2);

index of f is i-2, so i should not go below 3

for i = 3:n
    %backward difference - first and second derivative 
    y3(i) = (f(i)- f(i-1))/deltaX;
    y4(i) = (f(i) -- 2*f(i-1) + f(i-2))/(deltaX.^2);
end

Best Answer

Related Solutions

MATLAB: How to iterate the Newton method so that I get an array of values instead of just the last number

MATLAB: Index exceeds the number of array elements at y2 for the finite difference formulas

Related Question