Dear Mr. Erik Postma, Mr. Georgios Kokoridis, Mr Patrick T and Mr. Dave Linder,

I would like to thank you for all your help and support.

i have Maple 9.5 version and it seams that pakages like Statistic do not function. so i tried to solve the problem starting fromm your given help and finaly i was able to pull it out. 

below is a solution that i obtained which has a graph fits Mr. Kokovidis solution.

Thank you all. God bless you.

Best regards

Tom

>> restart;
>>with(Optimization):
>>data:=[[0,5],[2,1],[3,0.6]];
>>eq := a*exp(k*cos((x)-mu))/(2*Pi*BesselJ(0,k));        
>>residues:=map( proc (d) options operator, arrow; eval(eq, x = d[1])-d[2] end proc,data );     
>>sol := LSSolve(residues);
 >> poly := eval(eq, sol[2]);
  >> with(plots):
> >p1 := pointplot(data):
>> p2 := plot(poly, x=0..3):
>> display(p1, p2);