Hi

I have an equation which i want solve it with DTM, but there are two problems:

1- How do I write a Dirak Delta function in maple for this equation?

2- I used pade approximation due to initial condition. Then how many sentenceses should be used to achieve the correct answer?

in this equation there are some constant parameter.

M:=0;A:=1;Phi:=0;r:=0.5;lambda:=0;phi[1] := 1.;  beta := .67;   phi[1] := 1.

then:

F[0]:=0;F[1]:=lambda;F[2]:=a/2;

That's my pleasure If there is somebody help me to code and solve this problem.

Thanks for every one very very much.

Hi

How can I write and solve a dirac_delta function in maple, for example in below equations:

F[0]:=0;
F[1]:=a;
F[2]:=b;

eq1:=sum((k-m+2)*(k-m+1)*F[k-m+2]*Dirac[m-1], m = 0 .. k);

eq2:=(k+1)*F[k+1]-Dirac[k];

for k from 0 to 10 do
F[k+3]:=solve(F[k+3]=(eq1+eq2)/((k+3)*(k+2)*(k+1)),F[k+3]);
end do;

hi

i want to solve a function which contains below series, but I can't.

SS:=sum(F[k-m]*sum(F[m-L]*sum(F[L-j]*F[j],j=0..L),L=0..m),m=0..k);

or

SS:sum(F[k-m],m=0..k)*sum(F[m-L],L=0..m)*sum(F[L-j]*F[j],j=0..L);

eq:=(-1/(k+1))*(F[k]+0.5*sum((k-m+1)*F[k-m+1]*F[m],m=0..k)+0.05*SS);

n:=8;
for k from 0 to n do
F[k+1]:=solve(eq);
end do;

with the first SS I have gotten a wrong nswer and with the second SS this error has been seen:

Error, (in solve) cannot solve expressions with sum(F[L-j]*F[j], j = 0 .. L) for F[j]

is there qny one hepl me please.

thanks