output.m

别人的边界元程序

function zx=output()

%% output the calculated data from these equations

global  tkt tqt
global nodbs 
global l  sol

%%
fidw=fopen('bem3.out','a');
fprintf(fidw,'******results of the boundary nodes****** ')
fprintf(fidw,'\n');
for nj=1:nodbs
    fprintf(fidw,'%4i %4.4e %4.4e \n',nj, tkt(nj),tqt(nj));
end

fprintf(fidw,'******results of the inner nodes****** ')
fprintf(fidw,'\n');
for nj=1:l
    fprintf(fidw,'%4i %4.4e',nj, sol(nj));
end

%%
zx=1;
fclose(fidw);
return