R/tobit2derivatives.R

Defines functions fl1 fl2 fl3 fl4

fl1 = function(y, mu1, mu2, sigma, rho, censoring, eps)
{
    .e1 <- -mu2
tt = matrix(    c(rep(0,length(y)), 2 * (dnorm(.e1, 0, 1)/pnorm(.e1)), rep(0,length(y)), 
        rep(0,length(y))), nrow = length(y))

##
    .e2 <- 1 + eps - rho^2
    .e3 <- y - mu1
    .e4 <- sqrt(.e2)
    .e5 <- mu2 + rho * .e3/sigma
    .e6 <- -(.e5/.e4)
    .e7 <- 1 - pnorm(.e6)
    .e8 <- dnorm(.e6, 0, 1)
    .e9 <- .e7 * .e4
    .e10 <- .e3/sigma
    .e11 <- rho * .e8
    .e12 <- sigma^2
    matrix(c((2 * (.e11/.e9) - 2 * .e10)/sigma, -(2 * (.e8/.e9)), 
        (2 * (.e11/(.e12 * .e7 * .e4)) - 8 * (sigma * 
            .e3/(2 * .e12)^2)) * .e3 + 2/sigma, -(2 * ((.e10 + 
            rho * .e5/.e2) * .e8/.e9))), nrow = length(y)) -> t2
t2[censoring,] = tt[censoring,]
t2

}
fl2 = function(y, mu1, mu2, sigma, rho, censoring, eps)
{
    .e1 <- -mu2
    .e2 <- dnorm(.e1, 0, 1)
    .e3 <- pnorm(.e1)
 tt = matrix(   c(c(rep(0,length(y)), rep(0,length(y)), rep(0,length(y)), rep(0,length(y))), c(rep(0,length(y)), 2 * (.e2 * 
        (.e2/.e3 - mu2)/.e3), rep(0,length(y)), rep(0,length(y))), c(rep(0,length(y)), 
    rep(0,length(y)), rep(0,length(y)), rep(0,length(y))), c(rep(0,length(y)), rep(0,length(y)), rep(0,length(y)), rep(0,length(y)))), nrow = length(y))



    .e1 <- rho^2
    .e3 <- 1 + eps - .e1
    .e4 <- y - mu1
    .e5 <- sqrt(.e3)
    .e6 <- mu2 + rho * .e4/sigma
    .e7 <- -(.e6/.e5)
    .e8 <- 1 - pnorm(.e7)
    .e9 <- dnorm(.e7, 0, 1)
    .e10 <- rho * .e6
    .e11 <- sigma^2
    .e12 <- .e8 * .e5
    .e13 <- .e4/sigma
    .e14 <- .e10/.e3
    .e15 <- .e13 + .e14
    .e16 <- .e12^2
    .e18 <- .e8 * .e3 * .e5
    .e19 <- .e11 * .e8
    .e21 <- .e6/.e18 + .e9/.e16
    .e22 <- .e19 * .e5
    .e23 <- .e15 * .e9
    .e24 <- .e15 * .e6
    .e25 <- .e23 - rho * .e8/.e5
    .e26 <- .e22^2
    .e27 <- rho * .e15
    .e28 <- rho * .e9
    .e29 <- (2 * .e11)^2
    .e31 <- rho * .e21 * .e9
    .e32 <- -(2 * (.e31/sigma))
    .e33 <- ((rho * (.e24 - rho)/.e3 - 1)/.e12 + .e27 * .e9/.e16) * 
        .e9
    .e34 <- 1 - .e27 * .e6/.e3
    .e35 <- 2 * (.e1 * .e21 * .e9)
    .e36 <- sigma^4
    matrix(c(c((2 + .e35)/.e11, .e32, (16 * 
        (sigma/.e29) + 2 * (.e1 * (.e6/(sigma^3 * .e8 * .e3 * 
        .e5) + sigma * .e9/.e26) * .e9)) * .e4 - 2 * (.e28/.e22), 
        -(2 * (.e33/sigma))), c(.e32, 2 * 
        (.e21 * .e9), -(2 * (rho * (.e6/(.e19 * .e3 * 
        .e5) + .e11 * .e9/.e26) * .e9 * .e4)), -(2 * (((rho - 
        .e24)/.e18 - .e23/.e16) * .e9))), c(((.e35 + 
        4) * .e4/sigma - 2 * (.e28/.e12))/.e11, -(2 * (.e31 * 
        .e4/.e11)), (2 * (.e28 * (.e10 * .e4/(.e36 * 
        .e8 * .e3 * .e5) - (2 * (sigma * .e8) - .e28 * .e4/.e5) * 
        .e5/.e26)) - 8 * ((1 - 16 * (.e36/.e29)) * .e4/.e29)) * 
        .e4 - 2/.e11, -(2 * (.e33 * .e4/.e11))), c(2 * 
        ((.e34/.e12 - rho * .e25/.e16) * .e9/sigma), 2 * 
        ((.e25/.e16 + .e24/.e18) * .e9), 2 * ((.e34/.e22 - 
        rho * .e11 * .e25/.e26) * .e9 * .e4), -(2 * (((mu2 + 
        rho * (2 * .e13 + 2 * .e14) - .e15^2 * .e6)/.e18 - .e25 * 
        .e15/.e16) * .e9)))), nrow = length(y)) -> t2
t2[censoring,] = tt[censoring,]
t2

}


