我有一个带上下界约束的运输问题模型，问题比较复杂，无法用表上作业法求解，想用单纯形法，不知怎么写算法和编程实现（matlab）。其中，ijt是下标。
min CijtXijt
    Xijt=dj
    a=<Xijt=<b
    Xijt=<Si
    Xijt>=0