특정 잠재력을 가진 은하 궤도의 플롯을 생성하려고합니다. 내 코드는 function Eulersystem_MNmodel()
parsec = 3.08*10^18;
r_1 = 8.5*1000.0*parsec; % This converts our value of r into cm.
z_1 = 1.0;
theta_1 = 0.0; %Initial Valu
저는 RK4를 사용하여 ODE 시스템을 해결하고 있습니다. 저는 k3_1이 -3.1445e + 24로 제한되어 있기 때문에 직선 플롯을 생성합니다. 나는 그것이 왜 모자란 지 이해하지 못한다. function RK4system_MNModel()
parsec = 3.08*10^18;
r_1 = 8.5*1000.0*parsec; % in cm
z_1
Neville의 메소드를 구현 중이며 오류가 발생합니다. 유사한 오류에 대한 다른 응답을 보았지만 아무 소용이 없습니다.이 def f(x):
return x**2
n = 2
x = 1.5
xa = [1.0,2.0,3.0]
func = [f(k) for k in xa ]
P = [([0] * (n+1)) for i in range(n+1)
이분법을 사용하여 주어진 도메인 내의 방정식의 모든 근원을 찾고 싶습니다. 내 논리는 매우 간단합니다. 주어진 도메인을 하위 도메인으로 나눕니다. 주어진 하위 도메인에 대해 하나의 루트 만 존재합니다. f (a) * f (b) <이 0이면 다음 간격으로 건너 뜁니다. 그렇지 않으면 루트 (이등분)를 찾습니다. 내 문제 그러나 어떻게 적절한 하위 도메인 크기
친구는 현재 MLE (Maximum Likelihood Estimation) 방법을 사용하여 시계열 모델 SARIMAX (계절 ARIMA 외생 적)의 매개 변수에 대한 과제를 수행하고 있습니다. 그에게 사용 된 데이터는 외생 변수로 인도양 다이폴 (IOD) 지수가있는 2000 ~ 2012 년의 월별 강우량에 관한 것입니다. 그것은 SARIMAX 모델 분석을
C++에서 부드러운 함수의 수치 그래디언트를 계산하려고합니다. 그리고 매개 변수 값은 0에서 매우 큰 수 (아마도 1e10에서 1e20까지) 테스트 벤치로 함수 f (x, y) = 10 * x^3 + y^3을 사용했지만, x 또는 y가 너무 큰 경우 정확한 그래디언트를 얻을 수 없다는 것을 알았습니다. #include <iostream>
#include
배경 : 저는 Matlab에서 유한 차이 스키마를 실행하고 있습니다. 두 변수의 함수를 구별하고 싶습니다 f(x,y). 대신 격자 점 F(:,:)에 매트릭스 표현이 있습니다. F이 1 차원 배열이면 차별화 연산자는 행렬 Delta(:,:)이고 차별화는 F_tag = Delta*F입니다. 그러나 F이 매트릭스이므로 더 복잡해집니다. 기본 접근법은 n^2 1
Matlab의 erf/erfc 함수의 정확성에 문제가 있습니다. 우리가 알고 있듯이, x이 무한대이면 erf(x)은 1과 같습니다. 그러나 Matlab에서 나는 erf(6)이 이미 1과 같고 6도 상당하지 않다는 사실에 놀랐습니다! erfc(x)은 이 0이 아닌 반면에 erfc(28)은 0입니다. 이 함수의 수치 적 성능을 향상시킬 방법이 있습니까? 즉,
Matlab에서 polyeig 명령을 사용하여 Matlab에서 2 차 주문의 다항식 고유 값 문제를 해결했습니다. 나는 시스템이 단일 0 고유치를 갖는다는 것을 알고있다. (이것은 0 계수 행렬의 형태에 기인한다. 각 대각선 요소는 벡터의 동일한 행에있는 요소의 합을 -1 배로하여 1 1 1 ... 1)는 0 고유치를 가짐). I는 I 얻을 최소 고유치 주
복잡한 정수를 계산할 필요가 없으며이 문제 때문에 Fortran 77로 작성된 관리자 프로그램에서 얻은 기존 프로그램을 가져 왔습니다. 그러나 몇 가지 문제가 있습니다. 주로 DATA Statement의 구문 오류와 관련이 있습니다. FUNCTION CAUSSA(F,A,B,EPS)
IMPLICIT DOUBLE PRECISION (A-H,O-Z)