si=(((k 1((di/D))+((((ni/ei)))^(k 2))))^(-(Subscript[k, 3])))