function [ u, v, w ] = plot_channel_velocity( width, u0 )

M = 15;
yvec = -( width / M ) * (0:M);
zvec = -( width / M ) * (0:M);
[ y, z ] = meshgrid( yvec, zvec );
theta = 0;

[ u, v, w ]= channel_velocity( y, z, width, theta, u0 );

surf( y,z,u )
xlabel( 'y' )
ylabel( 'z' )
zlabel( 'u' )
shading interp