Advertisement
Guest User

Untitled

a guest
Jun 19th, 2019
75
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.45 KB | None | 0 0
  1. import scipy
  2. from scipy.optimize import minimize
  3.  
  4.  
  5. lista=[0.002,0.006,0.003,0.02,0.008,0.006,0.05]
  6.  
  7. def assign(k):
  8. return [0.01 if 0<=x< k*0.005 else 0.02 if k*0.005 <= x <k*0.01
  9. else 0.05 if k*0.01<=x<k*0.025 else 0.1 for x in lista]
  10.  
  11. def constraint(k):
  12. return sum(assign(k))-0.2
  13.  
  14. def fun(k):
  15. return k
  16.  
  17.  
  18. k0=0
  19. bnds=[(0,100)]
  20. cons={'type':'eq','fun':constraint}
  21. res=minimize(fun,k0,bounds=bnds,method='SLSQP',constraints=cons)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement