A graph for an approximation of the initial value

> u0:=(r,z)->sum(sum(A[n,m]*BesselJ(0,mu[m]*r)*sin(lambda[n]*z),n=1..40),m=1..30);

> plot(u0(r,3/4)+550,r=-5..5);

>