# -*- coding: utf-8 -*-
"""
Created on Thu Nov  3 04:45:41 2016

@author: пользователь
"""
from sympy import Symbol, solve, lambdify, Matrix
k1 = Symbol("k1")
k1m = Symbol("k1m")
k3 = Symbol("k3")
k2 = Symbol("k2")
k2m = Symbol("k2m")
k2y=Symbol("k2y")
x = Symbol("x")
y = Symbol("y")
eq1 = k1*(1 - x - y) - k1m*x - k3*x*y
eq2 = k2*(1 - x - y)**2 - k2m*y**2 - k3*x*y
res = solve([eq1, eq2],k2,x)
#print(res[0][0])
eqprek2=res[0][0]
#print(res[0][1])
#print(jacA)
#print('PPPPPP')
A = Matrix([eq1, eq2])
var_vector = Matrix([x, y])
jacA = A.jacobian(var_vector)
#print(det_jacA)
eqxy=res[0][1]
det_jacA = jacA.det()
A = Matrix([eq1, eq2])
var_vector = Matrix([x, y])
jacA = A.jacobian(var_vector)

eqk2=solve(det_jacA.subs(x,eqxy),k2)
#print('TTTTTTT')
#print(eqprek2)
#print(type(eqprek2))
eqprek1=eqprek2-eqk2[0]
#print(eqprek1)
eqeq=y*(k1+k1m+k3*y)*(k1*k2m*y-k1*k3*y+k1*k3+k1m*k2m*y+k2m*k3*y**2)/((k1m+k3*y)**2*(y-1)**2)-(k1+k1m+k3*y)**2*(k1**2*k3*y/2-k1**2*k3/2+k1*k1m*k3*y/2-k1*k1m*k3/2-k1*k2m*y*(k1+k1m+k3*y)+k1*k3*y*(k1+k1m+k3*y)/2-k1m*k2m*y*(k1+k1m+k3*y)-k2m*k3*y**2*(k1+k1m+k3*y))/(k1**3*k3*(y-1)**2+k1**2*k1m*k3*(y-1)**2+k1**2*k3**2*y*(y-1)**2+k1*k1m*y*(k1+k1m+k3*y)**2-k1*k1m*(k1+k1m+k3*y)**2+k1*k3*y*(k1+k1m+k3*y)**2-k1*k3*(k1+k1m+k3*y)**2-k1m*y*(k1+k1m+k3*y)**3+k1m*(k1+k1m+k3*y)**3-k3*y**2*(k1+k1m+k3*y)**3+k3*y*(k1+k1m+k3*y)**3)
eqk1=solve(eqeq,k1)
#print(eqk1)