def sb03md_example():
from numpy import zeros
A = array([[3, 1, 1], [1, 3, 0], [0, 0, 3]])
C = array([[25, 24, 15], [24, 32, 8], [15, 8, 40]])
out = slycot.sb03md57(A, C=C, dico='D')
print('--- Example for sb03md ---')
print('The solution X is')
print(out[2])
print('scaling factor:', out[3])
def sb03md(n, C, A, U, dico, job='X',fact='N',trana='N',ldwork=None):
ret = sb03md57(A, U, C, dico, job, fact, trana, ldwork)
return ret[2:]