Hamilton-Jacobi 源代码

//$$jacobi.cxx jacobi eigenvalue analysis // Copyright (C) 1991,2: R B Davies #define WANT_MATH …

@cindex Jacobi elliptic functions @cindex elliptic functions (Jacobi) The Jacobian Elliptic functions are defined in Abramowitz &…

* Document Jacobi eigen function, in particular that it only works for symmetric matrices.

function [y,n]=jacobi(A,b,x0,eps) if nargin==3 eps=1.0e-6; elseif nargin=eps x0=y; y=B*x0+f; n=n+1; end

/* * Compute the Jacobi symbol (small prime case only). * * Copyright (c) 1995 Colin Plumb. All rights reserved. * For licen…

/* * For a small (usually prime, but not necessarily) prime p, * Return Jacobi(p,bn), which is -1, 0 or +1. * bn must be odd. …

function [x,n]=jacobi(A,b,x0,eps,varargin) if nargin==3 eps= 1.0e-6; M = 200; elseif nargin=eps x0=x; x=B*…

/*源文件jacobi_test.c*/ #include #include"jacobi_fun.h" void main() { int i; double eps; static double a[10][10]={{4,-1,0,0…

function X = jacobi(A,B,X0,kmax) %This function finds a soltuion to Ax = B by Jacobi iteration. if nargin < 4 tol = 1e-6; kmax =…

%do_gauss A = [3 2;1 2]; b = [1 -1]'; % x = gauss(A,b) % x = jacobi(A,b) % x = gauseid(A,b) x = SOR(A,b,1.0) x1 = A\b %for compari…