fl3 = function(y, mu1, mu2, sigma, rho, censoring, eps)
{
    .e1 <- -mu2
    .e2 <- dnorm(.e1, 0, 1)
    .e3 <- pnorm(.e1)
    .e4 <- .e2/.e3
tt=matrix(    c(rep(0,21*length(y)) , 
        2 * (((2 * .e4 - mu2) * (.e4 - mu2) - 1) * .e2/.e3), 
       
rep(0,42*length(y))), nrow = length(y))


    .e1 <- rho^2
    .e3 <- 1 + eps - .e1
    .e4 <- y - mu1
    .e5 <- sqrt(.e3)
    .e6 <- rho * .e4
    .e7 <- .e6/sigma
    .e8 <- mu2 + .e7
    .e9 <- -(.e8/.e5)
    .e10 <- pnorm(.e9)
    .e11 <- 1 - .e10
    .e12 <- dnorm(.e9, 0, 1)
    .e13 <- rho * .e8
    .e14 <- .e4/sigma
    .e15 <- .e13/.e3
    .e16 <- .e14 + .e15
    .e17 <- .e11 * .e5
    .e18 <- .e17^2
    .e19 <- sigma^2
    .e21 <- .e11 * .e3 * .e5
    .e22 <- rho * .e11
    .e23 <- .e19 * .e11
    .e24 <- .e23 * .e5
    .e25 <- .e16 * .e12
    .e26 <- .e24^2
    .e27 <- .e16 * .e8
    .e28 <- .e25 - .e22/.e5
    .e29 <- rho * .e16
    .e30 <- .e21^2
    .e31 <- .e12/.e18
    .e33 <- .e8/.e21 + .e31
    .e34 <- 2 * .e15
    .e35 <- sigma^4
    .e36 <- rho * .e12
    .e37 <- sigma^3
    .e38 <- sigma * .e11
    .e39 <- .e8/.e3
    .e41 <- .e29 * .e8/.e3
    .e42 <- .e27 - rho
    .e43 <- 1/.e17
    .e44 <- 2 * .e14
    .e45 <- .e3 * .e8
    .e47 <- mu2 + rho * (.e44 + .e34)
    .e48 <- ((.e39 + 2 * (.e11 * .e12 * .e5/.e18))/.e18 + .e45/.e30) * 
        .e12
    .e49 <- .e33 * .e8
    .e50 <- .e36 * .e4
    .e51 <- .e37 * .e11
    .e52 <- .e35 * .e11
    .e53 <- .e50/.e5
    .e54 <- .e1/.e3
    .e55 <- .e27/.e3
    .e56 <- .e16^2
    .e57 <- 1 - .e41
    .e58 <- .e56 * .e8
    .e60 <- 2 * .e38 - .e53
    .e61 <- rho * .e42
    .e62 <- .e51 * .e3
    .e63 <- (1 + .e54) * .e8
    .e65 <- .e52 * .e3 * .e5
    .e66 <- .e16 * .e3
    .e67 <- .e47 - .e58
    .e69 <- .e61/.e3 - 1
    .e71 <- (.e49 - .e43)/.e3 + .e48
    .e73 <- (.e66 * .e12/.e5 - 2 * .e22) * .e5 - .e22 * .e3/.e5
    .e77 <- .e63 + .e29
    .e78 <- 2 * (.e28 * .e11 * .e5/.e18)
    .e83 <- .e69/.e17
    .e84 <- 2 * (.e16 * .e11 * .e12 * .e5/.e18)
    .e85 <- .e23 * .e3
    .e86 <- .e77/.e17
    .e87 <- rho - .e27
    .e88 <- .e13 * .e4
    .e89 <- .e85 * .e5
    .e90 <- .e62 * .e5
    .e91 <- .e28 * .e16
    .e92 <- .e55 + .e78
    .e93 <- .e83 + .e29 * .e12/.e18
    .e94 <- .e14 + .e34
    .e95 <- 2 * .e19
    .e96 <- rho * .e71
    .e97 <- .e19 * .e3
    .e98 <- .e19 * .e12
    .e99 <- .e95^2
    .e100 <- sigma * .e3
    .e101 <- .e98/.e26
    .e102 <- .e73 * .e8
    .e103 <- .e42/.e3
    .e105 <- .e92 * .e12/.e18
    .e106 <- .e93 * .e8
    .e107 <- .e25/.e18
    .e108 <- .e57/.e17
    .e109 <- .e57/.e24
    .e111 <- (.e43 - .e49)/.e3 - .e48
    .e115 <- .e8/.e89 + .e101
    .e116 <- .e4/(.e38 * .e3 * .e5)
    .e118 <- .e96 * .e4/sigma
    .e119 <- rho * .e28
    .e120 <- rho * .e94
    .e122 <- .e88/.e65 - .e60 * .e5/.e26
    .e125 <- .e91 * .e11
    .e126 <- .e28/.e18
    .e128 <- .e16 * .e4/sigma
    .e129 <- .e108 - .e119/.e18
    .e130 <- .e109 - rho * .e19 * .e28/.e26
    .e131 <- .e47 * .e8
    .e133 <- .e8/.e90 + sigma * .e12/.e26
    .e135 <- 2 * .e55
    .e136 <- 2 * (.e62 * .e60/.e26)
    .e137 <- rho * (.e116 - (.e102/.e30 + .e105))
    .e139 <- .e35 * .e28 * .e11
    .e141 <- (.e106 - .e86)/.e3 + .e12 * (rho * (2 * .e103 + 
        .e84) - 2)/.e18
    .e145 <- .e126 + .e27/.e21
    .e146 <- (.e39 + 2 * (.e52 * .e12 * .e5/.e26))/.e26
    .e150 <- .e67/.e21 - .e91/.e18
    .e152 <- .e87/.e21 - .e107
    .e153 <- .e89^2
    .e154 <- .e90^2
    .e155 <- .e65^2
    .e156 <- 16 * (.e35/.e99)
    .e157 <- 2 * (.e125 * .e5/.e18)
    .e159 <- 2 * (.e139 * .e5/.e26)
    .e161 <- 2 * mu2 + rho * (.e34 + 3 * .e14)
    .e163 <- .e12 * .e67/.e3
    .e164 <- rho * (.e25/.e5 + .e22/.e3)
    .e165 <- rho * .e111
    .e166 <- -(2 * (.e36 * (.e118 - .e33)/.e19))
    .e167 <- .e28^2
    .e168 <- .e129 * .e8
    .e169 <- .e130 * .e8
    .e170 <- ((.e43 - rho * .e33 * .e16) * .e8/.e3 + .e31 + .e137) * 
        .e12
    .e171 <- .e150 * .e8
    .e174 <- .e3 * .e12 * .e67/.e30
    .e175 <- .e131/.e3
    .e176 <- 1 - .e156
    .e177 <- 1/.e24
    .e178 <- 2 - .e136
    .e179 <- 2 - .e41
    .e181 <- 2 * .e54
    .e182 <- .e165 * .e12
    .e183 <- .e120/.e17
    .e184 <- .e41 - 1
    .e185 <- rho * (.e78 + .e135)
    .e186 <- .e51 * .e5
    .e187 <- -(2 * (.e170/sigma))
    .e188 <- -(2 * (.e12 * (rho * (.e141 * .e4/sigma - .e107) - 
        .e83)/.e19))
    .e189 <- -(2 * (.e182/sigma))
    .e190 <- -(2 * (.e1 * .e71 * .e12/.e19))
    .e191 <- (((.e27 + rho * ((.e131 + 2 * .e61)/.e3 + .e128 - 
        2))/.e17 - .e93 * .e16 * .e8)/.e3 + ((.e14 + rho * ((.e161 - 
        .e58)/.e3 - .e157)) * .e12 - .e28 * .e69)/.e18) * .e12
    .e192 <- ((.e86 + rho * .e152 * .e8)/.e3 + .e12 * (rho * 
        .e3 * .e87/.e30 - (rho * (.e103 + .e84) - 1)/.e18)) * 
        .e12
    .e193 <- .e73/.e30
    .e194 <- ((.e16 * (2 * .e63 + .e29) - rho * (2 + .e181))/.e21 + 
        rho * (.e171/.e3 + .e174) - (.e28 * (rho * (.e84 - rho/.e3) - 
        1) + .e25 * .e184)/.e18) * .e12
    .e196 <- (.e168 + .e86)/.e3 + (2 - .e185) * .e12/.e18
    .e197 <- (((rho * .e92 - 1)/.e18 + .e29 * .e3 * .e8/.e30) * 
        .e12 + (rho * .e145 * .e8 - .e86)/.e3) * .e12
    .e199 <- (.e55 + .e159) * .e12/.e26
    .e201 <- .e28 * .e57 + .e25
    .e205 <- .e169 + .e77/.e24
    .e206 <- ((1 + rho * ((2 * rho - 2 * .e27)/.e3 - .e84)) * 
        .e12/.e18 + (.e183 - .e106)/.e3) * .e12
    .e209 <- .e146 + .e45/.e153
    .e210 <- .e146 + .e97 * .e8/.e154
    .e211 <- .e115 * .e8
    .e212 <- .e133 * .e8
    .e213 <- .e94/.e17
    .e215 <- (.e14 + rho * .e161/.e3) * .e8 + .e120 * .e16
    .e216 <- .e128 + .e175
    .e218 <- (1 + .e164 - .e10)/.e5
    .e219 <- (2 + .e164 - 2 * .e10)/.e5
    .e220 <- .e8 * .e122
    .e221 <- 1 + .e88/.e100
    .e222 <- 1/.e186
    .e223 <- 16 * (.e176/.e99)
    .e224 <- 2 * (.e167 * .e11 * .e5/.e18)
    .e226 <- 2 * .e33
    .e230 <- 2 * .e7 + mu2
    .e231 <- 2 * (.e36 * .e122)
    .e232 <- 2 * (.e1 * .e33 * .e12)
    .e234 <- .e12 * (.e88/.e3 + sigma * .e178)/.e26
    .e235 <- .e96 * .e12
    .e236 <- rho * .e141
    .e237 <- rho * (.e33 * .e179 + .e137)
    .e238 <- rho * .e115
    .e239 <- rho * (.e135 + .e159)
    .e242 <- .e1 * (2 * .e118 - 4 * .e33) * .e12 - 4
    .e245 <- rho^3 * .e71 * .e12
    matrix(c(c(c(2 * (.e245/.e37), .e190, .e1 * 
        (2 * (rho * (.e210 * .e12 + (.e212 - .e222)/.e100) * 
            .e4) - 4 * .e133) * .e12 - 16 * (sigma/.e99), -(2 * 
        (.e236 * .e12/.e19))), c(.e190, 2 * 
        (.e235/sigma), -(2 * (.e36 * (rho * ((.e211 - 
        .e177)/.e100 + sigma * .e209 * .e12) * .e4 - .e115))), 
        -(2 * (.e192/sigma))), c(.e242/.e37, 
        .e166, (.e223 + 2 * (.e1 * ((.e220 - .e230/.e186)/.e100 + 
            .e12 * (rho * .e37 * .e3 * .e8 * .e4/.e155 - (.e136 - 
                .e221)/.e26)) * .e12)) * .e4 - .e231, .e188), 
        c(2 * (rho * .e196 * .e12/.e19), 2 * 
            (.e197/sigma), 2 * (.e12 * (rho * ((.e205/.e100 + 
            sigma * (2 - .e239) * .e12/.e26) * .e4 + .e19 * .e28/.e26) - 
            .e109)), -(2 * (.e194/sigma)))), c(c(2 * 
        (.e1 * .e111 * .e12/.e19), .e189, rho * 
        (2 * .e115 + 2 * (rho * ((.e222 - .e212)/.e3 - sigma * 
            .e210 * .e12) * .e4)) * .e12, -(2 * (.e206/sigma))), 
        c(.e189, 2 * (.e111 * .e12), -(2 * 
            (rho * ((.e177 - .e211)/.e3 - .e19 * .e209 * .e12) * 
                .e12 * .e4)), 2 * (((.e152 * .e8 + .e213)/.e3 + 
            ((.e87/.e3 - .e84)/.e18 + .e3 * .e87/.e30) * .e12) * 
            .e12)), c(rho * (.e226 + 2 * (.e165 * 
            .e4/sigma)) * .e12/.e19, -(2 * (.e182 * .e4/.e19)), 
            2 * (.e36 * (rho * (1/.e65 - .e35 * .e3 * 
                .e12 * .e8/.e155) * .e4 - (.e220/.e3 + .e234)) * 
                .e4), -(2 * (.e206 * .e4/.e19))), c(2 * 
            (.e12 * (.e12 * (.e185 - 1)/.e18 - (.e168 + .e183)/.e3)/sigma), 
            2 * (((.e213 - .e145 * .e8)/.e3 - (.e92/.e18 + 
                .e66 * .e8/.e30) * .e12) * .e12), 2 * 
                (.e12 * (.e98 * (.e239 - 1)/.e26 - (.e169 + .e120/.e24)/.e3) * 
                  .e4), -(2 * ((((1 + .e181 - (.e14 + 3 * 
                .e15) * .e16)/.e17 - .e171)/.e3 - (((.e119 - 
                .e56 * .e12 * .e8)/.e3 - 2 * (.e125 * .e12 * 
                .e5/.e18))/.e18 + .e174)) * .e12)))), c(c((2 * 
        (.e245 * .e4/sigma) - 2 * (2 + .e232))/.e37, .e166, 
        (.e223 + 2 * (.e1 * ((.e221 - .e136) * .e12/.e26 + 
            .e8 * (rho * .e133 * .e4/.e97 - .e100 * (3 * .e38 - 
                .e53) * .e5/.e154) - .e6/(sigma^5 * .e11 * .e3 * 
            .e5)) * .e12)) * .e4 - .e231, .e188), c(.e166, 
        2 * (.e235 * .e4/.e19), -(2 * (rho * (.e8 * 
            (.e238 * .e4/.e97 - .e3 * .e60 * .e5/.e153) + .e234 - 
            .e6/.e65) * .e12 * .e4)), -(2 * (.e192 * .e4/.e19))), 
        c((.e242 * .e4/sigma - 2 * ((.e232 + 4) * 
            .e4/sigma - 2 * (.e36/.e17)))/.e37, -(2 * (.e36 * 
            (.e118 - .e226) * .e4/.e37)), (2 * (.e36 * 
            (rho * ((.e122/.e97 - .e97 * (4 * .e38 - .e53) * 
                .e5/.e155) * .e8 - .e6/(sigma^6 * .e11 * .e3 * 
                .e5)) * .e4 - (2 * (1 - (.e10 + .e50/(sigma * 
                .e5))) - (2 * (.e85 * .e60^2/.e26) + .e1 * .e12 * 
                .e8 * .e4^2/(.e97 * .e5))) * .e5/.e26)) + 8 * 
            (.e37 * (16 * .e176 + 16 * (4 - .e156)) * .e4/.e95^4)) * 
            .e4 + 4/.e37, -(2 * (.e12 * (.e236 * .e4/sigma - 
            2 * .e93) * .e4/.e37))), c(2 * (.e12 * 
            (rho * (.e196 * .e4/sigma + .e126) - .e108)/.e19), 
            2 * (.e197 * .e4/.e19), 2 * (.e12 * 
                (rho * (.e205 * .e4/.e97 - (.e12 * .e184 * .e4 + 
                  sigma * .e28 * .e178)/.e26) - .e57 * .e60 * 
                  .e5/.e26) * .e4), -(2 * (.e194 * .e4/.e19)))), 
        c(c(2 * (.e237 * .e12/.e19), .e187, 
            (2 * (rho * (.e133 * .e179 + rho * (.e4/.e65 - 
                sigma * (.e199 + .e19 * .e73 * .e8/.e154))) * 
                .e4) - 2 * .e130) * .e12, -(2 * (.e191/sigma))), 
            c(.e187, 2 * ((.e116 - ((.e193 + 
                .e33 * .e16/.e3) * .e8 + .e105)) * .e12), -(2 * 
                (((.e177 - .e238 * .e16) * .e8/.e3 + rho * (.e4/.e90 - 
                  .e19 * (.e102/.e153 + .e199)) + .e101) * .e12 * 
                  .e4)), -(2 * ((((1 - .e216)/.e17 - .e152 * 
                .e16 * .e8)/.e3 - (.e73 * .e87/.e30 + (.e67/.e3 - 
                .e157) * .e12/.e18)) * .e12))), c((2 * 
                (.e237 * .e4/sigma) - 2 * .e129) * .e12/.e19, 
                -(2 * (.e170 * .e4/.e19)), 2 * 
                  ((.e57 * .e122 + rho * ((.e230/.e65 - rho * 
                    .e35 * .e73 * .e8/.e155) * .e4 - ((2 * (sigma * 
                    .e16) - (1 + rho * .e87/.e3) * .e4) * .e12 - 
                    .e60 * (2 * (.e139 * .e3/.e26) + rho/.e5))/.e26)) * 
                    .e12 * .e4), -(2 * (.e191 * .e4/.e19))), 
            c(-(2 * (((.e215/.e17 + .e129 * .e16 * 
                .e8)/.e3 + (.e201 + rho * (.e163 - (.e219 + .e224)))/.e18) * 
                .e12/sigma)), 2 * ((.e216/.e21 + (.e163 - 
                (.e218 + .e224))/.e18 - (.e193 + .e145/.e3) * 
                .e16 * .e8) * .e12), -(2 * (((.e215/.e24 + 
                .e130 * .e16 * .e8)/.e3 + .e19 * (.e201 + rho * 
                (.e163 - (.e219 + 2 * (.e35 * .e167 * .e11 * 
                  .e5/.e26))))/.e26) * .e12 * .e4)), -(2 * 
                ((((.e44 + rho * (2 * .e47 + 2 * .e8)/.e3 - (.e128 + 
                  2 * .e175) * .e16)/.e17 - .e150 * .e16 * .e8)/.e3 - 
                  (.e73 * .e67/.e30 + (.e28 * (.e47/.e3 - .e157) + 
                    .e16 * (.e163 - .e218))/.e18)) * .e12))))), nrow = length(y)) -> t2
t2[censoring,] = tt[censoring,]
t2

}
fl4 = function(y, mu1, mu2, sigma, rho, censoring, eps)
{
    .e1 <- -mu2
    .e2 <- dnorm(.e1, 0, 1)
    .e3 <- pnorm(.e1)
    .e4 <- .e2/.e3
    .e5 <- .e4 - mu2

tt = matrix(c(rep(0, 85*length(y)),
2 * ((((4 * .e4 - 
            mu2) * .e5 - 2) * .e5 + (2 * .e4 - mu2) * (.e2 * 
            .e5/.e3 - 1)) * .e2/.e3), rep(0, 170*length(y))), nrow = length(y))




    .e1 <- rho^2
    .e2 <- 1 + eps
    .e3 <- .e2 - .e1
    .e4 <- y - mu1
    .e5 <- sqrt(.e3)
    .e6 <- rho * .e4
    .e7 <- .e6/sigma
    .e8 <- mu2 + .e7
    .e9 <- -(.e8/.e5)
    .e10 <- pnorm(.e9)
    .e11 <- 1 - .e10
    .e12 <- .e4/sigma
    .e13 <- rho * .e8
    .e14 <- .e13/.e3
    .e15 <- dnorm(.e9, 0, 1)
    .e16 <- .e12 + .e14
    .e17 <- .e11 * .e5
    .e18 <- .e17^2
    .e19 <- rho * .e11
    .e20 <- sigma^2
    .e21 <- .e16 * .e15
    .e22 <- .e11 * .e3
    .e23 <- .e22 * .e5
    .e24 <- .e20 * .e11
    .e25 <- .e19/.e5
    .e26 <- .e21 - .e25
    .e27 <- .e24 * .e5
    .e28 <- .e16 * .e8
    .e29 <- .e27^2
    .e30 <- rho * .e16
    .e31 <- 2 * .e14
    .e32 <- .e23^2
    .e33 <- .e1/.e3
    .e34 <- 2 * .e12
    .e35 <- sigma^4
    .e36 <- .e34 + .e31
    .e37 <- rho * .e36
    .e38 <- sigma * .e11
    .e39 <- mu2 + .e37
    .e40 <- 1 + .e33
    .e41 <- .e8/.e3
    .e42 <- rho * .e15
    .e43 <- sigma^3
    .e44 <- .e42 * .e4
    .e45 <- .e44/.e5
    .e46 <- .e15/.e18
    .e47 <- .e11 * .e15
    .e50 <- 2 * (.e47 * .e5/.e18)
    .e51 <- .e41 + .e50
    .e53 <- .e30 * .e8/.e3
    .e54 <- .e28 - rho
    .e55 <- .e43 * .e11
    .e56 <- .e40 * .e8
    .e57 <- .e16^2
    .e58 <- .e26 * .e11
    .e59 <- .e16 * .e3
    .e60 <- .e3 * .e8
    .e62 <- 2 * .e38 - .e45
    .e63 <- .e8/.e23
    .e65 <- .e59 * .e15/.e5
    .e66 <- .e63 + .e46
    .e67 <- .e35 * .e11
    .e68 <- .e65 - 2 * .e19
    .e69 <- .e57 * .e8
    .e70 <- .e68 * .e5
    .e71 <- .e55 * .e3
    .e72 <- .e28/.e3
    .e74 <- .e51/.e18 + .e60/.e32
    .e75 <- .e70 - .e19 * .e3/.e5
    .e76 <- .e56 + .e30
    .e77 <- .e39 - .e69
    .e79 <- .e58 * .e5/.e18
    .e80 <- 2 * .e79
    .e81 <- .e16 * .e11
    .e82 <- rho * .e54
    .e83 <- rho - .e28
    .e84 <- .e74 * .e15
    .e85 <- .e11^2
    .e86 <- .e24 * .e3
    .e89 <- .e81 * .e15 * .e5/.e18
    .e90 <- .e67 * .e3
    .e91 <- .e90 * .e5
    .e92 <- .e12 + .e31
    .e93 <- 2 * .e89
    .e94 <- sigma * .e3
    .e95 <- 1 - .e53
    .e96 <- rho * .e26
    .e97 <- .e86 * .e5
    .e98 <- .e20 * .e3
    .e99 <- .e71 * .e5
    .e100 <- .e26 * .e16
    .e101 <- .e13 * .e4
    .e102 <- .e82/.e3
    .e103 <- .e72 + .e80
    .e104 <- .e66 * .e8
    .e105 <- 1/.e17
    .e106 <- .e102 - 1
    .e108 <- .e11 * .e8/.e3
    .e109 <- .e15/.e5
    .e110 <- .e54/.e3
    .e111 <- .e35 * .e26
    .e112 <- .e39 * .e8
    .e113 <- .e3^2
    .e114 <- 3 * .e12
    .e115 <- 2 * mu2
    .e116 <- .e111 * .e11
    .e118 <- rho * (.e31 + .e114)
    .e119 <- .e53 - 1
    .e121 <- .e16 * .e4/sigma
    .e122 <- .e19/.e3
    .e123 <- 2 * .e33
    .e124 <- .e30 * .e15
    .e125 <- .e21/.e5
    .e126 <- .e11 * .e113
    .e127 <- .e126 * .e15
    .e128 <- 2 * (.e71 * .e62/.e29)
    .e129 <- .e106/.e17
    .e130 <- .e15 * .e77
    .e131 <- .e115 + .e118
    .e133 <- rho * (.e125 + .e122)
    .e134 <- .e97^2
    .e135 <- .e99^2
    .e136 <- .e91^2
    .e137 <- rho/.e5
    .e138 <- .e130/.e3
    .e139 <- .e100 * .e11
    .e140 <- .e41 + 2 * (.e67 * .e15 * .e5/.e29)
    .e141 <- 2 * rho
    .e142 <- .e129 + .e124/.e18
    .e144 <- .e38 * .e3 * .e5
    .e145 <- 2 * .e72
    .e146 <- rho * .e92
    .e147 <- .e20 * .e15
    .e148 <- .e75 * .e8
    .e149 <- .e76/.e17
    .e153 <- 2 * (.e127 * .e8 * .e5/.e32)
    .e155 <- .e116 * .e5/.e29
    .e156 <- .e104 - .e105
    .e157 <- .e112/.e3
    .e158 <- 2 * .e155
    .e160 <- .e103 * .e15/.e18
    .e163 <- 2 * (.e139 * .e5/.e18)
    .e164 <- 2 * .e46
    .e165 <- 2 * .e110
    .e166 <- .e109 - .e108
    .e167 <- .e26/.e18
    .e170 <- rho * (.e165 + .e93) - 2
    .e171 <- .e140/.e29
    .e172 <- .e108 - .e109
    .e173 <- .e147/.e29
    .e174 <- .e148/.e32
    .e175 <- (1 + .e133 - .e10)/.e5
    .e176 <- 2 + .e123
    .e178 <- .e156/.e3 + .e84
    .e179 <- 2 - .e128
    .e180 <- 1/.e3
    .e181 <- .e8/.e97
    .e182 <- .e62 * .e5
    .e183 <- .e4/.e144
    .e185 <- .e74 * .e8/.e3
    .e191 <- .e182/.e29
    .e192 <- .e181 + .e173
    .e193 <- 2 * (.e51 * .e11 * .e5)
    .e194 <- 2 * (.e85 * .e3 * .e15/.e18)
    .e195 <- 3 * .e14
    .e196 <- .e101/.e3
    .e197 <- sigma * .e15
    .e198 <- .e153 - 1
    .e199 <- .e101/.e94
    .e200 <- .e197/.e29
    .e201 <- .e21/.e18
    .e203 <- rho * .e131/.e3
    .e204 <- rho * .e20
    .e205 <- .e55 * .e5
    .e206 <- .e8/.e99
    .e207 <- 2 - .e53
    .e208 <- 2 * .e56
    .e211 <- .e95/.e17
    .e212 <- .e206 + .e200
    .e214 <- 2 * (.e58 * .e3/.e18) + .e137
    .e215 <- .e138 - .e175
    .e216 <- .e174 + .e160
    .e218 <- rho * (.e80 + .e145)
    .e220 <- .e204 * .e26/.e29
    .e221 <- .e95/.e27
    .e222 <- .e12 + .e203
    .e223 <- 2 * .e39
    .e224 <- rho * .e103
    .e225 <- .e83/.e23
    .e226 <- .e208 + .e30
    .e227 <- .e211 - .e96/.e18
    .e228 <- .e112 + 2 * .e82
    .e230 <- .e101/.e91 - .e191
    .e231 <- .e26^2
    .e232 <- .e142 * .e8
    .e233 <- .e222 * .e8
    .e234 <- .e221 - .e220
    .e235 <- .e77/.e23
    .e236 <- .e75 * .e11
    .e237 <- .e172 * .e5
    .e239 <- .e228/.e3 + .e121
    .e240 <- .e233 + .e146 * .e16
    .e241 <- .e183 - .e216
    .e242 <- 2 * .e84
    .e243 <- .e131 - .e69
    .e244 <- .e92/.e3
    .e245 <- .e16 * .e226
    .e246 <- .e28/.e23
    .e247 <- 1 + .e199
    .e248 <- 2 * .e8
    .e249 <- .e196 + sigma * .e179
    .e250 <- .e236 * .e3
    .e251 <- .e243/.e3
    .e252 <- .e100/.e18
    .e253 <- .e167 + .e246
    .e254 <- .e225 - .e201
    .e256 <- .e15 * .e170/.e18
    .e257 <- (.e185 + ((2 * (.e237 + .e194) + .e193) * .e15/.e18 - 
        .e180)/.e18 + .e3 * .e198/.e32) * .e15
    .e259 <- .e12 + rho * (.e251 - .e163)
    .e260 <- 2 * .e28
    .e261 <- .e235 - .e252
    .e262 <- rho * .e66
    .e264 <- .e28 + rho * (.e239 - 2)
    .e265 <- .e72 + .e158
    .e266 <- .e11 * .e119
    .e267 <- 1 + .e123
    .e269 <- 2 * .e7 + mu2
    .e270 <- 3/.e17
    .e271 <- rho * .e83
    .e272 <- .e16 * .e166
    .e273 <- .e121 + .e157
    .e274 <- 3 * rho
    .e275 <- .e96/.e5
    .e276 <- sigma * .e5
    .e279 <- rho * .e178 * .e4/sigma
    .e280 <- rho * .e3
    .e281 <- .e116 * .e3
    .e283 <- .e171 + .e60/.e134
    .e284 <- .e171 + .e98 * .e8/.e135
    .e285 <- .e21 * .e119
    .e286 <- (1 + .e271/.e3) * .e4
    .e287 <- 1 - (.e10 + .e44/.e276)
    .e288 <- 1/.e23
    .e290 <- 2 * (.e281/.e29) + .e137
    .e291 <- 2 * .e10
    .e292 <- .e141 - .e260
    .e293 <- rho * (.e110 + .e93)
    .e294 <- .e224 - 1
    .e295 <- rho/.e3
    .e296 <- sigma * .e16
    .e297 <- .e26 * .e106
    .e298 <- (2 + .e133 - .e291)/.e5
    .e299 <- .e105 - .e104
    .e303 <- 3 * .e38 - .e45
    .e304 <- rho * (.e145 + .e158)
    .e306 <- rho * (.e223 + .e248)/.e3
    .e307 <- .e292/.e3
    .e308 <- .e4^2
    .e309 <- 2 - .e218
    .e310 <- 2 * .e296
    .e312 <- 4 * .e38 - .e45
    .e313 <- .e293 - 1
    .e314 <- .e1 * .e15
    .e315 <- (.e266 - .e275) * .e5
    .e317 <- .e259 * .e15 - .e297
    .e318 <- .e34 + .e195
    .e319 <- 2 * (.e96 * .e85 * .e3/.e18)
    .e320 <- .e310 - .e286
    .e321 <- .e262 * .e16
    .e322 <- .e98 * .e5
    .e323 <- .e35 * .e85
    .e325 <- .e26 * .e95 + .e21
    .e330 <- 1 + rho * (.e307 - .e93)
    .e331 <- 2 * (.e231 * .e11 * .e5/.e18)
    .e334 <- rho * (.e93 - .e295) - 1
    .e338 <- .e76/.e23
    .e339 <- .e12 - 2 * (.e250 * .e8 * .e5/.e32)
    .e340 <- .e232 - .e149
    .e341 <- .e12 + .e195
    .e344 <- .e314 * .e8 * .e308/.e322
    .e346 <- .e265 * .e15/.e29
    .e348 <- .e100 * .e15/.e5
    .e351 <- .e299/.e3
    .e352 <- .e176 * .e8
    .e353 <- .e144^2
    .e354 <- 2 * (.e315 + .e319)
    .e355 <- .e141 - .e28
    .e356 <- rho * (.e28 + .e274)
    .e357 <- .e124/.e5
    .e358 <- rho * .e176
    .e359 <- .e75/.e32
    .e360 <- .e26 * .e334
    .e361 <- .e341 * .e16
    .e362 <- .e121 + 2 * .e157
    .e363 <- .e62^2
    .e365 <- .e294/.e18 + .e30 * .e3 * .e8/.e32
    .e367 <- 1 + .e357 - .e10
    .e368 <- .e34 + .e306
    .e369 <- 2 * .e287
    .e371 <- .e280 * .e83/.e32
    .e372 <- sigma * .e26
    .e373 <- .e272 * .e5
    .e375 <- .e3 * .e15
    .e376 <- .e320 * .e15
    .e377 <- .e166 * .e5
    .e379 <- (rho * .e172 * .e16 - .e11 * .e40) * .e5 + 2 * (.e30 * 
        .e85 * .e3 * .e15/.e18)
    .e380 <- 2 * .e76
    .e382 <- .e15 * .e249/.e29
    .e383 <- rho * .e192
    .e384 <- rho * .e241
    .e385 <- rho * .e35
    .e386 <- sigma^6
    .e387 <- (((.e270 - .e104)/.e3 - .e242) * .e8 + .e164)/.e3
    .e389 <- .e340/.e3 + .e256
    .e390 <- .e360 + .e285
    .e391 <- .e234 * .e8
    .e393 <- 1 - .e153
    .e394 <- 2 * .e40
    .e395 <- 2 * (.e86 * .e363/.e29)
    .e396 <- 2/.e23
    .e398 <- .e15 * .e119 * .e4
    .e399 <- .e371 - .e313/.e18
    .e400 <- .e84 - .e288
    .e401 <- .e317/.e18
    .e402 <- .e227 * .e8
    .e403 <- .e76 * .e15
    .e405 <- .e81 * .e8/.e3
    .e406 <- .e11 * .e215
    .e408 <- .e375 * .e77
    .e409 <- .e3^3
    .e410 <- 1/.e27
    .e411 <- .e369 - (.e395 + .e344)
    .e412 <- 2 * .e13
    .e413 <- 2 * .e20
    .e414 <- .e356 - .e3 * .e40
    .e415 <- sigma^5
    .e416 <- .e362 * .e16
    .e417 <- .e26/.e5
    .e418 <- .e245 - .e358
    .e419 <- .e105 - .e321
    .e420 <- 2 * .e149
    .e421 <- .e386 * .e11
    .e422 <- ((.e339 * .e3 - .e412)/.e32 + (.e244 + (2 * ((.e373 - 
        .e11 * .e214) * .e15) - 2 * (.e26 * .e51 * .e11 * .e5))/.e18)/.e18 - 
        .e74 * .e16 * .e8/.e3) * .e15
    .e423 <- (((2 * (.e377 - .e194) - .e193) * .e15/.e18 + .e180)/.e18 + 
        .e393 * .e3/.e32 - .e185) * .e15
    .e425 <- .e66 * .e4/sigma
    .e426 <- .e408/.e32
    .e427 <- .e352 + .e115
    .e429 <- (2 + .e33) * .e4/sigma
    .e430 <- 1/.e205
    .e431 <- 2 * .e379
    .e432 <- rho * .e253
    .e433 <- .e96 * .e16
    .e434 <- .e415 * .e11
    .e435 <- .e348 + .e406
    .e436 <- .e142 * .e16
    .e437 <- .e16 * .e215
    .e441 <- .e309 * .e15/.e18
    .e443 <- .e77/.e3 - .e163
    .e445 <- .e83/.e3 - .e93
    .e446 <- 2 * .e16
    .e447 <- .e128 - .e247
    .e448 <- .e96 - .e57 * .e15 * .e8
    .e449 <- rho * .e254
    .e450 <- .e435 * .e5
    .e451 <- (.e417 - .e405) * .e5
    .e452 <- .e76/.e27
    .e453 <- .e261 * .e8
    .e454 <- .e66 * .e207
    .e455 <- .e212 * .e8
    .e457 <- .e16 * .e83 + mu2
    .e459 <- .e419 * .e8/.e3
    .e460 <- .e376 - .e62 * .e290
    .e461 <- .e355/.e3
    .e462 <- .e413^2
    .e464 <- .e429 + rho * (.e427 + .e118)/.e3
    .e465 <- 1 + .e394
    .e466 <- .e247 - .e128
    .e467 <- rho * .e43
    .e468 <- .e38 * .e5
    .e469 <- .e98 * .e312
    .e470 <- .e418/.e23
    .e471 <- .e264/.e17
    .e472 <- .e365 * .e15
    .e473 <- .e3 * .e62
    .e474 <- .e351 - .e84
    .e475 <- .e39/.e3
    .e476 <- .e8 * .e230
    .e477 <- 2 - .e304
    .e478 <- .e372 * .e179
    .e479 <- (.e84 - .e396) * .e8
    .e481 <- .e389 * .e4/sigma
    .e482 <- .e58 * .e214
    .e486 <- .e325 + rho * (.e138 - (.e298 + .e331))
    .e487 <- .e167 + .e425
    .e488 <- .e192 * .e8
    .e489 <- .e28 + .e141
    .e490 <- .e475 - .e163
    .e491 <- 2 * (.e26 * .e85 * .e3/.e18)
    .e492 <- 2 * (.e35 * .e231 * .e11 * .e5/.e29)
    .e493 <- 4 * .e66
    .e494 <- .e398 + .e478
    .e495 <- rho * .e75
    .e496 <- rho * .e318
    .e497 <- .e467 * .e3
    .e498 <- .e323 * .e3
    .e500 <- .e421 * .e3 * .e5
    .e501 <- .e390/.e18
    .e502 <- .e325 + rho * (.e138 - (.e298 + .e492))
    .e503 <- .e391 + .e452
    .e504 <- .e454 + .e384
    .e507 <- .e330 * .e15/.e18
    .e508 <- .e269/.e91
    .e510 <- 1 - .e273
    .e511 <- 1/.e91
    .e512 <- 2 * .e66
    .e514 <- .e368 - .e416
    .e515 <- 2/.e17
    .e516 <- .e138 - (.e175 + .e331)
    .e518 <- .e15 * .e8 * .e414
    .e519 <- .e15 * .e399
    .e520 <- .e449 * .e8
    .e521 <- .e146/.e17
    .e522 <- .e304 - 1
    .e523 <- rho * (.e36 - .e446)
    .e527 <- .e497 * .e8 * .e4/.e136 - .e447/.e29
    .e528 <- .e20 * .e75
    .e530 <- .e469 * .e5/.e136
    .e532 <- .e434 * .e3 * .e5
    .e534 <- (((.e104 - .e270)/.e3 + .e242) * .e8 - .e164)/.e3 + 
        .e257
    .e537 <- .e148/.e134
    .e538 <- .e450 - .e482
    .e539 <- .e451 - .e491
    .e541 <- (.e28 + rho) * .e16
    .e543 <- (.e457 + .e523) * .e15/.e5
    .e545 <- .e26 * .e490 + .e437
    .e546 <- .e487 + .e241 * .e8
    .e548 <- (.e56 + 2 * .e30) * .e15/.e18
    .e549 <- .e403/.e18
    .e550 <- .e76/.e3
    .e551 <- .e460/.e29
    .e553 <- (.e461 - .e80) * .e8 + .e12
    .e554 <- .e464/.e17
    .e555 <- .e140 * .e11
    .e563 <- .e318 * .e15
    .e564 <- .e8 * (rho * (.e110 + .e80) - 1)
    .e565 <- .e8 * .e170
    .e566 <- .e94^2
    .e567 <- .e205^2
    .e568 <- .e12 + 2 * (.e127 * .e83 * .e5/.e32)
    .e569 <- .e267 - .e361
    .e570 <- 16 * (.e35/.e462)
    .e571 <- 2 * (.e127 * .e77 * .e5/.e32)
    .e572 <- 2 * .e367
    .e574 <- 2 * .e279 - .e493
    .e575 <- 2 * (rho * .e156/.e3)
    .e576 <- 2 * (.e323 * .e62 * .e5/.e29)
    .e577 <- .e256 - .e338
    .e580 <- rho * ((.e479 - .e164)/.e3 + .e257) * .e4/sigma
    .e581 <- rho * (.e65 - 3 * .e19)
    .e582 <- .e218 - 1
    .e583 <- .e385 * .e75
    .e585 <- sigma * .e284 * .e15
    .e586 <- .e94 * .e303
    .e587 <- .e20 * .e26
    .e588 <- .e537 + .e346
    .e589 <- .e359 + .e253/.e3
    .e590 <- ((.e57 * .e166 + .e11 * .e39/.e3) * .e5 - .e81 * 
        .e214) * .e15
    .e591 <- .e423 + .e387
    .e592 <- .e284 * .e15
    .e594 <- (.e272 + .e122) * .e5 - 2 * (.e16 * .e85 * .e3 * 
        .e15/.e18)
    .e595 <- .e471 - .e436 * .e8
    .e597 <- .e103/.e18 + .e59 * .e8/.e32
    .e601 <- .e402 + .e149
    .e602 <- (.e3 * .e367 + .e581)/.e5
    .e603 <- (.e508 - .e583 * .e8/.e136) * .e4
    .e604 <- .e453/.e3
    .e606 <- .e445/.e18 + .e3 * .e83/.e32
    .e607 <- .e95 * .e62
    .e609 <- .e448/.e3 - 2 * (.e139 * .e15 * .e5/.e18)
    .e610 <- 2 * .e244
    .e611 <- .e432 * .e8
    .e612 <- .e383 * .e4
    .e614 <- rho * .e212 * .e4
    .e615 <- rho^3
    .e617 <- sigma * .e283 * .e15
    .e619 <- .e528 * .e8/.e135
    .e621 <- .e20 * .e283 * .e15
    .e622 <- .e35 * .e3
    .e624 <- (.e543 - .e572) * .e5
    .e625 <- .e486/.e18
    .e627 <- .e387 - .e257
    .e628 <- .e489 * .e15
    .e631 <- .e26 * .e76/.e18
    .e633 <- (.e285 - .e26 * .e40) * .e11
    .e635 <- .e443 * .e15/.e18
    .e637 <- .e607 * .e5/.e29
    .e639 <- .e473 * .e5/.e134
    .e640 <- .e56 + .e496
    .e642 <- .e411 * .e5/.e29
    .e643 <- (.e518 + .e495 * .e198)/.e32
    .e644 <- .e494/.e29
    .e645 <- .e230/.e98
    .e646 <- .e288 - .e84
    .e647 <- 2 * (.e590/.e18)
    .e648 <- 2 * (.e433 * .e85 * .e3 * .e15/.e18)
    .e649 <- .e123 + 3
    .e650 <- rho * (.e422 + (.e546 + .e575)/.e3)
    .e651 <- sigma * .e75
    .e652 <- sigma * (.e346 + .e619)
    .e655 <- sigma * .e477 * .e15/.e29
    .e656 <- .e20 * .e588
    .e658 <- .e595/.e3 + .e401
    .e659 <- .e624 - .e602
    .e661 <- .e601/.e3 + .e441
    .e662 <- (.e149 + .e520)/.e3
    .e663 <- .e401 + .e264/.e23
    .e664 <- .e603 - .e551
    .e665 <- .e92/.e17
    .e667 <- .e466 * .e15/.e29
    .e668 <- .e269/.e205
    .e669 <- .e516/.e18
    .e670 <- .e4/.e468
    .e672 <- .e4/.e91 - .e652
    .e673 <- 4 * .e89
    .e674 <- .e15 * .e527
    .e675 <- .e279 - .e66
    .e676 <- .e383 * .e16
    .e677 <- .e146/.e23
    .e678 <- .e6/.e500
    .e680 <- .e20 * .e502/.e29
    .e682 <- .e147 * .e522/.e29
    .e685 <- .e622 * .e15 * .e8/.e136
    .e688 <- .e75 * .e77/.e32 + .e545/.e18
    .e690 <- .e75 * .e83/.e32
    .e694 <- (.e633 - .e433 * .e15/.e5) * .e5 + .e648
    .e696 <- ((.e565 - .e380)/.e3 + (2 * (.e11 * .e170 * .e5) + 
        .e431) * .e15/.e18) * .e15/.e18
    .e697 <- (.e92 * .e16 + .e157)/.e3
    .e699 <- .e317 * .e11 * .e5
    .e700 <- .e628 * .e8
    .e701 <- .e227 * .e16
    .e702 <- .e234 * .e16
    .e703 <- .e403 + .e96
    .e704 <- .e640/.e17
    .e706 <- .e66 * .e16/.e3
    .e707 <- .e66 * .e40
    .e709 <- .e455 - .e430
    .e711 <- .e142 * .e4/sigma
    .e712 <- .e259 * .e8
    .e714 <- (.e465 - .e102)/.e17
    .e715 <- .e267/.e17
    .e716 <- .e56 + rho * (.e568 + .e14)
    .e719 <- (.e354 + 2 * (.e224 * .e11 * .e5)) * .e15/.e18
    .e720 <- (rho * (.e28 - .e274)/.e3 - 2)/.e17
    .e721 <- .e645 - .e530
    .e722 <- .e98^2
    .e723 <- .e465 - .e53
    .e724 <- .e430 - .e455
    .e725 <- 2 * (.e594 * .e15/.e18)
    .e726 <- 2 * (.e555 * .e5)
    .e727 <- 2 * (.e40/.e17)
    .e729 <- 2 * (.e498 * .e15/.e29)
    .e730 <- 2 + rho * .e355/.e3
    .e731 <- .e39 - .e541
    .e734 <- rho * .e534 * .e4/sigma
    .e736 <- .e96 * .e76/.e3
    .e737 <- rho * (.e34 + .e14)
    .e739 <- .e586 * .e5/.e135
    .e740 <- sigma^7
    .e741 <- sigma^8
    .e742 <- .e400 * .e16
    .e744 <- .e597 * .e15
    .e746 <- .e606 * .e15
    .e747 <- .e240/.e17
    .e748 <- .e240/.e27
    .e754 <- .e273/.e23
    .e755 <- .e76 * .e4
    .e756 <- .e261 * .e16
    .e758 <- .e254 * .e16
    .e759 <- .e609/.e18
    .e760 <- .e721 * .e8
    .e761 <- .e92/.e23
    .e762 <- .e95 * .e230
    .e763 <- .e459 + .e46
    .e764 <- .e724/.e3
    .e765 <- (2 * .e538 - 2 * (.e103 * .e26 * .e11 * .e5))/.e18
    .e767 <- (2 * .e539 - 2 * (.e103 * .e11 * .e5)) * .e15/.e18
    .e768 <- .e563 - .e25
    .e769 <- (.e611 - .e149)/.e3
    .e770 <- (.e521 - .e232)/.e3
    .e771 <- .e4/.e99
    .e772 <- 1 - .e570
    .e774 <- 2 * .e54
    .e775 <- 2 * .e426
    .e776 <- 2 * .e63
    .e777 <- 2 * sigma
    .e778 <- rho * (.e481 - .e201)
    .e781 <- rho * (.e423 + ((.e396 - .e84) * .e8 + .e164)/.e3) * 
        .e4/sigma
    .e782 <- rho * .e627
    .e783 <- rho * (.e604 + .e426)
    .e784 <- .e612/.e98
    .e785 <- .e614/.e98
    .e786 <- .e280 * .e15
    .e788 <- rho * (.e511 - .e685) * .e4
    .e789 <- .e6/.e532
    .e790 <- .e94 * .e4
    .e791 <- .e587/.e29
    .e793 <- ((.e553 * .e16 + .e157)/.e3 + .e765) * .e15/.e18
    .e794 <- ((.e553/.e3 + .e767)/.e18 + .e790/.e353) * .e15
    .e795 <- .e659 * .e8
    .e797 <- .e661 * .e4/sigma
    .e800 <- .e690 + .e635
    .e801 <- .e359 + .e706
    .e805 <- ((.e11 * (2 - 2 * .e224) * .e5 - .e354) * .e15/.e18 + 
        (mu2 + .e737 - .e564)/.e3) * .e15/.e18
    .e806 <- .e703/.e18
    .e807 <- (.e239 - 1)/.e3
    .e809 <- ((.e564 - .e30)/.e3 + .e719) * .e15/.e18
    .e812 <- (.e707 + .e459 + .e164 + rho * (.e183 - (.e742 + 
        .e174 + .e160))) * .e8 + .e321
    .e813 <- .e592 + .e709/.e94
    .e814 <- .e68 * .e62
    .e816 <- .e264 * .e15/.e18
    .e817 <- .e26 * (2 * (.e498 * .e62/.e29) + .e45)
    .e819 <- .e716 * .e3/.e32
    .e820 <- .e338 + .e441
    .e823 <- (.e8 * (.e279 - (.e776 + .e164)) + .e515)/.e3 + 
        .e580 - .e242
    .e824 <- .e488 - .e410
    .e825 <- .e212 * .e207
    .e827 <- .e365 * .e8/.e3
    .e831 <- .e287 * .e3 * .e5/.e353
    .e832 <- .e762 + rho * .e664
    .e834 <- .e3 * (.e30 * .e198 - .e56)/.e32
    .e835 <- .e330 * .e8
    .e836 <- .e507 + .e677
    .e837 <- .e763 + .e384
    .e839 <- .e574 * .e8/.e3
    .e840 <- .e376/.e5
    .e842 <- .e8 * .e399/.e3
    .e843 <- .e8 * .e249
    .e844 <- .e476 - .e668
    .e845 <- .e476/.e3
    .e846 <- .e771 - .e656
    .e847 <- .e12 + .e306
    .e848 <- .e410 - .e676
    .e850 <- 2 * (.e390 * .e11 * .e5)
    .e851 <- 2 * .e694
    .e852 <- 2 * .e160
    .e853 <- 2 * .e125
    .e855 <- .e512 + 2 * (rho * .e474 * .e4/sigma)
    .e856 <- 2 * .e142
    .e858 <- 2 * (rho * .e504 * .e4/sigma) - 2 * .e227
    .e859 <- 2 * (rho * .e299/.e3)
    .e860 <- 2 * (.e116 * .e62 * .e5/.e29)
    .e861 <- 3 + .e199
    .e862 <- 4 * .e12
    .e863 <- 4 * .e372
    .e865 <- .e15 * .e582/.e18
    .e866 <- .e778 - .e129
    .e867 <- rho * .e591
    .e870 <- rho * (.e696 + (.e714 + .e577 * .e8 - .e548)/.e3) * 
        .e4/sigma
    .e872 <- rho * (.e760 - .e678) * .e4
    .e873 <- rho * (.e725 + .e610)
    .e876 <- rho * .e731 * .e8/.e3
    .e878 <- rho * .e77 * .e8
    .e879 <- rho * sigma
    .e880 <- .e651 * .e4
    .e881 <- .e24 * .e411
    .e882 <- .e35 * .e75
    .e883 <- .e741 * .e11
    .e887 <- .e794 + (.e70 - (.e700 + .e22 * (2 * (.e75 * .e3 * 
        .e15 * .e8 * .e5/.e32) + .e137)))/.e32
    .e891 <- (.e503/.e94 + .e655) * .e4
    .e892 <- ((.e704 - .e232)/.e3 + (3 + rho * ((.e141 - (.e774 + 
        .e260))/.e3 - .e673)) * .e15/.e18) * .e8
    .e894 <- (.e26 * .e39/.e3 + .e437) * .e11
    .e895 <- .e253 * .e8
    .e896 <- .e662 + .e519
    .e897 <- .e240/.e23
    .e898 <- .e747 + .e701 * .e8
    .e899 <- .e814 + (.e840 - .e860) * .e11 * .e3
    .e900 <- .e470 + .e783
    .e901 <- .e754 + .e669
    .e903 <- .e26 * .e92/.e18
    .e904 <- .e139 * .e214
    .e909 <- .e26 * .e57
    .e910 <- .e58 * .e290
    .e912 <- (.e11 * (4 * sigma + .e196) - .e45)/.e5 - .e576
    .e913 <- .e702 * .e8
    .e914 <- .e254 * .e8
    .e915 <- .e472 + .e769
    .e916 <- .e510/.e17
    .e918 <- .e11 * (.e863 + .e398) - .e817
    .e920 <- .e11 * .e448/.e3
    .e921 <- .e569/.e17
    .e922 <- .e507 + .e770
    .e923 <- .e764 - .e585
    .e924 <- (2 * .e266 - .e275) * .e5
    .e925 <- .e563/.e18
    .e926 <- .e514/.e17
    .e927 <- .e368 * .e8
    .e929 <- .e62 * .e303 + .e881
    .e930 <- .e39 * .e4
    .e931 <- .e845 + .e382
    .e933 <- 1/.e99
    .e934 <- 2 * (.e538/.e18)
    .e935 <- 2 * (.e539 * .e15/.e18)
    .e936 <- 2 * .e178
    .e937 <- 2 * .e697
    .e938 <- 2 * (.e100 * .e85 * .e3 * .e15/.e18)
    .e939 <- 2 * .e192
    .e940 <- 2 * .e129
    .e941 <- .e34 + rho * (2 * .e77 + .e223 + .e248)/.e3
    .e944 <- .e223 + mu2 + .e37 - .e69
    .e945 <- 2 * (.e385 * .e26 * .e85 * .e3/.e29)
    .e946 <- 2 * .e173
    .e947 <- 2 * (.e111 * .e85 * .e3/.e29)
    .e948 <- .e576 - (.e38 * .e861 - .e45)/.e5
    .e949 <- 2 * (.e434 * .e409 * .e303/.e135)
    .e950 <- 2 * (.e421 * .e113 * .e15 * .e8 * .e5/.e135)
    .e951 <- 2/.e27
    .e952 <- 4 - .e570
    .e953 <- rho * ((.e795 + .e75 * .e339)/.e32 + .e793 + .e880/.e353)
    .e954 <- rho * (.e481 - 2 * .e201)
    .e955 <- .e867 * .e15
    .e957 <- .e782 * .e4/sigma
    .e958 <- .e872 - .e642
    .e959 <- .e784 - .e639
    .e960 <- .e785 - .e739
    .e961 <- rho * .e672
    .e963 <- .e786 * .e4/.e353
    .e964 <- .e6/.e91
    .e965 <- .e43 * .e3
    .e968 <- .e589 * .e16
    .e970 <- (.e894 + .e909 * .e15/.e5) * .e5 - .e904
    .e971 <- .e400 * .e207
    .e973 <- ((.e142 * .e40 + rho * (.e658 - .e16 * .e577)) * 
        .e8 + rho * (.e816 + .e436) - (.e56 + rho * ((3 + .e33) * 
        .e4/sigma + rho * (.e352 + .e380 + 3 * mu2 + rho * (.e31 + 
        .e862))/.e3))/.e17)/.e3 + (.e15 * (rho * ((.e712 + .e245 - 
        rho * .e649)/.e3 + (2 * .e699 - .e851)/.e18) - (.e106 * 
        .e119 + 1)) + .e736)/.e18
    .e975 <- .e898/.e3 + .e625
    .e984 <- .e450 - .e910
    .e986 <- (.e348 + .e920) * .e5 - .e938
    .e987 <- .e625 + .e897
    .e988 <- .e451 - .e947
    .e990 <- .e503 * .e4
    .e997 <- ((.e309 * .e8 + .e380)/.e3 + (2 * (.e11 * .e309 * 
        .e5) - .e354) * .e15/.e18) * .e15/.e18
    .e1000 <- .e809 + .e643
    .e1002 <- .e283 * .e8/.e3
    .e1003 <- .e284 * .e8
    .e1004 <- .e592 - .e511
    .e1007 <- (.e244 + .e35 * (2 * ((.e373 - .e11 * .e290) * 
        .e15) - 2 * (.e26 * .e140 * .e11 * .e5))/.e29)/.e29
    .e1010 <- .e240/.e97 + .e680
    .e1011 <- .e748 + .e913
    .e1013 <- (.e16 * (.e50 - .e565/.e3) + .e165 + rho * (.e647 + 
        2 * .e807) - 2 * (.e58 * .e170 * .e5/.e18)) * .e15/.e18
    .e1017 <- .e265 * .e11
    .e1019 <- (.e457 + .e37) * .e15/.e5
    .e1025 <- .e315 + .e945
    .e1029 <- .e227 * .e4/sigma
    .e1034 <- (.e3 * (rho * (2 * (.e250 * .e5/.e32) - .e295) - 
        1) + .e356) * .e15/.e32
    .e1035 <- .e76 * .e66
    .e1038 <- .e338 + .e519
    .e1040 <- .e76/.e99 + .e655
    .e1041 <- .e768/.e18
    .e1043 <- (.e461 - .e158) * .e8 + .e12
    .e1049 <- .e844/.e94 + .e674
    .e1055 <- .e914 + .e665
    .e1057 <- .e254 * .e4/sigma
    .e1061 <- .e759 + .e426
    .e1063 <- (.e670 + .e262)/.e3 + .e16 * .e646
    .e1068 <- .e665 - .e895
    .e1070 <- .e712/.e3
    .e1071 <- .e16 * .e207
    .e1072 <- .e245 + rho * (.e571 - .e176)
    .e1073 <- .e245 + .e876
    .e1075 <- .e510/.e23
    .e1076 <- .e916 - .e758 * .e8
    .e1079 <- .e95 * .e119
    .e1080 <- .e723/.e17
    .e1082 <- .e723/.e27 - .e220
    .e1083 <- .e921 - .e453
    .e1086 <- .e473 * .e8 * .e5/.e134
    .e1088 <- .e667 - (.e789 + .e586 * .e8 * .e5/.e135)
    .e1089 <- .e667 + .e8 * .e960
    .e1091 <- (1 + .e1 * .e649/.e3) * .e8 + rho * (.e92 * .e40 + 
        .e16 * .e267 + .e12 + .e203)
    .e1092 <- (.e180 + .e35 * (2 * (.e377 - .e729) - .e726) * 
        .e15/.e29)/.e29
    .e1094 <- .e848 * .e8/.e3
    .e1096 <- (2 - .e199)/.e205 + .e191
    .e1098 <- (.e941 - (.e121 + .e944 * .e8/.e3) * .e16) * .e15/.e3
    .e1099 <- .e514/.e23
    .e1100 <- .e926 - .e756 * .e8
    .e1101 <- .e839 + 2 * (.e580 - .e178)
    .e1103 <- (.e859 - ((.e474 * .e16 + .e183 - .e216) * .e8 + 
        .e167 + .e425))/.e3 - .e422
    .e1104 <- .e413^4
    .e1105 <- .e730/.e17
    .e1106 <- .e730/.e27
    .e1108 <- .e8 * (rho * (.e110 + .e158) - 1) - .e30
    .e1109 <- .e8 * .e959
    .e1110 <- (.e20 * (2 * ((.e11 * (.e196 + sigma * (4 - .e128)) - 
        .e45) * .e15 * .e5) - 2 * (.e555 * .e3 * .e62))/.e29 - 
        .e6/.e98)/.e29
    .e1111 <- (.e43 * (2 * (.e237 + .e729) + .e726) * .e15/.e29 - 
        1/.e94)/.e29
    .e1112 <- .e532^2
    .e1113 <- .e500^2
    .e1114 <- .e568 + .e31
    .e1115 <- 1 - (.e273 + 2 * (.e250 * .e83 * .e5/.e32))
    .e1116 <- .e267 - (.e361 + .e571)
    .e1117 <- .e410 - .e488
    .e1119 <- 16 * .e772 + 16 * .e952
    .e1120 <- 2 * (.e912 * .e3)
    .e1121 <- 2 * .e400
    .e1122 <- 2 * (.e699/.e18)
    .e1123 <- 2 * .e899
    .e1124 <- 2 * (.e26 * (.e924 + .e319))
    .e1125 <- 2 * .e21
    .e1127 <- 2 * (.e11 * .e249 * .e5)
    .e1129 <- 2 * (.e3 * .e948)
    .e1130 <- 2 * .e929
    .e1131 <- .e939 + 2 * (rho * .e923 * .e4)
    .e1132 <- .e368 - (.e416 + 2 * (.e250 * .e77 * .e5/.e32))
    .e1135 <- 2 * (rho * .e813 * .e4) - 4 * .e212
    .e1137 <- 2 * (rho * (.e825 + .e961) * .e4) - 2 * .e234
    .e1138 <- 2 * (rho * .e340/.e3)
    .e1139 <- 2 * (rho/.e23)
    .e1140 <- 2 * (.e740 * .e11 * .e409 * .e312/.e136)
    .e1142 <- 2/.e99
    .e1143 <- 4 * .e400
    .e1145 <- .e15 * (.e873 - .e51 * .e170)/.e18
    .e1146 <- rho * ((.e892 + .e720 + .e496 * .e15/.e18)/.e3 + 
        ((.e835 + .e380)/.e3 + (2 * (.e11 * .e330 * .e5) - .e431) * 
            .e15/.e18) * .e15/.e18)
    .e1147 <- .e734 - .e178
    .e1148 <- rho * ((.e819 + .e842 - ((2 * (.e11 * .e313 * .e5) + 
        .e431) * .e15/.e18 - .e550)/.e18) * .e15 + (.e549 + ((.e420 + 
        .e520)/.e3 + 2 * .e519) * .e8 - (.e727 + .e449))/.e3)
    .e1149 <- rho * ((((.e354 + 2 * (.e11 * .e294 * .e5)) * .e15/.e18 - 
        .e550)/.e18 + .e827 + .e834) * .e15 + (((.e611 - .e420)/.e3 + 
        2 * .e472) * .e8 + .e727 - (.e549 + .e432))/.e3)
    .e1151 <- rho * (.e824/.e94 + .e617) * .e4
    .e1153 <- rho * ((.e3 * (2 * .e245 + rho * (.e571 - (4 + 
        4 * .e33))) + .e878) * .e15/.e32 + ((2 * .e470 + rho * 
        (.e604 + .e775) - 2 * .e501) * .e8 - rho * .e261)/.e3) - 
        (((rho * ((.e28 - .e141)/.e3 + .e93) - 2) * .e119 + rho * 
            ((.e850 + 2 * (.e26 * .e379))/.e18 - .e76 * .e16/.e3)) * 
            .e15/.e18 + .e40 * (.e208 + 4 * .e30)/.e23)
    .e1154 <- rho * ((.e8 * (.e954 - .e940) + .e420)/.e3 + .e870 - 
        .e15 * (rho * (4 * .e110 + .e673) - 4)/.e18)
    .e1155 <- .e96 * .e92
    .e1156 <- rho * .e846
    .e1158 <- .e146/.e27
    .e1159 <- .e788 - .e931
    .e1160 <- rho * (.e853 + 3 * .e122)
    .e1163 <- .e1 * .e574 * .e15 - 4
    .e1164 <- .e1 * (2 + 2 * .e176)
    .e1165 <- .e682 - .e146/.e97
    .e1166 <- .e43 * .e26
    .e1167 <- .e882 * .e11
    .e1168 <- .e883 * .e409
    .e1169 <- -(2 * ((.e95 * .e675 + rho * (((.e104 - .e515)/.e3 + 
        .e84 + .e650) * .e4/sigma + .e174 + .e160)) * .e15/.e20))
    .e1170 <- -(2 * (rho * (.e84 + .e781 - (.e8 * .e675 + .e105)/.e3) * 
        .e15/.e20))
    .e1171 <- -(2 * (rho * (.e84 + .e957 - .e351) * .e15/.e20))
    .e1172 <- -(2 * (.e42 * (rho * .e823 * .e4/sigma - 2 * .e675)/.e43))
    .e1173 <- -(2 * (.e1 * .e591 * .e15/.e20))
    .e1174 <- -(2 * (.e1 * .e823 * .e15/.e43))
    .e1178 <- .e973 * .e4/sigma
    .e1180 <- (.e663 * .e8 + .e631 + .e711 + .e1138 - .e554)/.e3 + 
        .e1013
    .e1185 <- .e801 * .e8
    .e1190 <- (.e178 + .e650) * .e4/sigma
    .e1191 <- ((.e92 * .e419 + .e670 - (.e167 + (.e66 + .e384) * 
        .e16 + .e262 * .e39/.e3) * .e8)/.e3 - ((.e359 + .e837 * 
        .e16/.e3) * .e8 + .e852 + .e953)) * .e15
    .e1193 <- .e1011/.e3 + .e680
    .e1196 <- .e75 * .e16 * .e8/.e32
    .e1198 <- (.e541 - .e39) * .e8/.e3
    .e1199 <- (.e348 + 2 * .e406) * .e5
    .e1201 <- .e486 * .e11 * .e5
    .e1203 <- .e253 * .e4/sigma
    .e1204 <- (.e417 - 2 * .e405) * .e5
    .e1206 <- (.e402 + .e521)/.e3
    .e1209 <- (.e391 + .e1158)/.e3
    .e1211 <- (.e836 * .e8 + .e548 + .e720)/.e3 + .e1145
    .e1212 <- .e504 * .e8
    .e1213 <- .e1003/.e3
    .e1215 <- .e1055/.e3 + .e746
    .e1216 <- .e746 + .e761
    .e1217 <- .e1068/.e3
    .e1218 <- .e240 * .e26
    .e1220 <- .e240 * .e15/.e18
    .e1221 <- .e900 - .e501
    .e1223 <- .e901 - .e968 * .e8
    .e1225 <- .e1076/.e3 - .e800
    .e1226 <- .e832 * .e8
    .e1227 <- .e234 * .e4
    .e1228 <- .e1083/.e3
    .e1230 <- (.e3 * .e1132 - .e489 * .e77) * .e15/.e32
    .e1233 <- ((.e105 - ((.e459 + .e164 + rho * (.e1063 - .e216)) * 
        .e8 + .e321))/.e3 - (.e51 * .e15/.e18 + rho * .e887)) * 
        .e15
    .e1234 <- (.e351 + rho * .e1103 - .e84) * .e15
    .e1235 <- .e837/sigma
    .e1237 <- .e1100/.e3 - .e688
    .e1239 <- ((2 * .e218 - 3) * .e15/.e18 - (.e402 + .e704)/.e3) * 
        .e8 + .e1105
    .e1240 <- .e443 * .e8
    .e1241 <- .e192 * .e4
    .e1242 <- .e847/.e17
    .e1245 <- .e16 * .e393 * .e3 + .e13
    .e1247 <- .e510 * .e15/.e18
    .e1248 <- .e1075 - .e800
    .e1253 <- .e1116 * .e3 - .e77 * .e8
    .e1254 <- .e569 * .e15
    .e1256 <- .e569/.e23 - .e1061
    .e1258 <- .e267 * .e8 + mu2
    .e1260 <- .e1117/.e3
    .e1262 <- .e855 * .e8/.e3
    .e1263 <- .e1131 * .e8
    .e1265 <- .e514 * .e15/.e18
    .e1266 <- .e1099 - .e688
    .e1267 <- .e1135 * .e8
    .e1269 <- .e858 * .e8/.e3
    .e1270 <- .e1137 * .e8
    .e1271 <- .e269/.e94
    .e1272 <- .e930/sigma
    .e1273 <- .e1109 + .e382
    .e1274 <- .e8 * .e527
    .e1276 <- .e8^2/.e3
    .e1278 <- .e1153 * .e4/sigma
    .e1279 <- .e12 - 2 * (.e1167 * .e3 * .e8 * .e5/.e134)
    .e1280 <- .e2 - 3 * .e1
    .e1281 <- 16 * (.e772/.e462)
    .e1282 <- 16 * (.e43 * .e1119/.e1104)
    .e1284 <- 2 * (.e986/.e18)
    .e1286 <- 2 * .e452
    .e1287 <- 2 * (.e11 * .e494 * .e5)
    .e1288 <- 2 * .e474
    .e1289 <- 2 * .e256
    .e1290 <- 2 * (.e955/sigma)
    .e1291 <- 2 * .e781
    .e1292 <- 2 * .e521
    .e1293 <- 2 * .e271
    .e1294 <- 2 * (.e42 * .e958)
    .e1295 <- .e950 - 1
    .e1300 <- .e777 + .e196
    .e1301 <- .e1142 - .e585
    .e1302 <- 4 * .e72
    .e1303 <- .e865 - .e677
    .e1304 <- .e382 - (.e1086 + .e964)
    .e1306 <- .e1146 * .e4/sigma
    .e1308 <- rho * (((.e812 - .e105)/.e3 + .e805 - .e643)/.e20 - 
        .e831) * .e4
    .e1309 <- .e734 - .e936
    .e1311 <- rho * (.e797 + .e167) - .e211
    .e1313 <- rho * (.e891 + .e791) - .e221
    .e1315 <- rho * (.e990/.e98 - .e644) - .e637
    .e1317 <- .e1148 * .e4/sigma
    .e1319 <- .e1149 * .e4/sigma
    .e1320 <- .e279 - .e512
    .e1321 <- .e1151 - .e192
    .e1325 <- rho * .e389 * .e4/sigma - .e856
    .e1326 <- rho * (.e755/.e91 - .e644)
    .e1327 <- rho * (.e934 + .e937)
    .e1328 <- rho * (.e935 + .e610)
    .e1329 <- rho * .e230
    .e1331 <- .e385 * .e15 * .e4
    .e1332 <- .e1 * .e20
    .e1333 <- .e617 - .e933
    .e1334 <- .e38 * .e466
    .e1335 <- .e24 * .e409
    .e1338 <- .e469 * .e8 * .e5/.e136
    .e1339 <- y - (2 * (.e740 * .e75 * .e11 * .e3 * .e8 * .e5/.e135) + 
        mu1)
    .e1340 <- y - (2 * (sigma^9 * .e75 * .e11 * .e3 * .e8 * .e5/.e136) + 
        mu1)
    .e1341 <- -(2 * ((.e1178 - .e658) * .e15/.e20))
    .e1342 <- -(2 * ((.e481 + (.e297 - .e21)/.e18 + rho * (.e1180 * 
        .e4/sigma - .e635) - (.e471 + .e28 * .e866)/.e3) * .e15/.e20))
    .e1343 <- -(2 * ((.e501 + .e1278 - .e900) * .e15/.e20))
    .e1344 <- -(2 * (.e1191/sigma))
    .e1345 <- -(2 * (.e1233/sigma))
    .e1346 <- -(2 * (.e1234/sigma))
    .e1347 <- -(2 * (.e15 * (.e15 * .e106/.e18 + rho * (.e1211 * 
        .e4/sigma - (.e445 * .e15/.e18 + .e761)) - .e8 * .e866/.e3)/.e20))
    .e1348 <- -(2 * (.e15 * (.e1306 - .e922)/.e20))
    .e1349 <- -(2 * (.e15 * (.e1308 - .e1235)/sigma))
    .e1350 <- -(2 * (.e15 * (.e1317 - .e896)/.e20))
    .e1351 <- -(2 * (.e15 * (.e1154 * .e4/sigma - 2 * .e866)/.e43))
    .e1352 <- -(2 * (rho * (((.e812 - .e515)/.e3 + .e805 - .e643)/sigma + 
        .e963) * .e15/sigma))
    .e1353 <- -(2 * (rho * (.e178 * .e207 + .e650) * .e15/.e20))
    .e1354 <- -(2 * (.e1154 * .e15/.e43))
    .e1355 <- -(2 * (.e1 * .e627 * .e15/.e20))
    .e1356 <- -(2 * (.e314 * .e1309/.e43))
    .e1357 <- -(2 * (.e615 * .e534 * .e15/.e43))
    .e1358 <- (((((.e149 + (.e105 - .e432) * .e8)/.e3 - .e472) * 
        .e16 + .e167 + rho * (.e901 - .e1196)) * .e8 + .e631 + 
        rho * (.e1203 + 2 * .e769) - .e554)/.e3 + (((.e222 * 
        .e3 - 2 * (.e1 * .e16)) * .e8 + rho * .e339 * .e16 * 
        .e3)/.e32 + (.e58 * (2 - 2 * .e294) * .e5/.e18 + .e72 + 
        rho * (.e697 + .e934))/.e18 - .e365 * .e16 * .e8/.e3) * 
        .e15) * .e15
    .e1359 <- ((.e589 * .e40 + rho * .e1223/.e3 - ((.e472 - .e338)/.e3 + 
        .e1034) * .e16) * .e8 + ((.e1073/.e3 + (2 * (.e19 * .e516 * 
        .e5) - .e1124)/.e18)/.e18 + rho * .e273 * .e3/.e32) * 
        .e15 + rho * .e589 * .e16 - .e464/.e23) * .e15
    .e1360 <- (((.e501 + rho * (.e1228 - (.e759 + .e775)) - .e470) * 
        .e8 + (.e341 * .e40 + .e16 * (1 + 3 * .e33))/.e17 + rho * 
        ((.e1254 - .e100)/.e18 + .e235))/.e3 - ((.e1072 * .e3 + 
        .e878)/.e32 + ((2 * (rho * .e609 * .e11 * .e5) - 2 * 
        ((.e633 + .e96 * .e172 * .e16) * .e5 + .e648))/.e18 + 
        (rho * .e119 - (.e8 * (.e53 - .e394) - .e30) * .e16)/.e3)/.e18) * 
        .e15) * .e15
    .e1364 <- (((.e362 * .e40 + (.e429 + rho * (2 * (.e352 + 
        mu2 + .e37) + mu2 + .e7)/.e3) * .e16 - (2 + .e1164/.e3))/.e17 + 
        (.e261 * .e40 + rho * .e1237 - (.e470 + .e786 * .e77/.e32 - 
            .e501) * .e16) * .e8 + rho * (.e756 + .e1265))/.e3 - 
        ((.e75 * .e1072 + .e130 * .e414)/.e32 + ((.e1073 * .e16/.e3 + 
            .e490 * .e119 + 2 * (rho * .e545 * .e11 * .e5/.e18)) * 
            .e15 - (.e26 * (2 * (.e694/.e18) + .e358/.e3) + .e40 * 
            .e215))/.e18)) * .e15
    .e1367 <- (((.e402 + .e420)/.e3 + 2 * .e441) * .e8 + .e806 - 
        .e1080)/.e3 + .e997
    .e1368 <- (((.e1038 + rho * .e1215) * .e8 + rho * (.e925 - 
        .e225) - .e715)/.e3 + ((.e550 + (2 * (rho * .e445 * .e11 * 
        .e5) - .e431) * .e15/.e18)/.e18 + .e819 + rho * .e606 * 
        .e8/.e3) * .e15) * .e15
    .e1369 <- (((.e338 + rho * (.e1217 - .e744) - .e472) * .e8 + 
        rho * (.e1041 + .e246) - .e715)/.e3 - ((.e719 - .e550)/.e18 + 
        .e834 + rho * .e597 * .e8/.e3) * .e15) * .e15
    .e1370 <- ((((.e226 * .e39 - .e1164)/.e3 + (.e12 + 2 * (.e40 * 
        .e4/sigma + .e358 * .e8/.e3) + .e203) * .e16 - 2)/.e17 + 
        (.e235 - ((.e1071 * .e15 - (.e360 + .e25))/.e18 + .e470 + 
            .e783) * .e16) * .e8)/.e3 + (.e408 - .e75 * .e418)/.e32 + 
        rho * ((.e1266 * .e8 + .e261 * .e92)/.e3 + .e1230) - 
        ((.e240 * .e16 + .e77 * .e119) * .e15/.e3 + .e26 * (.e11 * 
            (.e1125 - 2 * .e390) * .e5/.e18 + rho * (.e647 - 
            .e176/.e3)) + .e215 * .e334)/.e18) * .e15
    .e1373 <- ((.e1043 * .e16 + .e157)/.e3 + .e35 * (2 * .e984 - 
        2 * (.e265 * .e26 * .e11 * .e5))/.e29) * .e15/.e29
    .e1374 <- (((.e254 * .e40 + rho * .e1225 - .e1038 * .e16) * 
        .e8 + .e554 + rho * (.e758 + .e1247))/.e3 - ((.e75 * 
        .e716 + .e15 * .e83 * .e414)/.e32 + ((.e245 + rho * (.e1240 - 
        .e176))/.e3 + (2 * (rho * .e443 * .e11 * .e15 * .e5) - 
        .e851)/.e18) * .e15/.e18)) * .e15
    .e1376 <- (((.e232 - .e420)/.e3 + .e1289) * .e8 + .e714 - 
        .e548)/.e3 + .e696
    .e1377 <- .e1185 + .e160
    .e1378 <- (((.e16 * (.e394 + .e33) + rho * .e226/.e3)/.e17 - 
        .e1221 * .e8)/.e3 + rho * ((.e1256 * .e8 + .e235 - .e252)/.e3 + 
        .e1253 * .e15/.e32) - (((.e16 * (.e146 - .e8 * .e334) + 
        .e83 * .e119)/.e3 + (2 * (rho * .e594 * .e26) - .e850)/.e18)/.e18 + 
        .e418 * .e3/.e32) * .e15) * .e15
    .e1382 <- (.e1198 - (.e361 + 1))/.e3
    .e1383 <- ((.e543 - (2 * (.e75^2 * .e11 * .e3/.e32) + .e572)) * 
        .e5 - .e602)/.e32
    .e1384 <- (.e1204 - .e491) * .e26
    .e1386 <- (((.e1280 * .e83 + rho * .e1115 * .e3)/.e32 - ((.e110 + 
        .e11 * (.e1125 - 2 * (.e26 * .e313)) * .e5/.e18 + rho * 
        (.e807 + .e647))/.e18 + .e28 * .e399/.e3)) * .e15 + (.e554 + 
        (.e225 + rho * .e1248 - (.e662 + ((2 - .e293)/.e18 + 
            .e371) * .e15) * .e16) * .e8 + rho * (.e1057 + 2 * 
        .e662) - .e631)/.e3) * .e15
    .e1389 <- ((.e549 + rho * (.e167 + (.e318/.e23 - .e744) * 
        .e8) - (.e915 * .e8 + .e715))/.e3 + ((rho * (.e244 + 
        .e935) - 2 * (.e47 * .e294 * .e5/.e18))/.e18 + rho * 
        .e1245/.e32 - .e827) * .e15) * .e15
    .e1392 <- ((((2 * .e112 + rho * (.e927 + .e774 + 2 * (.e28 + 
        rho * (.e273 - 2)) + 2 * .e1272 + 2 * (rho * .e228/.e3)))/.e3 + 
        2 * .e121 - 2)/.e17 + 2 * (rho * .e595/.e3) - ((.e658 * 
        .e8 + .e711) * .e16 + (.e663 * .e16 + .e142 * .e39/.e3) * 
        .e8 + .e264 * .e26/.e18))/.e3 + ((.e251 + rho * ((.e114 + 
        rho * (2 * .e243 + .e223 + .e248)/.e3 - .e416)/.e3 - 
        2 * (.e970/.e18)) - (.e1070 + .e80) * .e16) * .e15 - 
        ((.e264/.e3 + .e1122) * .e26 + .e215 * .e106))/.e18) * 
        .e15
    .e1393 <- (.e1043/.e3 + .e35 * (2 * .e988 - 2 * (.e1017 * 
        .e5)) * .e15/.e29)/.e29
    .e1394 <- (((.e843 - 2 * .e6)/sigma - mu2)/.e3 + .e43 * (.e1127 - 
        .e1129) * .e15/.e29)/.e29
    .e1395 <- (.e178 + rho * (.e422 + (.e487 + (.e183 - (.e178 * 
        .e16 + .e174 + .e160)) * .e8 + .e575)/.e3)) * .e15
    .e1397 <- .e178 * (3 - .e53) + .e650
    .e1398 <- .e1004 * .e8
    .e1401 <- (((.e12 + rho * (.e649 * .e8 + mu2 + rho * (2 * 
        .e92 + .e31 + .e114))/.e3)/.e17 - (((.e471 + .e142 * 
        .e83)/.e3 + .e401 + .e836 * .e16) * .e8 + .e816 + .e436))/.e3 + 
        (.e15 * (rho * ((.e176 - .e361)/.e3 - .e1284) - (.e1070 + 
            .e1122)) - (.e1155 - .e21 * .e8 * .e106)/.e3)/.e18) * 
        .e15
    .e1404 <- ((.e1242 + rho * (2 * .e770 - .e903) - ((.e663 + 
        .e922 * .e16) * .e8 + .e711))/.e3 + (.e307 + rho * ((2 + 
        2 * (rho * .e292/.e3) - 2 * .e273)/.e3 - .e647) - ((.e835/.e3 + 
        .e50) * .e16 + 2 * (.e58 * .e330 * .e5/.e18))) * .e15/.e18) * 
        .e15
    .e1411 <- .e545 * .e11
    .e1412 <- .e1199 - .e482
    .e1414 <- (.e58 * (2 + 2 * .e309) * .e5/.e18 + .e16 * (4 - 
        .e218) * .e8/.e3 + .e1327) * .e15/.e18
    .e1416 <- ((.e1079 + rho * ((.e16 * (.e208 + mu2 + .e737) + 
        .e876)/.e3 + (2 * .e1201 - .e1124)/.e18) - 1) * .e15 + 
        .e736)/.e18 + (.e8 * (rho * (.e975 + .e820 * .e16) - 
        .e227 * .e40) + rho * (.e1220 - .e701) - .e1091/.e17)/.e3
    .e1419 <- (.e1227 + .e464/.e468)/sigma - (.e1010 * .e8 + 
        .e587 * .e76/.e29)
    .e1420 <- ((.e715 + rho * (.e225 - (.e1216 * .e8 + .e201)) - 
        (.e896 * .e8 + .e549))/.e3 - (.e842 + (rho * (.e244 + 
        .e725) - 2 * (.e47 * .e313 * .e5/.e18))/.e18 + rho * 
        .e1114 * .e3/.e32) * .e15) * .e15
    .e1423 <- .e1040 * .e8 + sigma * .e76 * .e15/.e29 - .e1082/sigma
    .e1427 <- ((.e930 - .e28 * .e249)/.e3 - .e35 * (.e1123 + 
        2 * (.e58 * .e249 * .e5))/.e29) * .e15/.e29
    .e1429 <- (.e1108 * .e4/.e98 + .e20 * (2 * (.e918 * .e5) - 
        2 * (.e1017 * .e3 * .e62))/.e29) * .e15/.e29
    .e1430 <- (.e1108/.e94 + .e43 * (2 * .e1025 + 2 * (rho * 
        .e265 * .e11 * .e5)) * .e15/.e29)/.e29
    .e1431 <- .e1212 + .e1035
    .e1432 <- .e504/sigma
    .e1436 <- (.e51 * .e309 + .e1328) * .e15/.e18
    .e1442 <- .e1002 + sigma * (.e3 * (2 * (.e55 * .e113 * .e15 * 
        .e8 * .e5/.e134) - 1/sigma)/.e134 + .e1111)
    .e1448 <- (.e1279 * .e3 - .e412)/.e134 + .e1007 - .e283 * 
        .e16 * .e8/.e3
    .e1450 <- .e1007 + sigma * (.e3 * .e1339 - 2 * (.e879 * .e8))/.e135 - 
        .e284 * .e16 * .e8/.e3
    .e1452 <- .e240 * .e62
    .e1454 <- (.e28 + 2 * (.e236 * .e113 * .e5/.e32) + .e141) * 
        .e15/.e32
    .e1462 <- (.e1079 + .e102 - 1) * .e15 + .e736
    .e1463 <- .e1226 + .e76 * .e230
    .e1465 <- .e1029 + .e554 + 2 * (rho * .e601/.e3)
    .e1466 <- .e234 * .e40
    .e1468 <- .e1088 * .e8 + sigma * .e303 * .e5/.e567
    .e1469 <- .e76 * .e212
    .e1470 <- .e76 * .e62
    .e1473 <- .e1096/sigma
    .e1476 <- (.e1271 + 2 * (.e1 * sigma * .e8/.e566)) * .e4
    .e1477 <- .e460 * .e11
    .e1480 <- ((.e115 + rho * (6 * .e12 + rho * (2 * .e131 + 
        .e223 + 4 * .e8)/.e3)) * .e8 + .e146 * .e39)/.e3 + .e847 * 
        .e16 + .e222 * .e4/sigma
    .e1482 <- .e261 * .e4/sigma
    .e1486 <- .e140 * .e477
    .e1487 <- .e140 * .e522
    .e1489 <- .e140 * .e15/.e29
    .e1490 <- .e192 * .e40
    .e1491 <- .e1241/sigma
    .e1492 <- .e192/sigma
    .e1496 <- .e212 * .e4/sigma + .e672 * .e8 + .e1166/.e567
    .e1498 <- .e212/sigma + .e147/.e567
    .e1499 <- ((rho * (.e1139 - .e925) - (((.e1292 - .e232)/.e3 + 
        2 * .e507) * .e8 + .e129))/.e3 - (.e51 * .e330 + .e873) * 
        .e15/.e18) * .e15
    .e1501 <- .e241 * .e207 + (2 * .e670 - .e240 * .e66)/.e3
    .e1503 <- .e245 + rho * (.e731 * .e8/.e3 + 1)
    .e1507 <- (1 - 2 * (.e67 * .e113 * .e15 * .e8 * .e5/.e134)) * 
        .e3/.e134 + .e1092 - .e1002
    .e1508 <- .e126 * .e5
    .e1515 <- .e95 * .e215 + (2 * .e130 - .e1218)/.e3
    .e1516 <- .e1258 + .e118
    .e1517 <- .e3 * .e269
    .e1520 <- .e466 * .e8
    .e1521 <- .e1089 - .e789
    .e1522 <- .e646 * .e207
    .e1524 <- .e1092 + .e20 * (1 - .e950) * .e3/.e135 - .e1213
    .e1527 <- (.e853 + rho * (.e1019 + 3 + .e1160 - 3 * .e10)/.e3)/.e5
    .e1529 <- .e1262 + .e1121 + 2 * (.e84 + rho * ((.e164 - .e479)/.e3 - 
        .e257) * .e4/sigma - .e351)
    .e1531 <- .e36 * .e39 + .e927
    .e1532 <- .e642 + rho * (.e678 + .e1338) * .e4
    .e1538 <- (.e1293 - .e112)/.e3
    .e1539 <- .e320 * .e8
    .e1540 <- .e191 + .e430
    .e1545 <- (2 + rho * (.e843 - .e6) * .e4/.e98 - .e20 * (.e11 * 
        (2 * .e249 + .e777) * .e62 + .e1130) * .e3/.e29) * .e15/.e29
    .e1546 <- .e1106 - .e220
    .e1548 <- .e1106 + .e8 * .e1165 - .e20 * .e703/.e29
    .e1550 <- (.e951 - .e488)/.e3 - .e621
    .e1552 <- .e1301 * .e8 + sigma * (1/.e29 + .e20/.e567) * 
        .e15
    .e1553 <- (3 * .e125 + rho * (.e1019 + 4 + .e1160 - 4 * .e10)/.e3)/.e5
    .e1554 <- (.e15 * (rho - .e1071 * .e8) - .e1155)/.e3
    .e1556 <- .e8 * .e1315
    .e1558 <- .e8 * (.e102 - 2) - .e30
    .e1560 <- .e1276 - 1
    .e1562 <- (.e199 - 2)/.e205 - .e191
    .e1564 <- (.e6 - .e843)/.e3 - .e35 * (.e1120 + .e1127) * 
        .e15/.e29
    .e1565 <- .e276^2
    .e1566 <- .e322^2
    .e1568 <- .e4/.e205
    .e1569 <- .e12 + rho * (.e427 + rho * (.e446 + .e862 + 4 * 
        .e14))/.e3
    .e1571 <- 1/.e97 - .e621
    .e1572 <- .e933 - .e585
    .e1573 <- 1/.e5
    .e1574 <- 2 * .e970
    .e1575 <- 2 * .e359
    .e1577 <- 2 * .e346
    .e1578 <- 2 * (.e755/.e98)
    .e1579 <- 2 * .e820
    .e1582 <- 2 * .e26
    .e1583 <- 2 * .e918
    .e1584 <- 2 * (.e11 * .e516 * .e5)
    .e1585 <- 2 * .e832
    .e1586 <- 2 * .e668
    .e1587 <- 2 * (.e39/.e17)
    .e1588 <- 2 * .e241
    .e1589 <- .e369 - .e344
    .e1590 <- 2 * .e247
    .e1591 <- 2 * .e646
    .e1592 <- 2 * .e855
    .e1593 <- 2 * .e858
    .e1594 <- 2 * (.e15 * (.e1319 - .e915)/.e20)
    .e1595 <- 2 * .e577
    .e1596 <- 2 * .e382
    .e1599 <- 2 * (.e955 * .e4/.e20)
    .e1600 <- .e1291 - (.e839 + 4 * .e646)
    .e1601 <- 2 * .e1159
    .e1602 <- 2 * (.e42 * .e1147/.e20)
    .e1603 <- 2 * (.e1 * .e66 * .e15)
    .e1605 <- 2 * (.e1334 * .e5)
    .e1606 <- 2 * (.e38 * .e447 * .e5)
    .e1609 <- 2 * (.e55 * .e409 * .e62/.e134)
    .e1610 <- 2 * (.e323 * .e5/.e29)
    .e1611 <- 2 * (.e883 * .e113 * .e15 * .e8 * .e5/.e136)
    .e1613 <- 2 * (.e1168 * .e312 * .e8/.e136) + .e6
    .e1616 <- 2 + 2 * eps - ((2 * (.e1168 * .e15 * .e5/.e136) + 
        mu2 + .e7) * .e8 + 2 * .e1)
    .e1617 <- 2/.e468
    .e1619 <- 2/.e97 - .e621
    .e1620 <- 2/.e91
    .e1622 <- 3 * .e16
    .e1623 <- 3/.e27
    .e1626 <- 4 + .e1 * (.e493 + rho * (.e1101 - (.e936 + .e1143)) * 
        .e4/sigma) * .e15 - 2 * .e1163
    .e1627 <- 5 * .e14
    .e1628 <- .e518 - .e495
    .e1630 <- .e15 * .e414 * .e4
    .e1631 <- .e674 - .e508
    .e1632 <- rho * (.e887 + .e66 * .e92/.e3)
    .e1633 <- rho * .e801
    .e1634 <- rho * ((.e1239 - rho * .e768/.e18)/.e3 + ((.e8 * 
        .e582 - .e380)/.e3 + (.e354 + 2 * (.e11 * .e582 * .e5)) * 
        .e15/.e18) * .e15/.e18)
    .e1635 <- .e957 - .e1288
    .e1636 <- .e782 * .e15
    .e1637 <- .e1326 - .e637
    .e1638 <- rho * ((.e8 * (rho * (.e797 + 2 * .e167) - 2 * 
        .e211) - .e420)/.e3 + .e15 * (rho * (4 * .e79 + .e1302) - 
        4)/.e18 + rho * ((.e806 + .e820 * .e8 - .e1080)/.e3 + 
        .e997) * .e4/sigma)
    .e1639 <- rho * .e283
    .e1642 <- rho * .e284 * .e8 * .e4
    .e1645 <- rho * .e1082 * .e4/.e20
    .e1646 <- rho * (.e207 * .e574 + rho * (2 * .e1190 - 4 * 
        .e241))
    .e1647 <- rho * ((.e859 - .e546)/.e3 - .e422)
    .e1648 <- .e614/.e20
    .e1649 <- .e146 - .e8 * .e119
    .e1650 <- .e788 - .e382
    .e1651 <- rho * (.e937 + 2 * (.e35 * .e984/.e29))
    .e1652 <- rho * (.e610 + 2 * (.e35 * .e988 * .e15/.e29))
    .e1653 <- .e6/.e965
    .e1654 <- rho * .e415
    .e1656 <- .e615 * (.e1101 - .e1143) * .e15
    .e1658 <- .e615 * (2 * .e1147 - 4 * .e178) * .e15
    .e1659 <- .e372 * .e11
    .e1660 <- sigma * (.e1110 - .e3 * (2 * (.e1335 * .e62 * .e8/.e134) + 
        .e6/.e20)/.e134)
    .e1662 <- .e94 * .e1295/.e135
    .e1663 <- .e94 * .e8
    .e1664 <- .e94 * .e15
    .e1665 <- sigma * .e269
    .e1666 <- sigma * (2 * (.e43 * .e1025 * .e15/.e29) - 2 * 
        (.e76/.e94))
    .e1667 <- sigma * (3 - .e949)
    .e1668 <- .e20 * .e240
    .e1669 <- .e682 - (.e391 + .e640/.e27)/.e3
    .e1670 <- .e386 * .e75
    .e1671 <- .e741 * .e75
    matrix(c(c(c(c(2 * (rho^4 * .e534 * .e15/.e35), 
        .e1357, .e615 * (.e1267/.e94 + 2 * (rho * 
            ((.e1398 - .e1498)/.e94 + (.e1003/.e94 + .e1111 + 
                .e1662) * .e15) * .e4 - .e813) - 4 * .e1004) * 
            .e15, -(2 * (.e1 * .e1376 * .e15/.e43))), c(.e1357, 
        2 * (.e1 * .e534 * .e15/.e20), -(2 * (.e1 * 
            ((.e8 * (.e1151 - (2 * .e181 + .e946)) + .e951)/.e94 + 
                rho * (.e1442 * .e15 + (.e8 * .e1333 - (.e1492 + 
                  .e200))/.e94) * .e4 - 2 * .e617) * .e15)), 
        -(2 * (.e1148 * .e15/.e20))), c(.e1656/.e35, 
        .e1174, .e1 * (2 * (rho * ((((.e476 - .e1586)/.e94 + 
            2 * .e674) * .e8 + .e1473 - .e20 * .e269 * .e15/.e567)/.e94 + 
            (.e1274/.e94 + .e98 * (sigma * (2 * (rho * .e740 * 
                .e11 * .e113 * .e15 * .e5 * .e4/.e136) - 1) * 
                .e8 - .e6)/.e136 - (.e1271 + .e20 * (.e1129 + 
                .e1606) * .e15/.e29)/.e29) * .e15) * .e4) - 4 * 
            .e1049) * .e15 - .e1281, .e1354), c(2 * 
        (.e1 * .e1367 * .e15/.e43), 2 * (.e1149 * .e15/.e20), 
        2 * (rho * ((.e8 * (rho * (.e891 + 2 * .e791) - 
            2 * .e221) - .e1286)/.e94 + rho * (.e1423/.e94 + 
            (.e1486 - .e1666) * .e15/.e29) * .e4 + .e197 * (rho * 
            (.e1302 + 4 * .e155) - 4)/.e29) * .e15), -(2 * 
            (.e15 * .e1153/.e20)))), c(c(2 * 
        (.e615 * .e627 * .e15/.e43), .e1355, .e1 * 
        (.e1263/.e94 + 2 * (rho * ((.e1498 - .e1398)/.e3 - (.e1213 + 
            sigma * (.e1111 + .e1662)) * .e15) * .e4 + .e585 - 
            .e764) + 2 * .e1333) * .e15, -(2 * (.e1146 * 
        .e15/.e20))), c(.e1355, 2 * (.e1636/sigma), 
        -(2 * (.e42 * (rho * (((.e1550/sigma - .e617) * 
            .e8 + .e1492 + .e200)/.e3 - sigma * .e1442 * .e15) * 
            .e4 + .e621 - .e1260))), 2 * (.e1368/sigma)), 
        c(.e1 * .e1529 * .e15/.e43, .e1171, 
            2 * (rho * (.e931 + rho * (((((.e668 - .e476)/.e3 + 
                .e788 - .e382)/sigma - .e674) * .e8 - .e1473)/.e3 + 
                .e15 * (.e43 * (2 * .e1663 + rho * .e1616 * .e4)/.e136 - 
                  .e1394)) * .e4) * .e15), .e1348), c(2 * 
            (.e1634 * .e15/.e20), 2 * (.e1369/sigma), 2 * 
            ((.e1209 + rho * ((.e1546/sigma + (.e1669/sigma - 
                .e655) * .e8 - .e879 * .e92 * .e15/.e29)/.e3 + 
                sigma * (.e1487 + .e1666) * .e15/.e29) * .e4 - 
                .e682) * .e15), -(2 * (.e1360/sigma)))), 
        c(c(.e1658/.e35, .e1174, .e1 * 
            (2 * (rho * ((((.e1520 - .e269)/.e94 + .e20 * (.e1605 - 
                .e1129) * .e15/.e29)/.e29 - .e385 * .e3 * .e4/.e1112) * 
                .e15 + ((.e1089 + (rho * (.e592 - .e1620) * .e4 - 
                .e212)/sigma)/.e94 - .e94 * (.e949 - (2 + .e199)) * 
                .e15/.e135) * .e8 + ((.e410 - .e614)/.e98 + .e739)/sigma) * 
                .e4) - (2 * .e1049 + 2 * .e1521)) * .e15 - .e1281, 
            .e1354), c(.e1174, .e1602, 
            -(2 * (.e42 * (rho * ((.e1394 - .e497 * .e4/.e136) * 
                .e15 + ((.e1109 + (rho * (.e617 - .e1142) * .e4 - 
                .e192)/sigma + .e382)/.e94 - .e3 * (.e1609 - 
                .e247) * .e15/.e134) * .e8 + (.e639 + (.e1617 - 
                .e612)/.e98)/sigma) * .e4 - .e1273))), .e1350), 
            c(.e1626/.e35, -(2 * (rho * (.e512 + 
                rho * ((.e8 * (.e279 - (.e512 + .e776 + .e164)) + 
                  .e515)/.e3 + .e580 - (.e1121 + .e242)) * .e4/sigma) * 
                .e15/.e43)), (2 * (.e1 * (((.e872 - (.e642 + 
                .e230/sigma))/.e94 + .e530) * .e8 + rho * ((.e1096/.e98 + 
                .e530)/sigma + (.e1631/.e98 - .e98 * (.e1140 - 
                .e861) * .e15/.e136) * .e8 - .e1654 * .e3 * .e15 * 
                .e4/.e1113) * .e4 - (2 * (.e71 * .e411/.e29) - 
                (2 * (.e94 * .e62 * (.e576 - (.e38 * (2 + .e1590) - 
                  .e45)/.e5)/.e29) + .e1 * .e1560 * .e308/(.e965 * 
                  .e5)) * .e5) * .e15/.e29) * .e15) - .e1282) * 
                .e4 - .e1294, -(2 * ((.e856 + rho * ((.e8 * 
                (.e954 - (.e856 + .e940)) + .e420)/.e3 + .e870 - 
                (.e1595 + .e1289)) * .e4/sigma) * .e15/.e43))), 
            c(2 * (.e1638 * .e15/.e43), .e1594, 
                2 * ((.e637 + rho * (((.e1556 + (rho * 
                  .e1423 * .e4 - (2 * .e391 + .e1286))/sigma)/.e94 - 
                  ((.e1470/.e94 - .e95 * .e247 * .e15/.e5) * 
                    .e5 + (.e179 * .e119 + 1 + 2 * (.e55 * .e95 * 
                    .e3 * .e62/.e29) + rho * ((.e1558 * .e4/sigma - 
                    .e28)/.e3 + .e43 * (.e1287 - 2 * (.e26 * 
                    .e3 * .e948))/.e29)) * .e15)/.e29) * .e4 + 
                  .e644)) * .e15), .e1343)), c(c(2 * 
            (.e1 * (((.e1431 - .e105)/.e3 + .e971 - .e1000)/sigma + 
                .e963) * .e15/.e20), .e1352, rho * 
            (.e1270/.e94 + 2 * (rho * ((.e1004 * .e207 + (.e1469 - 
                2/.e205)/.e94 + sigma * (.e15 * (.e204 * .e3 * 
                .e4/.e136 - .e1430) - sigma * (.e518 + .e495 * 
                .e1295)/.e135)) * .e4 + .e652) - .e825) - 2 * 
                .e1040) * .e15, -(2 * (.e973 * .e15/.e20))), 
            c(.e1352, 2 * (.e15 * (rho * ((.e183 - 
                .e1377) * .e8/.e94 + .e375 * .e4/.e353) - ((.e1034 + 
                (rho * .e400 * .e16 - .e707)/.e3) * .e8 + .e809 + 
                .e288 - .e1633)/sigma)), -(2 * (.e15 * 
                (rho * (((((.e1490 + .e1094 + .e1156 + .e173)/sigma + 
                  .e200 - .e30 * .e1333) * .e8 - (.e1623 - .e676)/sigma)/.e3 + 
                  sigma * (.e1489 + sigma * (.e15 * (.e280 * 
                    .e4/.e135 - .e1430) - (.e1628/sigma + 2 * 
                    (.e467 * .e75 * .e11 * .e113 * .e15 * .e8 * 
                      .e5/.e134))/.e134))) * .e4 + .e656) - (.e1094 + 
                  .e173)))), -(2 * (.e1374/sigma))), c(rho * 
                (.e1269 + 2 * (rho * (((.e971 + (.e1035 - .e515)/.e3 - 
                  .e1000)/sigma + .e963) * .e4 + .e174 + .e160) - 
                  .e454) - .e1579) * .e15/.e43, -(2 * (.e15 * 
                (rho * ((((.e812 - .e270)/.e3 + .e805 - .e643)/sigma + 
                  .e963) * .e4 + .e174 + .e160) - .e763)/.e20)), 
                2 * (.e15 * (rho * (((.e1463 - .e1586)/.e94 + 
                  .e95 * .e1631 + rho * ((.e43 * ((.e1517 - .e8 * 
                  .e414) * .e15 - .e495 * (.e1611 - 1))/.e136 - 
                  2/.e532) * .e4 + 2 * (.e882 * .e8/.e136)) - 
                  (1 + rho * ((((.e310 + mu1 - (.e286 + y))/sigma - 
                    .e14) * .e8 + rho - .e755/sigma)/.e3 + .e247 * 
                    .e290/.e5 + 2 * (.e43 * .e460 * .e11 * .e5/.e29)) - 
                    (.e394 + 2 * (.e43 * (.e266 + .e96 * (.e1610 - 
                      .e1573)) * .e3 * .e62/.e29))) * .e15/.e29) * 
                  .e4 + .e551) - .e762)), .e1341), c(-(2 * 
                (.e1416 * .e15/.e20)), 2 * (.e1359/sigma), 
                -(2 * ((((.e8 * (rho * (.e1193/sigma + 
                  .e1040 * .e16) - .e1466/sigma) + rho * (sigma * 
                  .e240 * .e15/.e29 - .e702/sigma) - .e1091/.e205)/.e3 + 
                  .e20 * (.e1462/sigma + rho * (.e1503/.e94 + 
                    .e43 * (2 * (.e502 * .e11 * .e5) - 2 * (.e26 * 
                      (.e924 + .e945)))/.e29) * .e15)/.e29) * 
                  .e4 - .e1193) * .e15)), -(2 * (.e1364/sigma))))), 
        c(c(c(2 * (.e615 * .e591 * .e15/.e43), 
            .e1173, .e1 * (2 * (rho * (.e1524 * 
                .e15 + .e1552/.e94) * .e4) - (.e1267/.e3 + 4 * 
                .e1572)) * .e15, -(2 * (rho * ((.e892 + 
                .e548 + .e720)/.e3 + .e1145) * .e15/.e20))), 
            c(.e1173, .e1290, -(2 * 
                (.e42 * (rho * ((.e1619 * .e8 + .e946)/.e94 + 
                  sigma * .e1507 * .e15) * .e4 + .e621 - (.e8 * 
                  .e1321 + .e410)/.e3))), -(2 * (.e1420/sigma))), 
            c(.e1 * .e1600 * .e15/.e43, .e1170, 
                rho * (2 * (rho * (((.e8 * (rho * (.e1620 - 
                  .e685) * .e4 - .e382) + .e43 * .e269 * .e15/.e567 - 
                  .e1540)/sigma - .e1049 * .e8)/.e3 + .e15 * 
                  (.e467 * (1 - .e1611) * .e3 * .e4/.e136 - (.e1274/.e3 + 
                    (.e43 * (.e1120 - .e1606) * .e15/.e29 - .e6/.e94)/.e29))) * 
                  .e4) - .e1601) * .e15, .e1347), c(2 * 
                (rho * ((.e1239 - .e806)/.e3 - .e1436) * .e15/.e20), 
                2 * (.e1389/sigma), 2 * (.e15 * 
                  (rho * ((.e1548/.e94 - sigma * (.e1486 + .e1652) * 
                    .e15/.e29) * .e4 + .e92/.e97 - .e20 * .e265 * 
                    .e15/.e29) + .e20 * .e95 * .e15/.e29 - .e8 * 
                    .e1313/.e3)), -(2 * (.e1378/sigma)))), 
            c(c(.e1173, .e1290, rho * 
                (2 * .e1571 - (.e1263/.e3 + 2 * (rho * (.e1552/.e3 + 
                  sigma * .e1524 * .e15) * .e4))) * .e15, -(2 * 
                (.e1499/sigma))), c(.e1290, -(2 * 
                (.e591 * .e15)), 2 * (rho * ((((.e1623 - 
                .e488)/.e3 - 2 * .e621) * .e8 + .e946)/.e3 + 
                .e20 * .e1507 * .e15) * .e15 * .e4), 2 * 
                ((((.e274 - .e28)/.e23 - (((.e914 + 2 * .e665)/.e3 + 
                  2 * .e746) * .e8 + .e925))/.e3 - (.e606 * .e8/.e3 + 
                  (.e244 + (2 * .e594 + 2 * (.e445 * .e11 * .e5)) * 
                    .e15/.e18)/.e18 + .e1114 * .e3/.e32) * .e15) * 
                  .e15)), c(rho * (.e1591 - (.e1262 + 
                .e1291)) * .e15/.e20, .e1599, -(2 * 
                (rho * ((.e1564/.e29 + .e385 * .e1616 * .e4/.e136) * 
                  .e15 + (.e8 * (rho * (3/.e91 - 2 * .e685) * 
                  .e4 - (.e845 + .e1596)) - .e191)/.e3) * .e15 * 
                  .e4)), -(2 * (.e1499 * .e4/.e20))), c(2 * 
                (.e15 * (.e15 * (.e1328 - .e51 * .e582)/.e18 - 
                  (.e211 + (2 * .e865 - (.e402 + .e1292)/.e3) * 
                    .e8 + rho * (.e1139 - .e1041))/.e3)/sigma), 
                2 * (((.e1139 - ((((.e114 + .e1627)/.e17 - 
                  .e895)/.e3 - 2 * .e744) * .e8 + .e1041))/.e3 - 
                  ((.e244 + .e767)/.e18 + .e1245/.e32 - .e597 * 
                    .e8/.e3) * .e15) * .e15), 2 * (.e15 * 
                  (.e147 * (.e1652 - .e1487)/.e29 - (.e221 + 
                    (2 * .e682 - (.e391 + 2 * .e1158)/.e3) * 
                      .e8 + rho * (2 * (rho/.e97) - .e20 * .e768/.e29))/.e3) * 
                  .e4), -(2 * ((((.e100 - .e1254)/.e18 - 
                  (((2 * .e921 - .e453)/.e3 - (2 * .e759 + .e775)) * 
                    .e8 + (mu2 + rho * (.e1622 + .e114 + .e1627) - 
                    .e69)/.e23))/.e3 - (.e1253/.e32 - (((4 * 
                  rho - .e28) * .e8/.e3 + .e12) * .e16/.e3 + 
                  (2 * ((.e100 * .e166 + .e920) * .e5 - .e938) + 
                    2 * (.e609 * .e11 * .e5))/.e18)/.e18) * .e15) * 
                  .e15)))), c(c(.e1 * (2 * 
                (.e867 * .e4/sigma) - 4 * .e474) * .e15/.e43, 
                .e1170, rho * (2 * (rho * (((rho * 
                  .e1301 * .e4/.e20 - .e1089)/.e3 - .e1664 * 
                  (.e196 + .e1667)/.e135) * .e8 + ((.e7 - .e1520)/.e3 - 
                  .e43 * (.e1120 + .e1605) * .e15/.e29) * .e15/.e29 + 
                  .e785 + .e94 * (.e1331/.e1112 - .e303 * .e5/.e135)) * 
                  .e4) - .e1601) * .e15, .e1347), c(.e1170, 
                .e1599, -(2 * (rho * (((rho * .e1619 * 
                  .e4/.e20 - .e1273)/.e3 - .e375 * (.e196 + sigma * 
                  (2 - .e1609))/.e134) * .e8 + .e1564 * .e15/.e29 + 
                  .e3 * (.e1331/.e136 - .e182/.e134) + .e784) * 
                  .e15 * .e4)), -(2 * (.e1420 * .e4/.e20))), 
                c(.e42 * (rho * .e1600 * .e4/sigma - 
                  .e1592)/.e43, -(2 * (.e42 * (.e781 - 
                  (.e8 * .e1320/.e3 + .e1591)) * .e4/.e43)), 
                  2 * (.e42 * (rho * ((.e1650/.e98 - 
                    .e98 * .e15 * (.e196 + sigma * (4 - .e1140))/.e136) * 
                    .e8 + .e645 + .e98 * (.e1331/.e1113 - .e312 * 
                    .e5/.e136)) * .e4 - ((((.e1590 - .e1 * (1 - 
                    .e1276) * .e308/.e98)/.e5 - 2 * (.e20 * ((.e11 * 
                    (2 * .e1300 + .e777) - .e45)/.e5 - .e576) * 
                    .e3 * .e62/.e29)) * .e5 - 2 * (.e90 * .e411/.e29)) * 
                    .e15/.e29 + .e8 * .e958/.e3)) * .e4), -(2 * 
                    (.e15 * (rho * .e1211 * .e4/sigma - (.e8 * 
                      .e1325/.e3 + 2 * .e836)) * .e4/.e43))), 
                c(2 * ((.e95 * .e15/.e18 + rho * 
                  (((.e1105 + .e1303 * .e8 - .e806)/.e3 - .e1436) * 
                    .e4/sigma + .e761 - .e160) - .e8 * .e1311/.e3) * 
                  .e15/.e20), 2 * (.e1389 * .e4/.e20), 
                  2 * (.e15 * (rho * (.e1548 * .e4/.e98 - 
                    (.e1649 * .e4/.e3 - sigma * (.e16 * .e179 * 
                      .e8/.e3 + .e43 * (2 * (.e912 * .e26 * .e3) + 
                      .e1287)/.e29)) * .e15/.e29) - (((.e95 * 
                    .e1300 * .e15/.e5 - .e146 * .e62/.e3) * .e5 - 
                    2 * (.e67 * .e95 * .e3 * .e62 * .e15/.e29))/.e29 + 
                    .e1556/.e3)) * .e4), -(2 * (.e1378 * 
                    .e4/.e20)))), c(c(2 * (rho * 
                (.e1522 - (.e1212/.e3 + .e1632)) * .e15/.e20), 
                .e1345, (2 * (rho * (.e1572 * .e207 - 
                  rho * (.e212 * .e92/.e3 + sigma * ((.e1393 + 
                    .e965 * .e4/.e136) * .e15 + .e20 * (.e70 - 
                    (.e700 + .e22 * (2 * (.e1670 * .e3 * .e15 * 
                      .e8 * .e5/.e135) + .e137)))/.e135))) * 
                  .e4) - (.e1270/.e3 + 2 * .e1165)) * .e15, -(2 * 
                  (.e1401/sigma))), c(.e1345, -(2 * 
                ((.e794 + ((.e1063 - .e1377)/.e3 - .e1454) * 
                  .e8 + .e359 + .e706) * .e15)), -(2 * 
                (((.e410 - ((.e1094 + .e946 + rho * ((.e1568 + 
                  .e383)/.e3 + .e16 * .e1571 - .e656)) * .e8 + 
                  .e676))/.e3 - .e20 * (.e1489 + rho * ((.e1393 + 
                  .e790/.e135) * .e15 + (.e70 - (.e700 + .e22 * 
                  (2 * (.e882 * .e3 * .e15 * .e8 * .e5/.e134) + 
                    .e137)))/.e134))) * .e15 * .e4)), -(2 * 
                (((.e75 * .e1114 + .e628 * .e83)/.e32 - ((((.e254 * 
                  .e83 + .e916)/.e3 - (.e1216 * .e16 + .e690 + 
                  .e635)) * .e8 + .e1516/.e23 + .e758 + .e1247)/.e3 + 
                  ((.e267 - (.e1240 + .e361))/.e3 - (2 * .e986 + 
                    2 * (.e443 * .e11 * .e15 * .e5))/.e18) * 
                    .e15/.e18)) * .e15))), c((2 * 
                (rho * (.e1522 - .e1632) * .e4/sigma) - (.e1269 + 
                2 * .e1303)) * .e15/.e20, -(2 * (.e1233 * 
                .e4/.e20)), 2 * ((.e95 * .e1650 + rho * 
                ((.e511 - .e35 * (.e1517 * .e15 + rho * (.e70 - 
                  (.e700 + .e22 * (2 * (.e1671 * .e3 * .e15 * 
                    .e8 * .e5/.e136) + .e137))))/.e136) * .e4 - 
                  (((rho * (.e92 * .e4 + .e777) - .e1539)/.e3 - 
                    (.e290 * .e1300/.e5 + .e35 * (2 * ((.e26 * 
                      (.e1573 - .e1610) - .e405) * .e3 * .e62) + 
                      2 * (.e1477 * .e5))/.e29)) * .e15/.e29 + 
                    .e92 * .e230/.e3)) - .e1226/.e3) * .e15 * 
                .e4), -(2 * (.e1401 * .e4/.e20))), c(-(2 * 
                (((((.e227 * .e83 - .e747)/.e3 + .e16 * .e1303 - 
                  .e625) * .e8 + .e701 + .e1569/.e17 - .e1220)/.e3 + 
                  (.e1554 + .e15 * (rho * (.e1382 - 2 * (.e1384/.e18)) - 
                    2 * (.e1201/.e18)))/.e18) * .e15/sigma)), 
                2 * (((((.e1198 - .e361)/.e3 - (2 * .e1384 + 
                  .e1584)/.e18)/.e18 - .e273 * .e3/.e32) * .e15 + 
                  .e1516/.e1508 - (((.e589 * .e83 + .e754 + .e669)/.e3 + 
                  ((.e761 - .e744)/.e3 - .e1454) * .e16) * .e8 + 
                  .e968)) * .e15), -(2 * (((((.e234 * 
                  .e83 - .e748)/.e3 + .e16 * .e1165 - .e680) * 
                  .e8 + .e702 + .e1569/.e27 - .e1668 * .e15/.e29)/.e3 + 
                  .e20 * (.e1554 + .e15 * (rho * (.e1382 - 2 * 
                    (.e35 * (.e1204 - .e947) * .e26/.e29)) - 
                    2 * (.e35 * .e502 * .e11 * .e5/.e29)))/.e29) * 
                  .e15 * .e4)), -(2 * ((((rho * (2 + 2 * 
                  .e267 - .e362) - .e16 * (2 * (.e1258 + .e37) + 
                  .e7))/.e23 - (((.e261 * .e83 + .e926)/.e3 + 
                  .e1256 * .e16 - .e688) * .e8 + .e756 + .e1265))/.e3 - 
                  ((((((.e541 - (2 * .e37 + .e115))/.e3 + .e163) * 
                    .e8 - .e361) * .e16 * .e15 + rho * .e215)/.e3 + 
                    .e26 * (.e267/.e3 - .e1284) - 2 * (.e1411 * 
                    .e15 * .e5/.e18))/.e18 + (.e75 * .e1116 - 
                    .e628 * .e77)/.e32)) * .e15))))), c(c(c(2 * 
            (.e615 * .e15 * (.e734 - 3 * .e178)/.e35), .e1356, 
            .e314 * (rho * (.e1267/.e98 + 2 * ((.e1468 - 
                .e1648)/.e94 + (.e3 * (sigma * (2 - .e949) * 
                .e8 - .e6)/.e135 + .e1110 + .e1642/.e98) * .e15 - 
                .e709 * .e3/.e566)) * .e4 - 4 * .e1088) - .e1281, 
            -(2 * (.e42 * (rho * .e1376 * .e4/sigma - 2 * 
                .e389)/.e43))), c(.e1356, .e1602, 
            -(2 * (rho * ((.e639 + rho * .e1321 * .e4/.e98) * 
                .e8 + rho * (((.e3/.e134 + .e1639 * .e4/.e94) * 
                .e8 + .e171 + .e1660) * .e15 + ((1/.e468 - .e612)/.e20 + 
                .e191 + .e1304 * .e8)/.e94 - .e824 * .e3/.e566) * 
                .e4 - .e382) * .e15)), .e1350), c((.e1656 * 
            .e4/sigma - 3 * .e1163)/.e35, .e1172, (2 * 
            (.e1 * ((.e8 * (rho * (((.e476 - (3 * .e7 + mu2)/.e205)/.e94 + 
                .e674)/sigma - .e1338) * .e4 - .e642) + rho * 
                .e1096 * .e4/.e20 + .e1665 * .e303 * .e5/.e567)/.e94 + 
                .e15 * (rho * (.e1274/.e98 + .e94 * (sigma * 
                  (3 - .e1140) * .e8 - .e6)/.e136) * .e4 - (rho * 
                  (.e60/.e566 + .e1653) * .e4 + .e94 * (.e1130 - 
                  2 * (.e38 * .e62 * .e447))/.e29)/.e29) - .e844 * 
                .e3/.e566) * .e15) - .e1282) * .e4 - .e1294, 
            .e1351), c(2 * (.e42 * (rho * .e1367 * 
            .e4/sigma - 2 * .e661)/.e43), .e1594, 2 * 
            ((.e637 + rho * ((((.e1313/sigma + .e1326 - .e637) * 
                .e8 - (.e76 * .e1540 + .e1645))/.e94 + (2 + rho * 
                ((.e477 * .e4/sigma - .e446) * .e8/.e3 - sigma * 
                  (.e20 * (.e1583 + 2 * .e1659) * .e5/.e29 - 
                    .e1578)) - 2 * (.e71 * .e477 * .e62/.e29)) * 
                .e15/.e29 - .e503 * .e3/.e566) * .e4 + .e644)) * 
                .e15), .e1343)), c(c(2 * 
            (.e314 * .e1635/.e43), .e1171, rho * 
            (2 * .e1304 + rho * (.e1263/.e98 + 2 * ((.e1648 - 
                .e1468)/.e3 - (.e171 + .e1642/.e94 + sigma * 
                (.e3 * (.e1667 * .e8 - .e6)/.e135 + .e1110)) * 
                .e15)) * .e4) * .e15, .e1348), c(.e1171, 
            2 * (.e1636 * .e4/.e20), -(2 * (rho * 
                (((.e1086 + rho * .e1550 * .e4/.e20 - .e382) * 
                  .e8 + .e612/.e20 - .e191)/.e3 - .e15 * (.e1639 * 
                  .e8 * .e4/.e3 + sigma * (2 * .e283 + .e1660))) * 
                .e15 * .e4)), 2 * (.e1368 * .e4/.e20)), 
            c(.e42 * (rho * .e1529 * .e4/sigma - 
                .e1592)/.e43, -(2 * (.e42 * .e1635 * .e4/.e43)), 
                2 * (rho * ((.e1532 * .e8 + .e1329 * 
                  .e4/.e20)/.e3 + rho * (.e8 * .e1159/.e98 - 
                  .e20 * (((4 * .e94 + .e101) * .e8 - .e3 * .e1613) * 
                    .e15 + .e3 * .e312 * .e5)/.e136) * .e4 - 
                  .e1545) * .e15 * .e4), -(2 * (.e15 * 
                  (.e1306 - 2 * .e922) * .e4/.e43))), c(2 * 
                ((.e1206 + .e1634 * .e4/sigma - .e865) * .e15/.e20), 
                2 * (.e1369 * .e4/.e20), 2 * ((((.e637 + 
                  rho * (.e644 + .e1669 * .e4/.e20)) * .e8 + 
                  rho * (.e1546 * .e4/.e20 + .e92 * .e62 * .e5/.e29))/.e3 + 
                  (.e522 * .e249 + .e204 * (2 * (.e20 * .e918 * 
                    .e5/.e29) - .e1578)) * .e15/.e29) * .e15 * 
                  .e4), -(2 * (.e1360 * .e4/.e20)))), c(c((.e1658 * 
            .e4/sigma - 3 * (2 * (.e615 * .e178 * .e15 * .e4/sigma) - 
            2 * (2 + .e1603)))/.e35, .e1172, (2 * 
            (.e1 * (.e8 * (rho * ((.e8 * (.e785 - 2 * .e739) + 
                2 * .e667 - 2 * .e789)/.e98 - 2 * (sigma * .e212 * 
                .e3/.e722)) * .e4 - .e3 * (sigma * (.e11 * (3 - 
                2 * (.e35 * .e409 * .e303^2/.e135)) + 3 * .e287 - 
                .e344) - .e45) * .e5/.e135) + .e15 * (rho * ((.e466/.e98 - 
                .e3/.e566) * .e8 - .e1653) * .e4 - .e94 * (.e1130 + 
                2 * (.e1334 * .e62))/.e29)/.e29 + rho * (.e965 * 
                (5 * .e38 - .e45) * .e5/.e1112 - .e960/.e20) * 
                .e4) * .e15) - .e1282) * .e4 - .e1294, .e1351), 
            c(.e1172, 2 * (.e42 * .e1309 * 
                .e4/.e43), -(2 * (rho * (.e1545 + .e8 * 
                (rho * ((.e8 * (.e784 - 2 * .e639) + .e1596 - 
                  2 * .e964)/.e98 - 2 * (sigma * .e192 * .e3/.e722)) * 
                  .e4 - .e3 * (.e369 - (2 * (.e1335 * .e363/.e134) + 
                  .e344)) * .e5/.e134) + rho * (.e530 - .e959/.e20) * 
                .e4) * .e15 * .e4)), -(2 * (.e15 * (.e1317 - 
                2 * .e896) * .e4/.e43))), c((.e1626 * 
                .e4/sigma - 3 * (.e1163 * .e4/sigma - 2 * ((.e1603 + 
                4) * .e4/sigma - 2 * (.e42/.e17))))/.e35, -(2 * 
                (.e42 * (rho * ((.e8 * (.e279 - (.e63 + .e512 + 
                  .e46)) + .e105)/.e3 + .e580 - (.e84 + .e1121)) * 
                  .e4/sigma - 3 * .e1320) * .e4/.e35)), (2 * 
                (rho * ((.e3 * (2 * (.e363 + .e24 * (2 * .e1589 - 
                  .e395)) + 2 * .e881) * .e62/.e29 + rho * (2 * 
                  (.e247/(.e20 * .e5) - .e5/.e1565) + rho * (rho * 
                  .e1560 * .e4/(.e622 * .e5) - 2 * (.e1663 * 
                  .e5/.e1566)) * .e4) * .e15 * .e4) * .e5/.e29 + 
                  rho * (((rho * ((.e645 - 2 * .e530) * .e8 - 
                    2 * .e678) * .e4 - 2 * .e642)/.e98 - sigma * 
                    (((2 - .e1140) * .e312 + sigma * (4 * .e287 - 
                      .e344)) * .e5/.e136 + 2 * (.e230/.e722)) * 
                    .e3) * .e8 + rho * (.e622 * (6 * .e38 - .e45) * 
                    .e5/.e1113 - .e721/.e20) * .e4) * .e4) * 
                  .e15) - 8 * (.e20 * (5 * .e1119 + 512 * (.e35 * 
                .e952/.e462)) * .e4/.e1104)) * .e4 - 12/.e35, 
                -(2 * (.e15 * (rho * ((.e149 + .e8 * (.e778 - 
                  (.e129 + .e856)))/.e3 + .e870 - (.e1595 + .e256)) * 
                  .e4/sigma - 3 * .e1325) * .e4/.e35))), c(2 * 
                (.e15 * (.e1638 * .e4/sigma - 2 * .e1311)/.e43), 
                2 * (.e15 * (.e1319 - 2 * .e915) * .e4/.e43), 
                2 * (.e15 * (rho * (((.e8 * (rho * ((.e391 + 
                  .e1286) * .e4/.e98 - 2 * .e644) - 2 * .e637) - 
                  (.e1470 * .e5/.e29 + .e1645))/.e98 - 2 * (sigma * 
                  .e503 * .e3/.e722)) * .e4 - (rho * (((.e121 + 
                  1) * .e15 * .e8/.e3 - .e11/.e5) * .e179 + .e1558 * 
                  .e15 * .e308/.e98) - .e98 * (2 * (.e26 * .e929) + 
                  2 * (.e11 * .e62 * .e494))/.e29)/.e29) - (.e95 * 
                  .e411 + rho * .e76 * .e62 * .e4/.e98) * .e5/.e29) * 
                  .e4), -(2 * ((.e1278 - 2 * .e1221) * 
                  .e15 * .e4/.e43)))), c(c(2 * 
            (.e42 * (rho * ((.e1431/.e3 + .e971 - .e1000)/.e20 - 
                .e831) * .e4 - 2 * .e1432)/.e20), .e1349, 
            .e15 * (rho * (.e1270/.e98 + 2 * (.e1088 * 
                .e207 + rho * (.e1469 * .e4/.e98 - (.e346 + sigma * 
                ((.e75 * (sigma * (1 - .e949) * .e8 - .e6) + 
                  (2 * .e651 + .e1630) * .e8)/.e135 + .e1429 + 
                  .e94 * .e312 * .e5 * .e4/.e136))))) * .e4 - 
                2 * .e1637), .e1341), c(.e1349, 
            2 * (((((.e707 + rho * (.e183 - (.e1185 + .e742 + 
                .e160)))/.e3 - .e1034) * .e8 + .e1633 - .e809)/.e20 - 
                .e831) * .e15 * .e4), -(2 * (.e15 * (.e382 + 
                rho * ((.e1094 + .e1156 + .e173) * .e8 * .e4/.e98 - 
                  sigma * (.e3 * .e303 * .e5 * .e4/.e135 + 2 * 
                    .e588 + sigma * (.e1429 + (.e1628 * .e4/.e20 - 
                    2 * (.e528 * .e11 * .e409 * .e62 * .e8/.e134))/.e134))) - 
                ((.e191 + rho * (.e16 * .e1304 - .e1490 * .e4/.e20)) * 
                  .e8 + rho * .e848 * .e4/.e20)/.e3) * .e4)), 
            -(2 * (.e1374 * .e4/.e20))), c(.e15 * 
            (rho * ((.e1269 - .e1579)/sigma + 2 * (rho * ((.e971 + 
                .e1035/.e3 - .e1000)/.e20 - .e831) * .e4 - .e1432)) * 
                .e4 - .e1593)/.e43, -(2 * (.e15 * (.e1308 - 
            2 * .e1235) * .e4/.e20)), 2 * (.e15 * (rho * 
            ((.e1463/.e98 - (rho * (2 * (.e4/.e500) + .e20 * 
                ((4 * .e651 + .e1630) * .e8 - .e75 * .e1613)/.e136) + 
                .e98 * .e269 * .e312 * .e5/.e136)) * .e4 - (rho * 
                ((.e1539 - .e755) * .e4/.e20 + .e115) * .e15/.e3 - 
                (.e1589 * .e290 + .e98 * (2 * .e1477 + .e1583) * 
                  .e62/.e29))/.e29) - .e1532 * .e95) * .e4), 
            -(2 * ((.e1178 - 2 * .e658) * .e15 * .e4/.e43))), 
            c(-(2 * ((.e1416 * .e4/sigma - .e975) * 
                .e15/.e20)), 2 * (.e1359 * .e4/.e20), -(2 * 
                ((((.e16 * .e1637 + (rho * .e1193 - .e1466) * 
                  .e4/.e20) * .e8 - ((.e1091/.e27 + rho * .e234 * 
                  .e16) * .e4/.e20 + .e1452 * .e5/.e29))/.e3 + 
                  sigma * (.e502 * .e179 + sigma * (.e1462 * 
                    .e4/.e20 + rho * (.e1503 * .e15 * .e4/.e98 - 
                    2 * (.e587 * (.e11 * (2 * .e398 + .e863) - 
                      .e817) * .e5/.e29))))/.e29) * .e15 * .e4)), 
                -(2 * (.e1364 * .e4/.e20))))), c(c(c(2 * 
            (.e1 * .e1397 * .e15/.e43), .e1353, rho * 
            (.e207 * .e1135 + rho * (2 * ((.e813 + rho * (.e1450 * 
                .e15 + .e1496/.e94 + 2 * (.e879 * .e709/.e566))) * 
                .e4) - 4 * .e672)) * .e15, -(2 * ((.e389 + 
            rho * (((.e663 - .e389 * .e16) * .e8 + .e631 + .e711 + 
                .e1138 - .e554)/.e3 + .e1013)) * .e15/.e20))), 
            c(.e1353, 2 * (.e1395/sigma), -(2 * 
                ((.e95 * .e1321 + rho * (((.e488 - .e951)/.e94 + 
                  rho * ((.e1491 + .e846 * .e8 + .e791)/.e94 + 
                    sigma * (.e1448 * .e15 + 2 * (rho * .e824/.e566))) + 
                  .e617) * .e4 + .e656)) * .e15)), -(2 * 
                (.e1386/sigma))), c(.e1646 * .e15/.e43, 
                .e1169, (2 * (rho * (.e1049 * .e207 + 
                  rho * ((.e664 * .e8 + .e1562 * .e4/sigma + 
                    .e1166 * .e269/.e567)/.e94 + 2 * (.e879 * 
                    .e844/.e566) + .e15 * (.e20 * (.e280 * .e1340 + 
                    sigma * .e1280 * .e8) * .e4/.e136 - (.e28 * 
                    .e527/.e3 + (.e43 * (.e1123 - 2 * (.e1659 * 
                    .e447 * .e5))/.e29 - .e1476)/.e29)))) * .e4) - 
                  .e1585) * .e15, .e1342), c(2 * 
                ((.e661 + rho * ((.e1465 - ((.e661 * .e16 + .e625 + 
                  .e897) * .e8 + .e631))/.e3 - .e1414)) * .e15/.e20), 
                2 * (.e1358/sigma), 2 * ((.e891 + 
                  (.e748 - .e28 * .e1313)/.e3 + rho * ((.e1419/.e94 + 
                  sigma * (2 * (rho * .e503/.e566) - (.e16 * 
                    (4 - .e304) * .e8/.e3 + .e1651 + .e116 * 
                    (2 + 2 * .e477) * .e5/.e29) * .e15/.e29)) * 
                  .e4 + .e20 * (.e138 - (.e175 + .e492))/.e29) + 
                  .e587 * .e207/.e29) * .e15), -(2 * (.e1370/sigma)))), 
            c(c(2 * (rho * (.e474 * .e207 + 
                .e1647) * .e15/.e20), .e1346, (.e939 + 
                rho * ((2 * .e923 + 2 * (.e764 + rho * ((2 * 
                  (rho * .e724/.e3) - .e1496)/.e3 - sigma * .e1450 * 
                  .e15) - .e585)) * .e4 + 2 * .e846 - .e16 * 
                  .e1131 * .e8/.e3)) * .e15, -(2 * (.e1404/sigma))), 
                c(.e1346, 2 * (.e1103 * .e15), 
                  -(2 * ((.e1260 + rho * ((2 * (rho * 
                    .e1117/.e3) - (((.e1260 - .e621) * .e16 + 
                    .e771 - .e656) * .e8 + .e1491 + .e791))/.e3 - 
                    .e20 * .e1448 * .e15) - .e621) * .e15 * .e4)), 
                  2 * (((((.e1538 + 1 - .e121)/.e3 - (2 * 
                    .e590 + 2 * (.e26 * .e445 * .e11 * .e5))/.e18)/.e18 + 
                    (.e1115 * .e3 - .e1293)/.e32 - .e606 * .e16 * 
                    .e8/.e3) * .e15 + ((.e1075 - (.e1215 * .e16 + 
                    .e690 + .e635)) * .e8 + .e1057 + (.e1587 + 
                    2 * (rho * .e1055))/.e3 - .e903)/.e3) * .e15)), 
                c((.e512 + rho * ((.e1288 + 2 * 
                  (.e351 + .e1647 - .e84)) * .e4/sigma + .e1588 - 
                  .e16 * .e855 * .e8/.e3)) * .e15/.e20, -(2 * 
                  (.e1234 * .e4/.e20)), 2 * ((.e95 * 
                  .e1159 + rho * ((.e511 - .e43 * (rho * ((.e3 * 
                  .e1340 - sigma * .e489 * .e8) * .e15 + .e651) + 
                  .e1664 * .e8)/.e136) * .e4 - (((.e603 + 2 * 
                  (.e1329/.e3) - .e551) * .e8 + .e230 * .e4/sigma)/.e3 + 
                  .e1427))) * .e15 * .e4), -(2 * (.e1404 * 
                  .e4/.e20))), c(2 * (((.e58 * (2 - 
                  2 * .e582) * .e5/.e18 + .e16 * (3 - .e218) * 
                  .e8/.e3 + .e1327) * .e15/.e18 - ((.e16 * (.e865 - 
                  .e1206) - .e987) * .e8 + .e1029 + .e1242 + 
                  rho * (2 * .e1206 - .e903))/.e3) * .e15/sigma), 
                  2 * ((((.e1587 + 2 * (rho * .e1068))/.e3 - 
                    (((.e1217 - (.e174 + .e744)) * .e16 + .e754 + 
                      .e669) * .e8 + .e1203 + .e903))/.e3 - ((.e697 + 
                    .e765)/.e18 + (.e339 * .e16 * .e3 + (mu2 + 
                    .e523) * .e8)/.e32 - .e597 * .e16 * .e8/.e3) * 
                    .e15) * .e15), 2 * (.e15 * (.e20 * 
                    (.e16 * (3 - .e304) * .e8/.e3 + .e1651 + 
                      .e116 * (2 - 2 * .e522) * .e5/.e29) * .e15/.e29 - 
                    ((.e1227 + .e847/.e468)/sigma + (.e16 * (.e682 - 
                      .e1209) - .e1010) * .e8 + rho * (2 * .e1209 - 
                      .e587 * .e92/.e29))/.e3) * .e4), -(2 * 
                    (((((2 * .e358 - (.e12 + .e1622 + .e195) * 
                      .e39)/.e17 + 2 * (rho * .e1083))/.e3 - 
                      (((.e1228 - .e1061) * .e16 + .e1099 - .e688) * 
                        .e8 + .e26 * .e569/.e18 + .e1482))/.e3 - 
                      (((.e16 * (1 - (.e121 + (.e223 - .e69) * 
                        .e8/.e3)) * .e15 + rho * ((2 * .e448 + 
                        .e130)/.e3 - .e298))/.e3 - (2 * (((.e894 + 
                        .e909 * .e166) * .e5 - .e904) * .e15) + 
                        2 * (.e26 * .e609 * .e11 * .e5))/.e18)/.e18 + 
                        .e1230)) * .e15)))), c(c(rho * 
                (2 * (rho * .e1397 * .e4/sigma) - 4 * .e504) * 
                .e15/.e43, .e1169, (2 * (rho * 
                (.e1521 * .e207 + rho * ((((.e206 + .e961 + .e200)/.e98 + 
                  2 * (.e1332 * .e212/.e722)) * .e4 - sigma * 
                  ((.e3 * (3 * .e296 - .e286) * .e15/.e5 - 2 * 
                    (rho * .e303)) * .e5 - .e3 * (2 * (.e1670 * 
                    .e11 * .e113/.e135) + .e137) * .e303)/.e135) * 
                  .e8 + (.e1476 - (.e265 * .e466 + 2 * (.e43 * 
                  .e899/.e29))) * .e15/.e29 + (((.e614 - .e410)/.e98 - 
                  .e739)/sigma + .e1654 * .e75/.e1112) * .e4)) * 
                .e4) - .e1585) * .e15, .e1342), c(.e1169, 
                2 * (.e1395 * .e4/.e20), -(2 * 
                  ((.e8 * (rho * (.e1241/.e20 - (.e1273 - .e964) * 
                    .e16)/.e3 - .e639) + .e382 + rho * ((((.e181 + 
                    .e1156 + .e173)/.e98 + 2 * (.e1332 * .e192/.e722)) * 
                    .e4 - ((.e3 * .e320 * .e15/.e5 - 2 * (rho * 
                    .e62)) * .e5 - .e473 * (2 * (.e1167 * .e113/.e134) + 
                    .e137))/.e134) * .e8 + .e1427 + (((.e612 - 
                    .e1617)/.e98 - .e639)/sigma + .e583/.e136) * 
                    .e4)) * .e15 * .e4)), -(2 * (.e1386 * 
                  .e4/.e20))), c(.e15 * (.e1646 * 
                .e4/sigma - .e1593)/.e43, -(2 * ((.e95 * 
                .e1320 + rho * (.e1190 - .e1588)) * .e15 * .e4/.e43)), 
                2 * ((.e95 * .e958 + rho * ((.e760 + 
                  rho * ((.e664/.e98 + .e20 * (2 * (.e1329/.e722) - 
                    ((.e3 * (4 * .e296 - .e286) * .e15/.e5 - 
                      2 * (rho * .e312)) * .e5 - .e3 * (2 * (.e1671 * 
                      .e11 * .e113/.e136) + .e137) * .e312)/.e136)) * 
                    .e8 + ((.e1562/.e98 - .e530)/sigma + rho * 
                    .e386 * .e75/.e1113) * .e4)) * .e4 + (.e411 * 
                  .e290 + (2 * (.e20 * (.e814 + .e22 * (2 * .e840 - 
                  .e860)) * .e62/.e29) + rho * (((.e207 * .e8 + 
                  .e7)/.e322 + .e1332 * (.e3/.e5 + 2 * .e5) * 
                  .e8/.e1566) * .e308 + 2 * (((.e879/.e1565 - 
                  .e28/.e94) * .e4 - .e41)/.e5)) * .e15) * .e5)/.e29)) * 
                  .e15 * .e4), -(2 * (((.e389 + rho * .e1180) * 
                  .e4/sigma - (.e28 * .e1325/.e3 + 2 * .e663)) * 
                  .e15 * .e4/.e43))), c(2 * ((.e797 + 
                (.e747 - .e28 * .e1311)/.e3 + .e26 * .e207/.e18 + 
                rho * (((.e1465 - (.e987 * .e8 + .e631))/.e3 - 
                  .e1414) * .e4/sigma + .e669)) * .e15/.e20), 
                2 * (.e1358 * .e4/.e20), 2 * (((.e990/.e20 - 
                  .e28 * .e1315)/.e3 + rho * ((.e1419/.e98 + 
                  2 * (.e204 * .e503/.e722)) * .e4 - ((.e233 + 
                  .e16 * .e1649) * .e15 * .e4/.e3 + sigma * (.e179 * 
                  .e215 - .e1166 * (.e1123 + .e1287)/.e29))/.e29) - 
                  ((.e95 * .e320 * .e15/.e5 - .e1452/.e3) * .e5 + 
                    .e398 + .e478 - .e607 * .e290)/.e29) * .e15 * 
                  .e4), -(2 * (.e1370 * .e4/.e20)))), c(c(2 * 
                ((.e454 + rho * (.e1501 - ((.e504 * .e16/.e3 + 
                  .e1575) * .e8 + .e852 + .e953))) * .e15/.e20), 
                .e1344, (2 * ((.e825 + rho * (.e672 * 
                  .e207 + (2 * (.e4/(.e67 * .e5)) - .e240 * .e212)/.e3 - 
                  sigma * (.e1577 + 2 * .e619 + rho * (.e1373 + 
                    sigma * ((.e75 * .e1339 + sigma * .e659 * 
                      .e8)/.e135 + .e528 * .e4/.e136))))) * .e4) + 
                  2 * .e1010 - .e16 * .e1137 * .e8/.e3) * .e15, 
                -(2 * (.e1392/sigma))), c(.e1344, 
                -(2 * (((.e801/sigma + .e651/.e353) * .e4 + 
                  (.e1383 + (.e66 * (mu2 + rho * (.e446 + .e34 + 
                    .e31))/.e3 + .e16 * (2 * .e183 - ((.e706 + 
                    .e1575) * .e8 + .e852)))/.e3) * .e8 + .e793) * 
                  .e15)), -(2 * (((.e92 * .e848 + .e1568 - 
                  (((.e951 - .e676) * .e8/.e3 + 2 * .e1156 + 
                    .e946) * .e16 + .e383 * .e39/.e3 + .e791) * 
                    .e8)/.e3 - .e20 * (.e537 + .e1577 + rho * 
                  ((.e795 + .e75 * .e1279)/.e134 + .e1373 + .e880/.e135))) * 
                  .e15 * .e4)), -(2 * ((((2 * (rho * .e1076) - 
                  .e1531/.e17)/.e3 - ((.e1225 * .e8 + .e1057) * 
                  .e16 + (.e1248 * .e16 + .e254 * .e39/.e3) * 
                  .e8 + .e26 * .e510/.e18))/.e3 - ((.e659 * .e83 + 
                  .e75 * .e1115)/.e32 + ((.e941 - ((.e944/.e3 - 
                  .e163) * .e8 + .e121) * .e16)/.e3 - (.e1574 + 
                  2 * (.e26 * .e443 * .e11 * .e5))/.e18) * .e15/.e18)) * 
                  .e15))), c((2 * .e987 + 2 * ((.e454 + 
                rho * (.e1501 - (2 * .e174 + .e852 + .e953))) * 
                .e4/sigma) - .e16 * .e858 * .e8/.e3) * .e15/.e20, 
                -(2 * (.e1191 * .e4/.e20)), 2 * 
                  ((.e664 * .e95 + .e603 + rho * ((2 * (.e4/.e532) - 
                    .e43 * (.e75 * (rho * .e1340 + .e1665) + 
                      sigma * (.e70 + rho * (.e624 - (.e3 * (2 + 
                        .e357 - .e291) + .e581)/.e5)) * .e8)/.e136) * 
                    .e4 - (((2 * (sigma * .e39) - (.e16 * .e320 * 
                    .e8 + (rho * (.e1538 + 2 - .e121) - .e28) * 
                    .e4))/.e3 - .e320 * .e290/.e5) * .e15 - ((.e40/.e5 + 
                    2 * (.e35 * (.e435 * .e3 - .e58 * (2 * (.e281 * 
                      .e5/.e29) + .e141))/.e29)) * .e62 + 2 * 
                    (.e111 * .e460 * .e11 * .e5/.e29)))/.e29) - 
                    ((.e240 * .e230 + .e832 * .e16 * .e8)/.e3 + 
                      .e551)) * .e15 * .e4), -(2 * (.e1392 * 
                  .e4/.e20))), c(-(2 * ((((.e227 * 
                .e39/.e3 - .e987 * .e16) * .e8 + (.e1029 - .e975 * 
                .e8) * .e16 + .e1480/.e17 + 2 * (rho * .e898/.e3) - 
                .e1218/.e18)/.e3 + (.e1515 + rho * (.e1098 - 
                (.e1553 + 2 * (.e1412 * .e26/.e18))) - (.e58 * 
                (2 * .e486 + .e1582) * .e5/.e18 + .e298))/.e18) * 
                .e15/sigma)), 2 * ((.e1531/.e1508 + (.e1098 - 
                (.e26 * (2 * .e1412 + .e1584)/.e18 + .e1527))/.e18 - 
                (((.e1383 + ((.e273/.e17 + 2 * .e432)/.e3 + .e669 - 
                  .e1196)/.e3) * .e16 + .e589 * .e39/.e3) * .e8 + 
                  (.e589 * .e4/sigma + .e1223 * .e8/.e3) * .e16 + 
                  .e75 * .e273/.e32)) * .e15), -(2 * 
                ((((.e1480/.e468 + .e702 * .e4)/sigma + (.e234 * 
                  .e39/.e3 - ((.e913 + 2 * .e748)/.e3 + 2 * .e680) * 
                  .e16) * .e8 + 2 * (rho * .e1011/.e3) - .e1668 * 
                  .e26/.e29)/.e3 + .e20 * (.e1515 + rho * (.e1098 - 
                  (.e1553 + 2 * (.e35 * (.e1199 - .e910) * .e26/.e29))) - 
                  (.e298 + .e116 * (2 * .e502 + .e1582) * .e5/.e29))/.e29) * 
                  .e15 * .e4)), -(2 * ((((((2 - .e362) * 
                .e39 + .e248 + rho * (.e34 + 2 * .e368 + 2 * 
                .e306) - (.e1272 + 2 * (.e92 * .e39 + .e927)) * 
                .e16)/.e17 + 2 * (rho * .e1100))/.e3 - ((.e1237 * 
                .e8 + .e1482) * .e16 + (.e1266 * .e16 + .e261 * 
                .e39/.e3) * .e8 + .e26 * .e514/.e18))/.e3 - ((.e659 * 
                .e77 + .e75 * .e1132)/.e32 + (.e26 * ((.e34 + 
                rho * (.e248 + 4 * .e39)/.e3)/.e3 - (.e1574 + 
                2 * (.e1411 * .e5))/.e18) + (.e1098 - .e1527) * 
                .e16 + (2 * .e475 - .e163) * .e215)/.e18)) * 
                .e15)))))), length(y)) -> t2
t2[censoring,] = tt[censoring,]
t2
}

Try the cenGAM package in your browser

Any scripts or data that you put into this service are public.

cenGAM documentation built on May 1, 2019, 6:29 p.m.