disp('What is the answer to ... everything?')
a = input('> ');
 
disp(['You said ' num2str(a) '.']);
 
if a == 42,
  disp('That is correct.');
else
  disp('That is not correct.');
end;

function area = traparea(a, b, h)
% traparea(a, b, h)  Computes the area of a trapezoid given
%                    the dimensions a, b and h, where a and b
%                    are the lengths of the parallel sides and
%                    h is the distance between these sides
 
%  Compute the area, but suppress printing of the result
area = 0.5 * (a + b) * h;
end

function main
    area1 = traparea(1,2,3)
    disp(area1)
end
 
function area = traparea(a, b, h)
% traparea(a, b, h)  Computes the area of a trapezoid given
%                    the dimensions a, b and h, where a and b
%                    are the lengths of the parallel sides and
%                    h is the distance between these sides
 
%  Compute the area, but suppress printing of the result
area = 0.5 * (a + b) * h;
end

function main
  [r, th] = cart2plr(30, 40);
  disp(r);
  disp(th);
end
 
function [r, theta] = cart2plr(x, y)
%   cart2plr  Convert Cartesian coordinates to polar coordinates
%
%   [r, theta] = cart2plr(x, y) computes r and theta with
%
%       r = sqrt(x^2 + y^2);
%       theta = atan2(y,x);
 
r = sqrt(x^2 + y^2);
theta = atan2(y,x);
end

for i = 1:10
    for j = 1:10
        if j <= i,
            val = 0;
        else
            val = 1;
        end
        my_matrix(i, j) = val;
    end;
end;
 
disp(my_matrix);