"Numerical optimization
The function f(a,b) = a^2*e^(a*b) + 3*a*ln(b) is
calculated by:
x1 = a^2
x2 = a*b
x3 = e^x2
x4 = ln(b)
x5 = a*x4
x6 = 3*x5
x7 = x1*x3
f = x6 + x7
Use forward-mode automatic differentiation to find f(0.3,0.5) and
D_p f(0.3,0.5), where p = [a b]^T = [1 2]^T.