function write(a,fid,buslist)% write polynomial loadsif ~a.n, return, end% filter loads using bus listidx = [];for i = 1:a.n  jdx = find(buslist == a.bus(i)*a.u(i));  if ~isempty(jdx), idx = [idx; i]; endendif isempty(idx), return, end% write Pl datafprintf(fid,'Pl.con = [ Pl.con; ...\n');fprintf(fid,['   ',a.format,';\n'],a.con(idx,:)');fprintf(fid,'   ];\n\n');