X = [41 5 2.5 -83.5 1000];
M = model_v3(X);