%anonymous function
% function_name = @ (args) expr
%cylinder = @(radius, height) pi*radius^2 *height;

% stored functions
function [area, volume] = cylinder(radius, height)

% default values
if nargin ==1
    height = 1;
end

base = pi* radius^2;
area = 2*pi*radius*height+2*base;

if nargout==2
    volume = base * height;
end