Question: Problem with for loop

Hello every one,

I am facing problem with for loop. Please have a look.

My code is

restart;with(PDEtools):

Eq1:=diff(u(eta,t),t)=C*diff(u(eta,t),eta$2):

BCs := {u(0,t)=sin(t), u(L,t)=0}:

ICs := {u(eta,0)=0}:

L:=600:

sol:= pdsolve(Eq1,ICs union BCs,numeric,time=t,range=0..L);

for i from 0 to 600 do 

for t from 0 to 3 by 1 do

for C from 0 to 1 by 0.1 do

p0[C,t]:=sol:-value (C,t):

A[i*0.01]:=p0[C,t](i*0.01) 

end do:

end do:

od:

for i from 0 to 600 do 
u0[i]:=subs(A[i*0.01], u(eta,t)) 
od:

u0:=<<(seq(u0[i], i=0..600 ))>>;

 

Thanks

Please Wait...