/* phi_tri_pts.h  generate points uniformly in a triangle       */
/*                input, np, is number of points+1 along a side */
/*                triangel x1,y1 x2,y2 x3,y3                    */
/*                return is the number of points in xs,ys       */

int phi_tri_pts(int np, double x1, double y1, double x2, double y2,
	        double x3, double y3, double xs[], double ys[]);