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
