Question: Fractional partial differential equations

I working on wavelets method to solve fractional partial differential equations. The problem when I solved and create a system of algebraic equation using Newton itertaion method the system not solved becuase its appear as amtrix forum when I substitute a collcetion points as follows :


for i from 1 to N do           # calculate left side of equ. at a collection points
    for j from 1 to N do
        ll[i,j]:=eval(q(T[i],T[j]));
    end do:
end do:
ll:
sys := []:
                   # create the system of equations
for i from 1 to N do
    for j from 1 to N do
    l[i,j]:= eval(ual(T[i],T[j])+ ubet(T[i],T[j])+uth(T[i],T[j])+ uwt(T[i],T[j]));
    end do:
end do:

for i from 1 to N do
    for j from 1 to N do
    sys:=[op(sys),(evalf(l[i,j]-ll[i,j]=0 ))]: 
    end do:
end do:
sys;

Coefficients:= fsolve(sys);

 

the results appears as:

[Vector[column](1, [0.92016313e-1+.38891735*c[1, 1]-.24496159000000003*c[2, 1]-.3889173499999999*c[1, 2]+.24496159000000003*c[2, 2]]) = 0., Vector[column](1, [-.9381081886+1.24805825*c[1, 1]-.66648505*c[2, 1]-.4702235500000001*c[1, 2]+.17656187000000007*c[2, 2]]) = 0., Vector[column](1, [2.131874035+.76110696*c[1, 1]-.12722801999999994*c[2, 1]-.76110696*c[1, 2]+.12722802*c[2, 2]]) = 0., Vector[column](1, [.355172452+1.5700572*c[1, 1]+.3444861000000001*c[2, 1]-0.47843280000000044e-1*c[1, 2]-.59894214*c[2, 2]]) = 0.]

 

please can you help me fix this problem 

Please Wait...