showgan@tg says to YSITD
可以用矩陣找r