当我想通过RCaller运行一些代码时,Java会抛出以下异常:
Exception in thread "JavaFX Application Thread" com.github.rcaller.exception.ExecutionException: Can not send the source code to R file due to: java.io.IOException: The pipe is being closed Maximum number of retries exceeded.
这是我的密码:
protected void initialize
我有一个问题,因为这个工作适用于许多函数,但当我试图绘制正弦的积分时(我正在使用matlab 2010),我遇到了麻烦:
clear all
close all
clc
x = linspace(-10, 10, 100);
f = @(x) sin(x);
I = arrayfun(@(x) quad(f, 0, x), x);
plot(x, f(x),'r', x, I, 'b')
我希望有一个-cos(x),但是我得到了偏移量为1的东西,为什么会发生这种情况?该如何解决这个问题呢?
我需要编写一个python代码来计算1/(1+x^2)的积分(-5,5)的确切值。我知道答案是2弧坦(5),大致相当于2.746801.
下面是我所写的代码,但是我得到了一个稍微不同的答案,我想知道我是否可以做些什么来使这段代码更准确?谢谢你的帮助!
## The function to be integrated
def func(x):
return 1/(1 + x**2)
## Defining variables
a = -5.0
b = 5.0
dx = 1.0
Area = 0
## Number of trapezoids
n = int((b-a)/dx)
我试图使用python/scipy计算两个极限之间的积分。
我正在使用在线计算器双重检查我的结果(,),当我设定了某些限制时,我的结果不同意。
所用的代码是:
import scipy as sp
import numpy as np
def integrand(x):
return np.exp(-0.5*x**2)
def int_test(a,b):
# a and b are the lower and upper bounds of the integration
return sp.integrate.quad(integrand,a,b)
当设定(a,
我想计算一下功率。在一致性研究中,我采用了Walter公式来计算样本。然而,当我有主语的数量时,我需要用它来计算功率。
calcul_f <- function(alpha,puissance,ccc0,ccc1,nb){
cat("ccc0 correpond to the concordance to compare et ccc1 the new concordance et nb the number of replicate")
U_alpha <- abs(qnorm(alpha))
U_beta <- qnorm(puissance
我想在下面的代码中找到输出功率Po的积分: Vo = 54.6
# defining a function for duty cycle, output current and output power
def duty_cycle(output_voltage, array_voltage):
duty_cycle = np.divide(output_voltage, array_voltage)
return duty_cycle
def output_current(array_current, duty_cycle):
output_current =
我正在绘制数据集并进行对数拟合,但由于某些原因,拟合似乎是,在某个时候,我得到了足够好的拟合,但后来我绘制了该数据集,结果出现了这种不合适的情况。一开始有一个0.0 0.0076,但是为了避免渐近线,我把它改成了0.001 0.0076。
我正在使用(上面的图像不是这个,但现在我正在用这个测试,而且也有这个坏的匹配)这是适合的
f(x) = a*log(k*x + b)
fit = fit f(x) 'R_B/R_B.txt' via a, k, b
输出为
而且,有时它说7次迭代和上面的屏幕截图一样,其他的只有1次,当它做了“正确”的匹配时,它就像35次迭代或者什么的,如果我
我正在寻找一个java库/实现,它支持计算β分布的逆累积分布函数(也称为分位数的估计),并具有合理的精度。
当然,我已经尝试过,但是在第3版中,似乎仍然有一些。下面对导致这个问题的问题作了广泛的描述。
假设我想用大量的试验来计算贝塔分布的可信区间。在阿帕奇公域的数学里。
final int trials = 161750;
final int successes = 10007;
final double alpha = 0.05d;
// the supplied precision is the default precision according to the source cod
我将查找$I_m=\int_0^{2\pi} \prod_{k=1}^m cos(kx){}dx$,其中$m=1,2,3\ldots$
简单的SAGE代码:
x=var('x')
f = lambda m,x : prod([cos(k*x) for k in range(1,m+1)])
for m in range(1,15+1):
print m, numerical_integral(f(m,x), 0, 2*pi)[0],integrate(f(m,x),x,0,2*pi).n()
输出:
1 -1.47676658757e-16 0.000000000000
我试图运行一个SymPy集成,在这里我创建了自己的变量,积分是
import numpy as np
import seaborn
import math
import pandas as pd
import matplotlib as plt
import scipy.integrate as integrate
import scipy.special as special
from sympy import init_session
init_session()
%matplotlib inline
dP, dT = symbols('dP dT', real = T
我试图构建一些简单的运算,比如一个导数和积分函数来操作lambda函数,因为渐近和was都在努力地集成一些我传递给它们的东西。
导数函数没有给出任何问题,而是在绘制时返回输入函数的导数,但积分函数不返回,也没有绘制输入的正确积分。
import matplotlib.pyplot as plt
import numpy as np
from phys_func import func
sr = [-10,10]
x = np.linspace(sr[0],sr[1], 100)
F = lambda x: x**2
f = func.I(F,x)
plt.plot(x,F(x), lab
我试图计算一些积分,例如:
a = 1/sqrt(2);
b = -5;
c = 62;
d = 1;
f = exp(-x^2-y^2)*(erfc((sym(a) + 1/(x^2+y^2)*(sym(b)*x+sym(d)*y))*sqrt((x^2+y^2)*sym(10.^(c/10))))...
+ erfc((sym(a) - 1/(x^2+y^2)*(sym(b)*x+sym(d)*y))*sqrt((x^2+y^2)*sym(10.^(c/10)))));
h = int(int(f,x,-Inf,Inf),y,-Inf,Inf);
我已经创建了一个python脚本,它可以找到christoffel符号并求解测地方程。今天我想,在弯曲的2D空间中添加计算弧长的功能会很有用。我已经创建了使用渐近的代码,一切似乎都是正确的,但由于某种原因,它只是没有计算积分!它只是简单地显示了积分。谁能帮帮我。我已经尝试了所有的方法,甚至使用dsolve,但都不起作用。这是我的代码的当前状态: import sympy as sp
from sympy.interactive import printing
printing.init_printing(use_latex=True)
g = [[[11], [12]], [[21],