求解方程组,编程中如何解决
我有种解方程组的思路,比如:
x+y=a
x*y=b
a,b是常数,求解这种类型的方程组,我这样写的:
double t1,t2,t3;
read(t1); /*输入步长,精度*/
read(t2); /*输入常数a*/
read(t3); /*输入常数b*/
for(x=0;x<100;x=x+t1)
for(y=0;y<100;y=y+t1)
if (x + y > t2 - t1 && x + y < t2 + t1 && x * y > t3 - t1 && x * y < t3 + t1 )
write x,y; /*输出x,y*/
但这样解极不精确,t1越精确的话,运行速度极慢,t1=1都有些慢,t1=0.1的话等了半天都没结果。
要怎么解决求解方程组的问题?