23.12.2012 Views

MATLAB Function Reference Volume 1: A - E - Bad Request

MATLAB Function Reference Volume 1: A - E - Bad Request

MATLAB Function Reference Volume 1: A - E - Bad Request

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

2ellipj<br />

Purpose Jacobi elliptic functions<br />

Syntax [SN,CN,DN] = ellipj(U,M)<br />

[SN,CN,DN] = ellipj(U,M,tol)<br />

Definition The Jacobi elliptic functions are defined in terms of the integral:<br />

u<br />

Then<br />

=<br />

∫<br />

dθ<br />

1 msin θ 2<br />

φ<br />

------------------------------------ 1<br />

0<br />

--<br />

( – ) 2<br />

sn( u)<br />

= sinφ, cn( u)<br />

= cosφ,<br />

dn( u)<br />

1 msin φ 2<br />

= ( – ) 2,<br />

am( u)<br />

= φ<br />

Some definitions of the elliptic functions use the modulus k instead of the<br />

parameter m . They are related by<br />

k2 m sin α 2<br />

= =<br />

The Jacobi elliptic functions obey many mathematical identities; for a good<br />

sample, see [1].<br />

Description [SN,CN,DN] = ellipj(U,M) returns the Jacobi elliptic functions SN, CN, and<br />

DN, evaluated for corresponding elements of argument U and parameter M.<br />

Inputs U and M must be the same size (or either can be scalar).<br />

[SN,CN,DN] = ellipj(U,M,tol) computes the Jacobi elliptic functions to<br />

accuracy tol. The default is eps; increase this for a less accurate but more<br />

quickly computed answer.<br />

Algorithm ellipj computes the Jacobi elliptic functions using the method of the<br />

arithmetic-geometric mean [1]. It starts with the triplet of numbers:<br />

1<br />

1<br />

--<br />

--<br />

2<br />

2<br />

a0 = 1, b0 = ( 1 – m)<br />

, c0 =<br />

( m)<br />

1<br />

--<br />

ellipj<br />

2-481

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!