Files
2025-09-12 18:55:25 +09:00

11 lines
198 B
C

extern unsigned long ija[];
extern double sa[];
void asolve(unsigned long n, double b[], double x[], int itrnsp)
{
unsigned long i;
for(i=1;i<=n;i++) x[i]=(sa[i] != 0.0 ? b[i]/sa[i] : b[i]);
}