val x by symbol
    val y by symbol

    val normal = DerivativeStructureExpression {
        val x = bind(x)
        val y = bind(y)
        exp(-x.pow(2)/2) + exp(-y.pow(2)/2)
    }

    val startingPoint: Map<Symbol, Double> = mapOf(x to 1.0, y to 1.0)

    @Test
    fun testOptimization() {
        val result = normal.optimize(startingPoint)
        println(result.point)
        println(result.value)
    }