%autosave 10

from pulp import *

x = LpVariable('x', -10, 10)
y = LpVariable('y', -10, 10)
prob = LpProblem("Toy problem", LpMinimize)
prob += 3*x - y
prob.solve()

prob

(x.value(), y.value())

import sympy as S
from math import *
S.var("x mu sigma", real=True)

f = 1/(sqrt(2*S.pi)*sigma)*S.exp(-(x-mu)**2/(2*sigma**2))
f

ll = log(f)