当前你的浏览器版本过低，网站已在兼容模式下运行，兼容模式仅提供最小功能支持，网站样式可能显示不正常。
请尽快升级浏览器以体验网站在线编辑、在线运行等功能。
3161：Numerical Integration
题目描述
In his freshman year, flymouse studied integral. Symbolic integration frustrated him a lot. He was often confused by those subtle techniques: alternately trying integration by substitution and integration by parts until the integrand appears in a tabulated form then employing integration by quadrature to have the work finally done. In contrast, numerical integration intrigued him much. By following some fixed procedures such as Newton-Cotes formulas, he could easily find the approximation to some definite integrals without too much effort. Given the exercises flymouse did, can your program do as nice a job as flymouse did?
输入解释
The input contains several test cases. Each test case consists of one line containing a univariate function f(x) and two real numbers a and b meaning the approximate value of 
is to be found. f(x) is expressed as a C-expression-like string. Allowed tokens are summarized in the following table:
| Token | Meaning | Example | Example in Mathematical Form |
+ | addition, positive | x+1, +x | x + 1, +x |
- | subtraction, negative | x-1, -x | x − 1, −x |
* | multiplication | x*x | x · x |
/ | division | x/2 | x⁄2 |
^ | power | x^-x | x−x |
( | left parenthesis | (x+1)*x | (x + 1) x |
) | right parenthesis | -(x+1) | −(x + 1) |
sin | sine | sin(x) | sin x |
cos | cosine | cos(x) | cos x |
tan | tangent | tan(x) | tan x |
log | natural logarithm | log(x) | ln x |
exp | exponential | exp(x) | ex |
asin | inverse sine | asin(x) | arcsin x |
acos | inverse cosine | acos(x) | arccos x |
atan | inverse tangent | atan(x) | arctan x |
abs | absolute value | abs(x) | |x| |
x | variable | x | x |
| constant | C floating-point constants without suffixes | 0, 1.5 | 0, 1.5 |
And the syntax is given below:
S | ::= | +S | -S | S+S | S-S | S*S | S/S | S^S | F(S) | (S) | x | constant |
F | ::= | sin | cos | tan | log | exp | asin | acos | atan |
Operator precedence and associativity are almost the same as those in C except that ‘^’ takes higher precedence than ‘*’ and ‘/’ and has right associativity.
f(x) is not necessarily continuous. It can have singularities which are always isolated if exist. It can be oscillatory, but it will never be oscillatory in the neighborhood of singularities where it blows up. The integral is guaranteed to converge.
Length of the interval over which f(x) is integrated is not longer than 10.
Process to end of file.
输出解释
For each test case, output one line with only the approximate value of the integral rounded to exactly four digits past the decimal point.
输入样例
x+1 0 1 +x 0 1 x-1 0 1 -x 0 1 x*x 0 1 x/2 0 1 x^-x 0 1 (x+1)*x 0 1 -(x+1) 0 1 sin(x) 0 1 cos(x) 0 1 tan(x) 0 1 log(x) 1 2 exp(x) 0 1 asin(x) 0 1 acos(x) 0 1 atan(x) 0 1 abs(x) 0 1 x 0 1 1 0 1 sin(x)/x 0 1
输出样例
1.5000 0.5000 -0.5000 -0.5000 0.3333 0.2500 1.2913 0.8333 -1.5000 0.4597 0.8415 0.6156 0.3863 1.7183 0.5708 1.0000 0.4388 0.5000 0.5000 1.0000 0.9461
提示
- Be cautious about outputting ‘
-0.0000’. - Make your integration module as robust as possible.
最后修改于 2020-10-29T06:54:35+00:00 由爬虫自动更新
共提交 0 次
通过率 --%
| 时间上限 | 内存上限 |
| 5000 | 131072 |