Mathematica, Matlab 및 Python과 같은 언어로 한 차원 (예 : Talbot 또는 Euler 알고리즘)으로 라플라스 변환을 수치 역으로 변환하는 알고리즘에 사용할 수있는 많은 구현이 있습니다. 여러 차원의 경우 또는 언어를 구현할 수있는 코드를 찾고 있습니다. 이 사실을 알고 계시다면 알려 주시기 바랍니다.
내가이 재발하고 있어요 "N + 1 = (3분의 13)이 N을 X X - (4/3)는 N-1을 X". x n의 주어진 값을 가진= 1 및 x = 1의 값을 가진 첫 번째 50 개의 값을 인쇄하는 파이썬 스크립트를 작성하려고합니다. import math
def printRecurrence():
x = [0]*51 #initialize list of
최근에 C++을 사용하기 시작했고 사용자 정의 된 ode 's 시스템의 통합을 허용하는 클래스를 만들었습니다. 성능을 비교하기 위해 두 개의 서로 다른 통합자를 사용합니다. 여기서, 코드의 일반 레이아웃은 : class integrators {
private:
double ti; // initial time
double *x
저는 이분법에서 몇 가지 코드를 수행하는 중입니다. 물론 스크립트를 작성하고 실행하는 것은 쉽습니다. 그러나 x^2와 같은 표현식을 입력으로 사용하는 함수를 정의하고자 할 때이 트릭이 제공됩니다. 나는 함께 무언가를 폐기했다. 일종의 일을 끝내지 만 첫 번째 if 문에서는 함수로 함수를 곱할 수 없다는 것에 대해 불평합니다. 어떻게하면이 문제를 해결할 수
나는 중력 포텐셜에 대한 식 (here에서 15)을 얻었고, 궤도를 계산하기 위해 국부적 인 그래디언트 인 중력을 평가할 필요가 있다는 것을 의미합니다. 단일 값으로 수만 번 P2, P4 및 P6의 Legendre polynomials 파생물을 평가합니다. 나는 this question의 표현을 사용하여 계산할 수 있지만, 명시 적으로 유한 차분으로 유도체
입니다. 최대 우도 추정량에 대해 연구 중이며 매개 변수 중 하나는 digamma 함수를 사용하여 추정됩니다. 나는 방정식을 풀기 위해 uniroot를 사용하려하지만 그렇게 할 수 없다.이 dig = function(alpha){
digamma(2 + alpha) - digamma(alpha) - (1/(2+alpha)) + (2/(2+alpha)
NumPy와 quad부터 scipy.integrate까지 다음과 같은 정수를 수치 적으로 풀려고합니다. 코드는 좀 작업,하지만 난 결과 데이터에 가짜 노치를 얻을 : 누구 그들이 발생하는 방법과 올바른 부드러운 결과를 얻을하는 이유 아이디어가있다? 여기 파이썬에서 원래의 코드입니다 : #!/usr/bin/env python
import numpy as
float 배열의 제곱의 합을 계산하려고합니다. 반올림 오류를 줄이는 방법은 무엇입니까? 실제 프로그램의 내부 루프에서 약 5,000,000 개의 수레를 합치려고합니다. Test.cpp에 :가 #include <iostream>
#include <stdint.h>
template <typename Sum, typename Element>
Sum sum
현재 C++의 유한 차분 방법을 사용하여 미분 방정식을 풀기위한 작은 프로그램을 작성하고 있습니다. 문제는 내가 tridiagonal 매트릭스를 해결하기 위해 lapack을 사용하고 있지만 루틴이 호출 된 후 Segmentation Fault가 발생합니다. 루틴은 여전히 info = 0로 종료되므로 문제가 무엇인지 알 수 없습니다. 여기에 프로그램이 있