numerical method using matlab