function cosine() x=%pi/3; y=cos(x); disp(y) endfunction