- 01
- 02
- 03
- 04
- 05
- 06
- 07
- 08
- 09
- 10
- 11
- 12
- 13
- 14
- 15
- 16
- 17
- 18
- 19
- 20
- 21
- 22
- 23
- 24
- 25
- 26
- 27
- 28
- 29
- 30
- 31
- 32
- 33
- 34
- 35
- 36
- 37
- 38
- 39
- 40
- 41
- 42
- 43
- 44
- 45
- 46
- 47
- 48
#include "math.h"
#include <iostream>
void Resolve(
float a, float b, float c, // коэффициенты уравнения
int& count, float x[2] // число корней, сами корни
// если корней нет, то значения x[0] и x[1] не определены
// если корень один, то значение x[1] не определено
// если корней бесконечно много, то вернуть INT_MAX
// предполагается, что a,b,c и корни по модулю не превосходят 1E10
)
{
const float MAXFREV = 1e-10f;
if (a==0.0 && b==0.0 && c==0.0)
{
count = INT_MAX;
} else
{
count = 0;
float d = b*b-4.0f*a*c;
if (d==0.0f)
{
if (2.0f*abs(a)>abs(b)*MAXFREV)
x[count++] = -b/(2.0f*a);
} else if (d>0.0f)
{
float num = -b - sqrt(d);
if (2.0f*abs(a)>abs(num)*MAXFREV)
x[count++] = num/(2.0f*a); // (-b-sqrt(d))/2a
if (abs(num)>2.0f*abs(c)*MAXFREV)
x[count++]=(2.0f*c)/num;
}
}
}
int main ()
{
float x[2];
int count;
Resolve (0.0001f,2.000f,0.001f, count, x);
std::cout << "count= " << count;
if (count<=2)
for (int i=0; i<count; ++i)
std::cout << ", x[" << i << "]= " << x[i];
std::cout<<std::endl;
return 0;
}
Просто решение квадратного уравнения.
На флоатах.
Числа подобраны так, что классическая формула лажает, выдавая второй корень -0.000596 вместо -0.0005