# Output for GMT plotting

function [ n ] = gmtout( x, y, name )
  	n = length(x);

	fout = fopen(name,"w");

	for i=1:n
		fprintf(fout,"%f %f\n",x(i),y(i));
	endfor

	fclose(fout);
end