Example:
x = [-10:0.1:10]';
Y = [x x.^2 x.^3 cos(x)];
location = 'top';
labels = ['x-axis '
'y=x '
'y=x^2 '
'y=x^3 '
'y=cos(x)'];
[hlines,haxes]=plots(x,Y,location,labels);

Plots the columns of the matrix Y versus the vector column x, using different axes at the top of the figure, with the specified axes-labels.

The location of the axes can be 'left' (default), 'right', 'top' or 'bottom'.

PLOTSES.M its the almost the same but with the axes moved from the fixed axis (useful?).

Both M-files come with the colored example of the screenshot.