Hey people,

I am trying to get the following code to run, but it keeps returning an error about too many arguments

restart:with(plots):with(PDEtools):

alpha_const := 0.5:gamma_const := 2.5: D1 := 0.05: D2 := 0.002:

A := diff_table(a(x,t)):B := diff_table(b(x,t)):

Selkov[1] := A[t] = 1 - A[]*B[]^(gamma_const) + D1*A[x,x]:

Selkov[2] := B[t] = alpha_const * ( A[]*B[]^(gamma_const) - B[]) + D2 * B[x,x]:

bc[1] := D[1](a)(0,t)=0: bc[2] := D[1](b)(0,t)=0: bc[3] := D[1](a)(4*Pi,t) = 0: bc[4]:=D[1](b)(4*Pi,t)=0:

ic[1] := eval(A[],t=0)=a_0:ic[2] := eval(B[],t=0)=b_0:

case1 := eval(ic,[a_0=1,b_0=1]):

case2 := eval(ic, [a_0=piecewise((x<2*Pi+1) and (x>2*Pi-1), 0.99, 1), b_0=piecewise((x<2*Pi+1) and (x>2*Pi-1), 0.99, 1)]):

Case1Default := pdsolve({Selkov[1],Selkov[2]},{bc[1],bc[2],bc[3],bc[4],case1[1],case1[2]},numerical);

Error, (in pdsolve/sys) too many arguments; some or all of the following are wrong: [{a(x, t), b(x, t)}, {a(x, 0) = 1, b(x, 0) = 1, (D[1](a))(0, t) = 0, (D[1](a))(4*Pi, t) = 0, (D[1](b))(0, t) = 0, (D[1](b))(4*Pi, t) = 0}, numerical]

**My code worked just earlier today, and now it wont. **If i try to run pdsolve({Selkov[1],Selkov[2]}) it says that there is an error with general case of floats.

You help is greatly appreciated!