ntriangles=length(t(1,:));
npoints=length(p(1,:));
nrand=length(e(1,:));

t1=t;
t1(4,:)=[];

fout=fopen('FEM_daten','w');
fprintf(fout,'%8i  %8i  %8i\n',ntriangles,npoints,nrand);
fprintf(fout,'           \n');
fprintf(fout,'    %5i  %5i  %5i\n',t1);
fprintf(fout,'%7.4f   %7.4f\n',p);

