解方程头文件

#ifndef EQUATION_H
#define EQUATION_H
#include<bits/stdc++.h>
using namespace std;
double one_equa(double x,double y){
	return -y/x;
}
double two_equax(double a,double b,double c,double d,double e,double f){
	double m=-d/a;
	double y=(f+m*c)/(e+m*b);
	double x=(c-b*y)/a;
	return x;
}
double two_equay(double a,double b,double c,double d,double e,double f){
	double m=-d/a;
	double y=(f+m*c)/(e+m*b);
	return y;
}
double one_twosolu_equax1(double a,double b,double c){
	double wa=-2100000.00;
	double delta=b*b-4*a*c;
    if(delta>0) {
        double x1=(-b+sqrt(delta))/(2*a);
        double x2=(-b-sqrt(delta))/(2*a);
        if(x1>x2) swap(x1,x2);
        return x1;
    }
	else if(delta==0){
        double x1=(-b+sqrt(delta))/(2*a);
        double x2=(-b-sqrt(delta))/(2*a);
        return x1;
    }
	else return wa;
}
double one_twosolu_equax2(double a,double b,double c){
	double wa=-2100000.00;
	double delta=b*b-4*a*c;
    if(delta>0) {
        double x1=(-b+sqrt(delta))/(2*a);
        double x2=(-b-sqrt(delta))/(2*a);
        if(x1>x2) swap(x1,x2);
        return x2;
    }
	else if(delta==0){
        double x1=(-b+sqrt(delta))/(2*a);
        double x2=(-b-sqrt(delta))/(2*a);
        return x2;
    }
	else return wa;
}
double one_inquar(double a,double b,double c){
	double wa=-2100000.00;
	double deeta=pow(b,2)-4*a*c;
	if(deeta<0) {
		return wa;
	}
	double x1=(-b+sqrt(deeta))/(2*a);
	double x2=(-b-sqrt(deeta))/(2*a);
	if(deeta==0) {
		return x1;
	}
	if(x2>x1) swap(x1,x2);
	return x1;
}
double one_inqual(double a,double b,double c){
	double wa=2100000.00;
	double deeta=pow(b,2)-4*a*c;
	if(deeta<0) {
		return wa;
	}
	double x1=(-b+sqrt(deeta))/(2*a);
	double x2=(-b-sqrt(deeta))/(2*a);
	if(deeta==0) {
		return x2;
	}
	if(x2>x1) swap(x1,x2);
	return x2;
}
#endif