with matlab from pk