!	Compare numerical velocity with analytical for Hele-Shaw flow
clear all
clf reset
array=load('u.dat');
n=10;
y=linspace(0,1,n);
u=-4.*(y.^2-y);
plot(y,u,y,array(:,1),y,array(:,5))