I am solving a PDE coupled with an ODE, as following:

pde1:=diff(RR(tau), tau) = dRdt                                                                               (1)

Anyone know how to debug FOR loops in MAPLE? Thanks in advance.

Luciane.

I am trying to solve a problem of liquid-liquid equilibrium based on the gibbs free energy equation:

DGdef := DG/RT = sum(n[i]*ln(phi[i]), i = 1 .. nc)+sum(sum(X[i, j]*n[i]*phi[j], j = i+1 .. nc), i = 1 .. nc)

DGdef1 := subs(n[i] = phi[i]*V/v[i], DGdef)

DG := eval(rhs(DGdef1), nc = 3)

Hallo Everybody,

I hope someone can help me.

I want to solve numerically the diffusion equation,

          diff(csh(rh, t), t) = Dif*(diff(rh^2*(diff(csh(rh, t), rh)), rh))/(r0*rh^2)

          where "Dif" is the diffusivity and "r0" is the initial value of rh (the equations were normalized).

subject to the following boundary and initial conditions: