图书馆( 千斤顶 )
日本石油公司 <- 杰克波尔( 2 , λ = c(c)( 三 , 1 ), 阿尔法 = " 2/5 " ) 日本石油公司
# #98/25*x^3.y+28/5*x^2.y^2+98/25*x.y^3
evalQ喷雾( 日本石油公司 ,c( " 2 " , " 3/2 " )) # #大理性(bigq):
# # [1] 1239/10
贾卡斯 <- as.功能( 日本石油公司 )
贾卡斯( 2 , " 3/2 " ) # # [1] "1239/10"
贾卡斯( " x个 " , " x个 " ) # #[1]英寸(336*x^4)/25英寸
贾卡斯( " 2+2*I " , " 2/3+I/4 " ) # #[1]“综合体(-158921)/2160101689/2160)”
evalQ喷雾( 日本石油公司 , 值_re = c(c)( " 2 " , " 2/3 " ), 值_im = c(c)( " 2 " , " 1/4 " )) # #长度为2的Big Rational(“bigq”)对象:
# # [1] -158921/2160 101689/2160
千斤顶(c( " 2 " , " 3/2 " ), λ = c(c)( 三 , 1 ), 阿尔法 = " 2/5 " ) # #大理性(bigq):
# # [1] 1239/10
( J型 <- JackSymPol公司( 2 , λ = c(c)( 三 , 1 ))) # #{[2*a^2+4*a+2]}*X^3.Y+{[4*a+4]}*X ^2.Y^2+{[2*a ^2+4*a+2]{*X.Y^3
( J5号机组 <- substitute参数( J型 , 5 ) ) # #72*X^3.Y+24*X^2.Y^2+72*X.Y^3
J5号机组 == 杰克波尔( 2 , λ = c(c)( 三 , 1 ), 阿尔法 = " 5 " ) # #[1]正确
显示SymbolicQsprayOption( J型 , " 一 " ) <- " 阿尔法 " 显示符号QsrayOption( J型 , " X(X) " ) <- " x个 "
J型
# #{[2*alpha^2+4*alpha+2]}*x1^3.x2+{[4*alpha+4]}*x1^2.x2^2+{[2*alpha^2+4*alpha+2]}*x1.x2^3
显示SymbolicQsprayOption( J型 , " show单项式 " ) <- 显示单项式XYZ(c( " x个 " , " 年 " )) J型
# #{[2*alpha^2+4*alpha+2]}*x^3.y+{[4*alpha+4]}*x^2.y^2+{[2*alpha^2+4*alfa+2]}*x.y^3
( J型 <- 杰克波尔( 三 , λ = c(c)( 4 , 三 , 1 ), 阿尔法 = " 2 " ) ) # #3888*x^4.y^3.z+2592*x^3.y^2.z^2+3888*x^2.y^3+3888*x^3.y ^4.z+4752*x|3.y^3.z ^2+4752*x^3.y^2.z ^3+3888*x ^3.y.z ^4+2592*x ^2.y ^4.z ^2+4752*x^2.y^3.y ^3+3888*x.y^3.z^4 紧凑对称Qspray( J型 ) | > cat() # #3888*M[4,3,1]+2592*M[4,2,2]+4752*M[3]
λ <- c(c)( 三 , 1 ) 阿尔法 <- " 三 "
J4号机组 <- 杰克波尔( 4 , λ , 阿尔法 ) 第9页 <- JackPol公司( 9 , λ , 阿尔法 ) 紧凑对称Qspray( J4号机组 ) | > cat() # #32*M[3]+16*M[2,2]+28*M[2,1,1]+24*M[1,1,1,1] 紧凑对称Qspray( J9号机组 ) | > 猫() # #32*M[3]+16*M[2,2]+28*M[2,1,1]+24*M[1,1,1,1]
阿尔法 <- " 三 "
J1号机组 <- 杰克波尔( 4 , λ = c(c)( 三 , 1 ), 阿尔法 , 哪一个 = " P(P) " ) J2型 <- 杰克波尔( 4 , λ = c(c)( 2 , 2 ), 阿尔法 , 哪一个 = " P(P) " ) HallInnerProduct公司( J1号机组 , J2型 , 阿尔法 ) # #大理性(bigq):
# # [1] 0 HallInnerProduct公司( J1号机组 , J1号机组 , 阿尔法 ) # #大理性(bigq):
# # [1] 135/8 HallInner产品( J2型 , J2型 , 阿尔法 ) # #大理性(bigq):
# # [1] 63/5
HallInner产品( J1号机组 , J1号机组 , 阿尔法 = 无效的 ) # #3/128*α^4+1/4*α^3+63/128*α^2+81/64*α
J型 <- JackSymPol公司( 4 , λ = c(c)( 三 , 1 ), 哪一个 = " P(P) " ) 显示SymbolicQsprayOption( J型 , " 一 " ) <- " t吨 " HallInnerProduct公司( J型 , J型 , 阿尔法 = 2 ) # #[20*t^4-24*t^3+92*t^2-48*t+104]%//%[t^4+4*t^3+6*t^2+4*t+1]
( 臀部 <- HallInnerProduct公司( J型 , J型 , 阿尔法 = 无效的 ) ) # #{[24]%//%[4*t^4+16*t^3+24*t^2+16*t1+4]}*α^4+{[9*t^2-6*t+1]%//%[t^4+4*t^3+6*t^2+4*t+1]}*γ^3+{[3*t^4-6*t^3+5*t^2]%//%[t^4+4*t63+6*t+2]}*β^2+{[4*t ^4]%//%[t^4+4*t^3+6*t^2+4*t+1]}*alpha
changeVariables(更改变量)( 臀部 , 列表 (qlone( 1 ))) # #[3*t^4+t^3]%//%[t^2+2*t+1]
Laplace Beltrami公司 <- 功能 ( q喷雾 , 阿尔法 ) { n个 <- 变量数( q喷雾 ) 衍生产品1 <- 重叠(seq_len( n个 ), 功能 ( 我 ) { 衍生Q喷雾( q喷雾 , 我 ) }) 衍生产品2 <- 重叠(seq_len( n个 ), 功能 ( 我 ) { 衍生Q喷雾( 衍生产品1 [[ 我 ]], 我 ) }) x个 <- 重叠(seq_len( n个 ), qlone公司 ) # x_1、x_2。。。, x个n
# 第一学期
输出1 <- 0升
对于 ( 我 在里面 seq_len( n个 )) { 输出1 <- 输出1 + 阿尔法 * x个 [[ 我 ]] ^ 2 * 衍生产品2 [[ 我 ]] } # 第二任期
输出2 <- 0升
对于 ( 我 在里面 序列号(_L)( n个 )) { 对于 ( j个 在里面 seq_len( n个 )) { 如果 ( 我 != j个 ) { 输出2 <- 输出2 + x个 [[ 我 ]] ^ 2 * 衍生产品1 [[ 我 ]] / ( x个 [[ 我 ]] - x个 [[ j个 ]]) } } } # 在这个步骤中,“out2”是一个“ratioOfQsprays”对象,因为存在除法
# 按`x[[i]]-x[[j]]`; 但实际上它的分母是1,因为
# 简化,然后提取其分子以获得“qspray”对象
输出2 <- 获取分子( 输出2 ) 输出1 / 2 + 输出2 }
阿尔法 <- " 三 "
J型 <- 杰克波尔( 4 ,c( 2 , 2 ), 阿尔法 ) 共线Q喷雾( q喷雾1 = Laplace Beltrami公司( J型 , 阿尔法 ), q喷雾2 = J型 ) # #[1]正确