agim

0 Reputation

2 Badges

14 years, 221 days

MaplePrimes Activity


These are replies submitted by agim

Thanks it works very fast and it gives the results that I hope to see

 

Again thanks for your kindly help

Thanks it works very fast and it gives the results that I hope to see

 

Again thanks for your kindly help



 

 

restart:
#with(plots):

lambda:=-1; alpha:=0.5;  N:=1;  m:=ceil(alpha);

(1)

u[0]:=0;      #initial conditions

(2)

for i from 0 to N do
DI1:=simplify((1/GAMMA(m-alpha))*int((x-t)^(m-alpha-1)*diff(u[i],[t$m]),t=0..x));
DI1:=subs(x=t,%);
DI2:=DI1-(t^(1/2)/evalf(Pi)^(1/2))*(u[i]+1)^2;
u[i+1]:=u[i]+simplify((1/GAMMA(alpha))*int(lambda*(x-t)^(alpha-1)*DI2,t=0..x));
u[i+1]:=subs(x=t,%);
end do:

v:=subs(t=x,u[i]);

(3)

 



Download example.mws

 

also when I run this file one after onather it create different results. What is the problem, I do not understand

I coded it in Classic Worksheet Maple 10

I coded it in Classic Worksheet Maple 10

Page 1 of 1