Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- from numpy import sqrt
- from scipy.optimize import minimize
- from numpy import pi
- V_0 = 200.0
- z=12
- y=15
- def V(x):
- r,h = x
- v = pi*r*r*h
- return v
- def lim1(x):
- return V(x) - V_0
- def lim2(x):
- return min(x)
- def S(x):
- r,h = x
- return 2*pi*r*h + 2*pi*r*r
- def P(x):
- r,h,y=x
- P=(4*pi*r+h)*y
- def Price(x):
- z,r,h,y,S=x
- Pr=(S*z)+((4*pi*r+h)*y)
- return Pr
- cons = ({'type': 'eq', 'fun': lim1},
- {'type': 'ineq', 'fun': lim2})
- r0 = 5
- h0 = V_0/(pi*r0*r0)
- x0 = (r0,h0)
- sol = minimize(S,x0, constraints = cons)
- #print sol.x
- #print V(sol.x)
- print S(sol.x)
- #print Price
- print P
Advertisement
Add Comment
Please, Sign In to add comment