agim

0 Reputation

2 Badges

15 years, 310 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