x=linspace(0,4,8);       %Dividimos en 8 partes iguales el dominio de definición de x (perfil longitudinal) 
y=linspace(0,1,10);      %Dividimos en 10 partes iguales el dominio de definición de y (perfil transversal) 
[X,Y]=meshgrid(x,y);     %Creamos la matriz que tiene por filas las coordenadas x y por columnas las coordenadas y
U=(Y.*(1-Y))/2;          %Determinamos la componente horizontal (i) de los vectores del campo
V=0*U;                   %Determinamos la componente vertical (j) de los vectores del campo
quiver(X,Y,U,V)          %Visualizamos el campo de velocidades del fluido
axis([0,4,-1,2])         %Restringimos el gráfico mostrado al intervalo que se nos solicita
print('velocidades.png') %Obtenemos una imagen del gráfico mostrado en formato .png