我想求兩個二元一次方程的解,這是我的代碼
方程是ax+by=c,dx+ey=f
#include
#include
#include
int main()
{
double a, b, c, d, e, f;
scanf("%lf %lf %lf %lf %lf %lf", &a, &b, &c, &d, &e, &f);
assert(fabs(a*e-d*b) < 1.0e-8);
printf("x=%lf\ny=%lf\n", (c*e-b*f)/(a*c-d*b), (d*c-a*f)/(d*b-a*c));
return 0;
}
我想當解不唯一時程序異常退出,唯一時算出結果,所以故意輸入 5,1,8,20,4,6 按理應該異常退出啊,但是為什麼還是求出結果了呢?
兄弟,你的程序有問題。
#include <stdio.h>
#include <math.h>
#include <assert.h>
int main()
{
double a, b, c, d, e, f;
scanf("%lf %lf %lf %lf %lf %lf", &a, &b, &c, &d, &e, &f);
assert(fabs(a*e-d*b) < 1.0e-8); // 這裡應該是assert(fabs(a*e-d*b)>1.0e-8)
//此處還可考慮添加 assert(fabs(a*c-d*b))
printf("x=%lf\ny=%lf\n", (c*e-b*f)/(a*c-d*b), (d*c-a*f)/(d*b-a*c));
//此處應該是printf("x=%lf\ny=%lf\n", (c*e-b*f)/(a*e-d*b), (d*c-a*f)/(d*b-a*e));
return 0;
}
void assert(expression)
assert的作用是現計算表達式 expression ,如果其值為假(即為0),那麼它先向stderr打印一條出錯信息,然後通過調用 abort 來終止程序運行。
你輸入5,1,8,20,4,6 ,此時
fabs(d*b-a*e)<1.0e-8
為true,故不會異常退出。