> y:='y';

> eta:=evalf(y/R*((2*vzmax*R^2/(9*Das*z))^(1/3))*cm); Defining eta in terms of y
f:=y->int(exp(-eta1^3),eta1=0...46.0204*y)/Gamma(4/3); When z is constant, concentration depends only on the distance from the wall, or y

Modified yxplot from CENG 303. Maple couldn't handle plotting f in integral form; hence values of y and f(y) were determined and plotted.

> yxplot2:=proc(xn,yn,x,y)
local n,z,k;
n:=nops(x);
z:=[[x[1],y[1]]];
for k from 2 to 100 do z:=[op(z[1..k-1]),[x[k],y[k]]] od;
plot(z, labels=[xn,yn])
end;

> y:=array(1...100);:y[1]:=0:
for i from 2 to 100 do
y[i]:=y[i-1]+.001:
od:

> fy:=array(1...100);

> for i from 1 to 100 do
fy[i]:=evalf(f(y[i]));
od;

> yxplot2('y','f',y,fy);