numcalu.hpp

#pragma once
#include "scmath.h"

class numcalu
{
public:
    static double solve(double a, double b, BasicNode* f, int n=1000)
    {
        double c;
        for (int i = 0; i < n; i++)
        {
            c = (a + b) / 2;
            if (scmath::caluExp(f,"x",c) == 0)
                return c;
            else if (scmath::caluExp(f,"x",c)*scmath::caluExp(f,"x",a) > 0) //选择异号的一边
                a = c;
            else
                b = c;
        }
        return c;
    }

    static double integ(double lower, double upper, BasicNode* f, double accuracy = 0.1)
    {
        double result = 0;

        for (; lower <= upper; lower += accuracy)
        {
            result += scmath::caluExp(f,"x",lower)*accuracy;
        }

        return result;
    }
};