#include?
#include?
#include
int?lud（double*?m?int?n?int*?index）?{
????int?i=0?j=0?k=0?i_max?=?0;
????double?d_max=0?d_temp=0?d_dum=0;
????double?pivot[10];
????for（i=0;i<10;i++）{
????????pivot[i]=0;
????}

????for（i=0;?i????????for（d_max?=?0.0?j=0;?j????????????d_temp?=?fabs（m[i*n+j]）;
????????????if（?d_temp?>?d_max?）
		d_max?=?d_temp;
????????}
????????if（?d_max?==?0.0?）{
????????????return?0;
????????}
????????pivot[i]?=?1.0?/?d_max;
????}
???
????for（j=0;?j????????for（i=0;?i????????????for（d_temp?=?m[i*n+j]?k=0;?k		d_temp?-=?m[i*n+k]?*?m[k*n+j];
????????????}
????????????m[i*n+j]?=?d_temp;
????????}
????????for（d_max=0.0?i=j;?i????????????for（d_temp?=?m[i*n+j]