def main():
print('Differentiation:')
x, y = map(Symbol, 'x y'.split())
expr = (x - acos(y))*exp(x + y)
Dexpr = expr.diff(y)
print(Dexpr)
print(Dexpr._print_ccode())
def _get_array():
X, Y, Z = inp = array.array('d', [1, 2, 3])
args = x, y, z = se.symbols('x y z')
exprs = [x+y+z, se.sin(x)*se.log(y)*se.exp(z)]
ref = [X+Y+Z, math.sin(X)*math.log(Y)*math.exp(Z)]
def check(arr):
assert all([abs(x1-x2) < 1e-13 for x1, x2 in zip(ref, arr)])
return args, exprs, inp, check
