Note

Go to the end to download the full example code.

Create a quadratic functionΒΆ

In this example we are going to create a quadratic function of the form

f : \vect{x} \mapsto \mat{A} ( \vect{x} - \vect{b} ) + \vect{c} + \frac{1}{2} \vect{x}^T\tens{M}\vect{x}

import openturns as ot
import openturns.viewer as otv

create a quadratic function

inputDimension = 3
outputDimension = 2
center = [1.0] * inputDimension
constant = [-1.0, 2.0]  # c
linear = ot.Matrix(inputDimension, outputDimension)  # A
quadratic = ot.SymmetricTensor(inputDimension, outputDimension)  # M
quadratic[0, 0, 1] = 3.0
function = ot.QuadraticFunction(center, constant, linear, quadratic)
x = [7.0, 8.0, 9.0]
print(function(x))

[-1,56]

draw y1 with x1=2.0, x2=1.0, x0 in [0, 2]

graph = (
    ot.ParametricFunction(function, [1, 2], [2.0, 1.0]).getMarginal(1).draw(0.0, 2.0)
)
view = otv.View(graph)
y1 as a function of x0

Display all graphs

otv.View.ShowAll()