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