x=linspace(0,4,400);    %Dividimos en 400 partes iguales el dominio de definición de x (perfil longitudinal) 
y=linspace(0,1,100);    %Dividimos en 100 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
Z=3-X;                  %Expresamos la función de la presión
h=pcolor(X,Y,Z);        %Visualizamos el campo de velocidades del fluido mediante una gama de colores
axis([0,4,-1,2])        %Restringimos el gráfico mostrado al intervalo que se nos solicita
print('presiones.png')  %Obtenemos una imagen del gráfico mostrado en formato .png


