@Carl Love you are right, but these functions is the elements of matrices and when I multiply it I got one element as a combination of functions as you can see. 

@Carl Love thanks alot for your answer. Unfortunatly, I can't do it becuase these numbers change in each trying. And it is result of multiplications of vectors and matrices , therefore I can't represnting as functions.

@tomleslie thanks alot for your response of my question. Is it ok in general but this not correct, because as you can see in the program the result of approximate not 0.24 exactly is   [.24000000001864324]. So there is an error as a result of the figure.

@acer unfortunetly not all examples make the same sense, I will try to fix it by searching another way because not all examples contains same number of conditions and same number of equations may be the degree of system will be 12 by 12 or bigger

@tomleslie thanks alooooot. I will check it 

thank you very much

@acer Hi, because the full code has alot of details related special work and not allowed to share. Also, the problem just as i mentioned. And this part gave the equations above 
for i from 1 to N-1 do
    l:= eval(1/sqrt(t+1) * yvv[1,i] +4*uu[1,i]+ 4*(t+1)*uu1[1,i], t = T[i]);
     ll:=eval(q, t = T[i]);
    sys:=[op(sys),(evalf(l-ll=0))];
end do;

ics:=eval(yvv[1,1]=sqrt(Pi),t=0),eval(uxx[1,1]=sqrt(Pi) /2 ,t=0),eval(yvv[1,1]=sqrt(2*Pi),t=1);

@Kitonum thanks alot for your response. I know if we use this way, but my problem is : 

1/ how to collect the equations together as you did by yourself because the result appeared separated

2/ the result is not like this if we deal with the conditions alone and get one equation only then put it with system it will be solve correctly. I did by my self because it is just three equations so it is easy. but if more it difficult so what we can do 

sys := [.8164965809*c[1]+101.3127271*c[3]-2.479560632 = 0., .7071067810*c[1]+1.414213562*c[2]+146.5538537*c[3]-2.772453851 = 0., 4/3*c[1]+c[2] = (1/24)*sqrt(Pi)*(19+16*sqrt(2))]; ics := c[1]-2*c[2]+3*c[3] = sqrt(Pi), 4*c[2]-16*c[3] = (1/2)*sqrt(Pi), c[1]+2*c[2]+3*c[3] = sqrt(2)*sqrt(Pi); 4/3*c[1]+c[2] = (1/24)*sqrt(Pi)*(19+16*sqrt(2)); fsolve(sys);
{c[1] = 2.187662474, c[2] = 0.1573948503, c[3] = 0.00684357949}

@Kitonum 

sys := [.8164965809*c[1]+101.3127271*c[3]-2.479560632 = 0., .7071067810*c[1]+1.414213562*c[2]+146.5538537*c[3]-2.772453851 = 0.]

ics := c[1]-2*c[2]+3*c[3] = sqrt(Pi), 4*c[2]-16*c[3] = (1/2)*sqrt(Pi), c[1]+2*c[2]+3*c[3] = sqrt(2)*sqrt(Pi)

and it can be bigger than this form nut this as a example