@Markiyan Hirnyk ...  should we follow that ? Surely not !

But what makes you think that what's done is incorrect?

In fact the term "p value" and some related notations are a kind of paradigm in statistics

Regards

@Carl Love
Your solution is undoubtedly more concise, "automatic" (and elagant) than my proposal and goes quickly to the essential.
When I read it I was debating whether I should publish "mine" ?

But I thought it could be interesting for Tazatel to understand why the problem is in fact easy to solve. Let say that I was more interested by explanation and pedagogy than by efficiency ... or elegance

I feel like you should have publish your whole approach, from the "manual" treatment to the final elegant solution, for your phrase "The implicit nature of the model can be handled by subtracting one side from the other and using 0 as the dependent variable values." might be thought of as "pulling a rabbit out of a hat" (although you declare "The model is linear in the parameters a and b")

I am not a Maple pro and I have often a lot of difficulties to understand some solutions that, sometimes, look to me as an exercise in style.
On the one hand I agree with you thet "the essence of elegance in programming is to reduce such redundancy".
But on the other hand the best is the enemy of good ... and elegance can be the enemy of readibility, maintenability and shareability (think to people that develop codes used by others).

Your last remark makes me smile : until I came here I had never use the tilde but only map and zip.
It is people like you who have showed me this "trick" I am now a fan of (but when I develop code for others, I keep prefering "map" ... )

@Carl Love 

First : In french the term is "p valeur", a word to word translation of "p value". But a lot of statistician use in fact "p value" even in discussions in french

Second : the summary Maple 2016 provides is an exact copy (even if it is not explicly mentioned) of the output of the "lm" function of the statistical language R. 
R being a world reference in Statistics, I guess we can accept blindly the structure of the summary Maple returns.
Note : even the notation P(>|t|) is used in R ; so I do not  think it would be suitable to replace it

kind regards

@Thomas Richard 
I forgot to say I work with Maple 2015.
Interesting to know this limitation has been lifted !

Thank you for the answer.

@Carl Love Sorry Carl for having posted "my" solution.
I had already prepared an answer but I had i little family emergency to deal with ... when I got back I posted the answer without looking if anyone else had already done the job.

In this case my answer presents no great interest.

@cskoog

I'm ashamed to see that everyone else but me seems to know this printlevel command 

@Carl Love 

A quick reply to a useful answer ; thanks

@acer 

from the "do help page" : 

The setting of printlevel causes the display of the results of all statements executed up to the level indicated by printlevel. By default, printlevel is initially set to 
1
. As a result, if nested loops are being used, it may be necessary to set printlevel to a higher value in order to see the results of all the statements in the loops.

Thank you for this information


postcsript : 
I was desperatly looking for the good help page. 
Couldn't it be judicious to refer this "do help page" in the "; help page" ?

@Carl Love  ... I do understand now a little bit more the roles of savelib and march.

I will dig a little bit deeper for my part.

Thanks again

@Carl Love 

thanks Carl.

I am not sure I have clearly understood all you have said : a lot of things are quite new for me.


Just a last question before I start digesting all this material : 

When I construct a module and I want to archive it, I write in a worksheet a code like the one below (it is probably a pitiful method)

# 1) definition of the module

MyModule := module ()
export F1, F2, ....F10;
option package;
F1 := proc(...)
...

end proc;
F2 := proc(...)

...

end proc;
...
end module;

# 2) Archiving

MyArchive := "/.../MyModule.mla";
march(`create`, MyArchive, 10);
savelib(MyModule, MyArchive);


I have no .m file here.
But you write : "A library is stored in an .mla file. A .mla file is constructed from .m files by using march."
Did I missed something ?

@Carl Love ... but  I am a little bit confused by your answer

I use to write modules that I save in  .mla files which I manage with the march command.

Of course I should have read more carefully the march help page and more specifically all which concerns .m files.
But, before your answer I thought naively that the concept of "archive" was related ONLY to mla files and the concept of "library"

Here is my first source of confusion : probably because some background in numerics I associate the idea of library to a collection of numerical routines or procedures (as for the NAG library for instance).
I naively thought that, in Maple, the term "library" had the same meaning  whereas it seems to have a broader sense and that a library may contain other things than just procedures  ... ???
I have to confess here that I am probably confused between the terms archive and library ?

An other source of mistake is my knowledge of Matlab, where .m files generally contain structures data, not code.


If it is not much to ask, could you shed some light on these issues for me ?

Thanks in advance

@tomleslie 
Huge thanks to you Tom 

@Carl Love I have just browsed all the questions I have asked to see that I hadn't gave you a reply :-(

With all my apologies here it is :

Quite strangely the set of instructions returns no error at all !!!

I'm not sure what to think of this ? It is likely that I did some mistake when I answered you for the first time.

Even if I am probably not the only one to "forget" about answering I'm sorry about that. I'll try to do better the next time.

Have a good day

@vv 

OSX El Capitan

@doduythao 

VeHinhNemXien := proc(Alpha,Vbd)
local Y,V0,alpha,X,ball,Xmax,bgr;
uses plottools, plots;
Y := unapply(V0*sin(alpha)*X/(V0*cos(alpha)) - 1/2*9.8*(X/(V0*cos(alpha)))^2,alpha,V0,X);
ball := proc(x,y) plots[pointplot]([[x,y]],color=red,symbol=solidcircle,symbolsize=40) end proc;
Xmax := 2*Vbd^2*sin(Alpha)*cos(Alpha)/9.8;
bgr := plot(Y(Alpha,Vbd,X),X=0..Xmax,linestyle=[2]);
animate(ball,[X,Y(Alpha,Vbd,X)],X=0..Xmax,scaling=constrained,labels=["Độ xa","Độ cao"],frames=60,background=bgr);
exportplot(FileTools:-JoinPath([FileTools:-TemporaryDirectory(), "dothi.gif"]), animate(ball,[X,Y(Alpha,Vbd,X)],X=0..Xmax,scaling=constrained,labels=["Độ xa","Độ cao"],frames=60,background=bgr), gif);
end proc:

VeHinhNemXien(1., 2.)

Here the result (screen capture of the gif file I wasn't capable to load)