Recurrence Relation

Solution 

mydeterm.m 

function s = mydeterm(A)

n=length(A);

if n==1

s=A;

else s=0;

for j=1:n

C=A;

C(1,:)=[];

C(:,j)=[];

s=s+A(1,j)*(-1)^(j+1)*mydeterm(C);

end

end 

myfac.m 

function s = myfac(n)

if n==0

s=1;

else s=n*myfac(n-1);

end 

myfib.m 

function s = myfib(n)

if n==1

s=1;

elseif n==2

s=1;

else s=myfib(n-1)+myfib(n-2);

end

Solutions.m 

% Question 1.

myfib(25)

% Question 2.

myfac(80)

% Question 3.

A=magic(5);

mydeterm(A)

B=magic(6);

mydeterm(B)