MATLAB Series Solutions with Loops & Arrays
Uploaded by
Hiba Al-QaisyMATLAB Series Solutions with Loops & Arrays
Uploaded by
Hiba Al-QaisyComputer Applications 2024 Zahraa Al-Milaji
Approaches to Solving Series in MATLAB
Loops: Iterating through each term of the series using a loop structure.
Arrays: Utilizing MATLAB's ability to perform operations on entire arrays of data simultaneously, often
referred to as vectorization.
Using Loops in MATLAB
Syntax of for loops: for index = start: increment: end
Example: Summing the first n terms of a series using a for loop
s=0;
n=5;
for i = 1:n
s = s + i;
end
disp(['Summation of numbers 1:5 is: ',num2str(s)]);
n=5
sum_of_squares = 0;
for i = 1:n
sum_of_squares = sum_of_squares + i^2;
end
disp(['Summation of squaring numbers 1:5 is: ',num2str(sum_of_squares)]);
% Using arrays
numbers = 1:n;
squares = numbers.^2;
sum_of_squares = sum(squares);
disp(['Summation of squaring numbers 1:5 is: ',num2str(sum_of_squares)]);
Computer Applications 2024 Zahraa Al-Milaji
Example 2:
n=5;
sum_of_reciprocal_squares = 0;
for i = 1:n
sum_of_reciprocal_squares = sum_of_reciprocal_squares + 1/(i^2);
end
disp(['Summation of squaring numbers 1:5 is: ',num2str(sum_of_reciprocal_squares)]);
%using arrays:
numbers = 1:n;
reciprocal_squares = 1./(numbers.^2);
sum_of_reciprocal_squares = sum(reciprocal_squares);
disp(['Summation of squaring numbers 1:5 is:',num2str(sum_of_reciprocal_squares)]);
Example 3:
n=5;
sum_of_reciprocal_factorials = 0;
for i = 0:n
f=1;
for j= 1:i
f=f*j;
end
sum_of_reciprocal_factorials = sum_of_reciprocal_factorials + 1/f;
end
disp(['Summation of squaring numbers 1:5 is:',num2str(sum_of_reciprocal_factorials)]);
Example 4:
n=5;
sum_of_series = 0;
for i = 1:n
sum_of_series = sum_of_series + i^(i+1);
end
sum_of_series
Computer Applications 2024 Zahraa Al-Milaji
% Using arrays:
indices = 1:n;
terms = indices.^(indices + 1);
sum_of_series = sum(terms)
Example 5: Find elements of factor of 5 in 5x5 arrays:
% Create a 5x5 array
a = magic(5)
% Find the factors of 5 using logical indexing
f = a(rem(a, 5) == 0);
% Display the factors of 5
disp('Factors of 5 in the array:');
disp(f);