Extract results held in Molpro variables

from pymolpro import Project

p=Project("molpro_variables")
p.write_input("""
do i=1,73
x(i)=i*10, cosine(i)=cos(x(i)), sine(i)=sin(x(i))
enddo
put,xml
""")
p.run(wait=True)
assert p.status=='completed' and not p.errors()

values={}
for name in ['x','sine','cosine']:
    # values[name]=[float(v.text) for v in p.xpath('//variable[@name="'+name.upper()+'"]/value')]
    values[name]=p.variable(name)

import matplotlib.pyplot as plt
plt.plot(values['x'],values['sine'],values['x'],values['cosine'])
plt.legend(['Sine','Cosine']);