numerical-methods

-4熱

2答えて

ので、行2で、A [0] [0] =次いで0 がどのように行列「A」の1行を入れ替えることができガウス消去法における行列の最初の値を仮定します正しい結果が得られるか？

2熱

1答えて

私は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

-6熱

1答えて

float 'オブジェクトに属性' __getitem__ 'がありません

私は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 i

0熱

1答えて

方程式が変わる速度をどのように変えて、私の二分法の適切な間隔が出てくるかを調べるにはどうすればいいですか？

私は二分法を使って所与の領域内の方程式の根をすべて探したいと思います。私のロジックは非常に単純です。与えられたドメインをサブドメインに分割して、与えられたサブドメインに対して最大でも1つのルートしか存在しないようにします。 f（a）* f（b）<が0ならば、次の区間にスキップし、そうでなければ根（二分法）を見つける。私の問題しかし、どのように私は適切なサブドメインのサイズを決定するのですか？頻繁に

0熱

1答えて

ガウス - ルジャンドル直交を使用した二重積分

二重積分を評価するためにガウス - ルジャンドル直交を使用したいと思います。私は、重みと横を生成するために、表記脚（X1、X2、m）を使用し、次のコード m=5000; %generate weights and abscissas [wx,xx]=leg(-1,1,m); [wy,xy]=leg(-1,1,m); %define function [email protected]

2熱

6答えて

数値勾配を計算するための「標準的な」方法はありますか？

私は、C++で滑らかな関数の数値勾配を計算しようとしています。パラメータ値はゼロから非常に大きな値（おそらく1e10から1e20まで）で変化する可能性があります。私はテストベンチとして関数f（x、y）= 10 * x^3 + y^3を使用しましたが、 xまたはyが大きすぎると、正しい勾配を得ることができないことがわかりました。ここではgraidientを計算するために私のコードです： #inc

0熱

1答えて

Matlabの行列演算子

背景：私はMatlabで有限差分スキームを実行しています。私は2つの変数f(x,y)の関数を区別したいと思います。代わりに、私はグリッドポイントF(:,:)上の行列表現を持っています。 Fが1次元配列の場合、微分演算子は行列Delta(:,:)となり、微分はF_tag = Delta*Fとなります。ただし、Fが行列であるため、これはより複雑になります。基本的なアプローチは、 1次元ベクトルとして

0熱

1答えて

Matlabのpolyeigの精度の向上

Matlabのコマンドを使用して、Matlabの2次の多項式固有値問題を解いています。私はシステムが単一の0固有値を持っていることを知っています（これは、各対角要素が同じ行の要素の合計の-1倍なので、ベクトル（1 1 1 ... 1 ）は0の固有値を持つ）。私は得る最小固有値が順1E-4（0固有値していると想定される）であるpolyeigコマンドを使用して、第二の最低がである場合、システムのサイ

0熱

2答えて

Fortran 90のDATA文の構文エラー

複雑な積分を計算する必要がありません。この目的のために、私は上司からFortran 77で書かれた古いプログラムを入手しました。主にDATA文の構文エラーに関連しています。これは、実際の積分を計算する機能を備えたコードの一部です： FUNCTION CAUSSA(F,A,B,EPS) IMPLICIT DOUBLE PRECISION (A-H,O-Z) external f

2熱

1答えて

次の積分を計算するためのMatlabの組み込み関数はありますか？

私は、次の積分の値を見つける必要がある。これを計算するためのMATLABで組み込み関数がある場合誰もが知っていますか？