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