#include <stdio.h>
#include <math.h>
int main(void) {
float a = 1.0f, b = 10.0f, c = 0.001f;
float r1, s1, r2, r3, r4;
r1 = (-b + sqrtf(b * b - 4 * a * c)) / 2 * a;
s1 = (-b - sqrtf(b * b - 4 * a * c)) / 2 * a;
r2 = c / (a * s1);
r3 = -c * 2 / (b + sqrtf(b * b - 4 * a * c));
r4 = -c / b * (1 + a * c / (b * b));
printf("r1 = %e, s1 = %e\nr2 = %e\nr3 = %e\nr4 = %e\n", r1
, s1
, r2
, r3
, r4
);
return 0;
}
I2luY2x1ZGUgPHN0ZGlvLmg+CiNpbmNsdWRlIDxtYXRoLmg+CiAKaW50IG1haW4odm9pZCkgewoJZmxvYXQgYSA9IDEuMGYsIGIgPSAxMC4wZiwgYyA9IDAuMDAxZjsKCWZsb2F0IHIxLCBzMSwgcjIsIHIzLCByNDsKCglyMSA9ICgtYiArIHNxcnRmKGIgKiBiIC0gNCAqIGEgKiBjKSkgLyAyICogYTsKCQoJczEgPSAoLWIgLSBzcXJ0ZihiICogYiAtIDQgKiBhICogYykpIC8gMiAqIGE7CgkKCXIyID0gYyAvIChhICogczEpOwogCglyMyA9IC1jICogMiAvIChiICsgc3FydGYoYiAqIGIgLSA0ICogYSAqIGMpKTsKCQoJcjQgPSAtYyAvIGIgKiAoMSArIGEgKiBjIC8gKGIgKiBiKSk7CiAKCXByaW50ZigicjEgPSAlZSwgczEgPSAlZVxucjIgPSAlZVxucjMgPSAlZVxucjQgPSAlZVxuIiwgcjEsIHMxLCByMiwgcjMsIHI0KTsKIAoJcmV0dXJuIDA7Cn0=