function y = symmetry(x)
n=length(x);
for i=1:n
if x(i)==x(n-i+1)
y=1;
else
y=0;
end
end
end
unction a = FillWithOnes(n,mat)
m=size(mat,1);
a = zeros(n);
for i=1:m
b=mat(i,:);
c=b(1,1);
d=b(1,2);
a(c,d)=1;
end
end
function m = timestables(x)
for i=1:x
m(i,:)=i:i:i*x;
end
end
function out = findIndices(vec, thresh)
out =find(vec>thresh);
end
function tf = mono_increase(x)
if isempty(find(diff(x)<=0))
tf = true;
else
tf = false;
end
function y = existsInVector(a,b)
if any(find(b(:)==a))==1
y=1;
else
y=0;
end
function B = swap_ends(A)
m=size(A,1);
n=size(A,2);
B=zeros(m,n);
for i=1:n
if i==1
B(:,i)=A(:,end);
elseif i==n
B(:,end)=A(:,1);
else
B(:,i)=A(:,i);
end
end
end
function [q,r] = swapInputs(a,b)
[q,r] = deal(b,a);
end
function B = column_removal(A,n)
A(:,[n]) = [];
B=A;
end
function c = hypotenuse(a,b)
c = sqrt(a^2+b^2);
end