Mesh, Octave and Freefem++

[fotios]

Once I wrote some simple interface to generate a mesh through Octave 
using Freefem++ [1], after reordering my code I ended up with the 
attached (non documented) functions and here is a simple example

b1 = boundary (1, "cos(t)", "sin(t)", [0, 2*pi], 50);
b2 = boundary (2, "0.5*cos(t)", "0.5*sin(t)", [0, 2*pi], -50);
s = fixmesh (5, 0.05);
[p, e, t, hmax] = buildmesh (s, b1, b2);
triplot (t, p(:,1), p(:,2))
axis equal;

If anybody is interested I could write documentation and make it a real 
package (maybe with some help)

/Fotis

[1] http://www.freefem.org/ff++/

boundary.m
Description: Text document

buildmesh.m
Description: Text document

fixmesh.m
Description: Text document

readmesh.m
Description: Text document