# Import numerical libraries
import numpy as N
from pylab import *
# Import the JModelica.org Python packages
from jmodelica.jmi import compile_jmu
from jmodelica.jmi import JMUModel
jmu_name = compile_jmu("opt","ocp.mop")
# Load the dynamic library and XML data
dyn = JMUModel(jmu_name)
dyn.initialize()
opts = dyn.optimize_options()
opts['n_e'] = 5 # number of elements
opts['n_cp'] = 3 # number of collocation points
opts['result_mode']='default'
opts['IPOPT_options']['max_iter']=100
opt = dyn.optimize(options=opts)
t=opt['time']
print t[-1]
plot(t,opt['sim.y'])
print opt['sim.y']
import os.path
# Import numerical libraries
import numpy as N
from pylab import *
# Import the JModelica.org Python packages
from jmodelica.jmi import compile_jmu
from jmodelica.jmi import JMUModel
jmu_name = compile_jmu("opt", "ocp.mop")
# Load the dynamic library and XML data
dyn = JMUModel(jmu_name)
dyn.initialize()
opts = dyn.optimize_options()
opts['n_e'] = 5 # number of elements
opts['n_cp'] = 3 # number of collocation points
opts['result_mode'] = 'default'
opts['IPOPT_options']['max_iter'] = 100
opt = dyn.optimize(options=opts)
t = opt['time']
print t[-1]
plot(t, opt['sim.y'])
print opt['sim.y']
print opt['sim.a']
