Dear all,

restart:with(plots):
eq1:=diff(f(y), y$4)-(diff(f(y), y$2));

bcs:=f(h1) = (1/2), f(h2) = -(1/2), (D(f))(h1) = -1, (D(f))(h2) = -1:

h1:= 1+cos(x):h2:=-1-cos(x+g):

db:=eq1,bcs:
d1 := subs(g=1,[db]):
P1:= eval(diff(diff(f(y),y$2)-f(y),y));

for x from 0 to 1 by 0.1 do
F2[x]:=dsolve(d1, numeric,maxmesh=25500,output=listprocedure): 
P2[x]:=subs(F2[x],P1); # subing values into P1 
end do:
Vls:=Vector([seq(P2[x],x=0..1,0.1)]):
XX := `<|>`(`<,>`(seq(x, x = 0..1, 0.1))):
plot(<<XX>|<Vls>>, color=red);

I'm trying to plot P1 vs x but getting empty plot. Please help me out. 

Thanks

restart:

Eq1:=(-1/k)*B*(diff(-f(r),r$3)-(1/(r+k)^2)*diff(f(r),r$1)+(1/(r+k))*diff(f(r),r$2)+(1/(r+k)^2))

-((r+k)/k)*B*(diff(-f(r),r$4)+(2/(r+k)^3)*diff(f(r),r$1)-(2/(r+k)^2)*diff(f(r),r$2)+(1/(r+k)^2)*diff(f(r),r$3)

-(2/(r+k)^3));

bcs:=f(-h)=1/2,(D@@1)(f)(-h)=1,f(h)=-1/2,(D@@1)(f)(h)=1;

dsolve({Eq1,bcs},f(r));

Please have a look.

Thanks

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

 pde:=1/r*diff(r*diff(U(r,z),r),r)+1/b^2*diff(U(r,z),z,z)=0;

ans := pdsolve(pde);

ics:=(D[1](U))(0, z) = 0,(D[2](U))(r, 0)-B*U(r,0) = 0;

bcs:=(D[2](U))(r, 1)+B*U(r,1) = B,(D[1](U))(1, z)+B*U(1,z) = 0;

B:=1:b:=1:

S:= pdsolve(pde, {bcs, ics}, numeric);

Error, (in pdsolve/numeric) unable to handle elliptic PDEs

anyway around this?

RIZPDE.mw

restart:
Eq1:=r^2*diff(w(r),r$2)+r*diff(w(r),r$1)-r^2*G5*w(r)-P*r^3*G6-G7*r^3-G8*r^5+C1*r/G3=0;

res1:=dsolve(Eq1);

bcs:=D(w)(0)=0,w(h)=0;

res2:=(dsolve({Eq1,bcs},w(r)));

match(rhs(res2)=rhs(res1),r,s);

eval(_C2, s);

Any idea?

Thanks

Hello everyone,

I came across an image/photo and thought, It will be fun to try it in maple.

Except plotting a few triangles and circles, I couldn't make it. 

Here is the image. 

Have a look please.

Cheers!