1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28
|
#! /usr/bin/env python
from __future__ import print_function
from openturns import *
from math import *
TESTPREAMBLE()
try:
# Analytical construction */
input = Description(1)
input[0] = "x"
output = Description(1)
output[0] = "y"
formulas = Description(output.getSize())
formulas[0] = "x-cos(x)"
analytical = NumericalMathFunction(input, output, formulas)
myAlgo = Secant()
print("myAlgo=", myAlgo)
value = 0.0
print("Solve ", formulas[0], "=%.5f" % value, " for ", input[0])
print("x=%.5f" % myAlgo.solve(analytical, value, 0.0, 3.0))
print("myAlgo=", myAlgo)
except:
import sys
print("t_Secant_std.py", sys.exc_info()[0], sys.exc_info()[1])
|