﻿功能:第三类不完全椭圆积分

格式:
EllipticPI(n, p, m)
EllipticPI(n, m)

n : 偏心参数,可以为矩阵变量也可以为数值.
p : 振幅,可以为矩阵变量也可以为数值. 一般此值在[0,1]之间取值
m : 为模数的平方,可以为矩阵变量也可以为数值,此参数值介于(0,1)之间. 

说明:
1. 这里的第三类不完全椭圆积分定义为
$$
\mathbf{EllipticPI}(n, p, m) = \int_0^p \dfrac{1}{(1-n\sin^2(t))\sqrt{1-m\sin^2(t)}}dt
$$

2. 这里的第三类完全椭圆积分定义为
$$
\mathbf{EllipticPI}(n, m) = \int_0^1 \dfrac{1}{(1-n\sin^2(t))\sqrt{1-m\sin^2(t)}}dt
$$

3. 在使用过程中, 注意定义中各个变量的正负号.

例子:
t = EllipticPI(1.5, 0.25, 0.36)//回车得到如下结果
t =
[ 0.25916284529762 ]