Apprendre Maple

[Bibinou]

Bonjour,
j'ai fais une procédure permettant de faire une régression linéaire.
Ce que je voudrais pouvoir faire c'est rajouter qq chose pour pouvoir estimer l'erreur de la mesure en prenant l'intervalle de confiance à 95% [/x - delta x, /x + delta x] puis qui m'affiche x = /x +- delta x

voilà ma procédure

regression_lineaire:=proc(L::list,h)
> local xk,yk,xkyk,xk2,a1,b1,n,a,b,C1,C2,graph,Lx,Ly,i,j,f,q,F,Q:
> global m:
> with(plots):
> n:=nops(L):
> m:=nops(L):
> Lx:=[]:Ly:=[]:
> Lx:=[op(Lx),seq(L[j,1],j=1..n)]:
> f:=min(op(Lx));
> q:=max(op(Lx));
> Ly:=[op(Ly),seq(L[j,2],j=1..n)]:
> F:=min(op(Ly));
> Q:=max(op(Ly));
> xk2:=sum(('L[i,1]')**2,'i'=1..n):
> yk:=sum('L[i,2]','i'=1..n):
> xk:=sum('L[i,1]','i'=1..n):
> xkyk:=sum(('L[i,1]')*('L[i,2]'),'i'=1..n):
> a1:=(xk2*yk-xk*xkyk)/(m*xk2-(xk)**2):
> b1:=(m*xkyk-xk*yk)/(m*xk2-(xk)**2):
> a:=evalf(a1);
> b:=evalf(b1);
> C1:=plot(a+b*x,x=f..q,y=F..Q,color=green):
> C2:=plot(L,x=f..q,y=F..Q,style=point):
> graph:=display(C1,C2);
> if nargs=2 then print(graph)
> fi:RETURN(a,b)
> end:

merci