@mmcdara  in the article they got the results in terms of bessels

@mmcdara   thanks 

yes u is a function of r and z, whre P is dp/dz and h(z) =1-3*delta*(11*(z-d)), where delta and d is constants , r  is a function of z.

this equations come from the article " https://www.banglajol.info/index.php/JNAME/article/view/31165"

i want to solve and plot it completely using maple,

i

@Carl Love  thanks

@tomleslie 

thanks a lot

@Carl Love still error in plot when subtracting

plot(exact-trial,x=0..1); 

@Fabio92 

if i am plotting the graph

plot(exact-trial,x=0..1);  

i am getting an error

Error, invalid input: plot expects its 1st argument, p, to be of type {set, array, list, rtable, algebraic, procedure, And(`module`, appliable)}, but received y(x)-(419847598778369106465896716800/5083592179-(1141264098481404792104878055040/5083592179)*exp(-1))*x-(3461556549687386008320/86421067043-(9409486267668445136640/86421067043)*exp(-1))*x^2-(297392012620392925990857231360/86421067043-(808395303834877146808816851840/86421067043)*exp(-1))*x^3-(183196962307265652345600/86421067043-(497980973668785808481280/86421067043)*exp(-1))*x^4-(3716895819325875576015225600/86421067043-(10103570363948922189840551040/86421067043)*exp(-1))*x^5-(80814587...
 

@tomleslie sir  in the piecewise graph, the red colour line start at 0and end at 0.2 , the next should start from 0.2 end at 0.4,but  in between 0 to 0.2 another graph starting at 0.1 and end at 0.3 , i dont require , see the sample (Note) from 0.8 to 1.0 the line start at 0  up till 0.8 then it take the shape,  see the dark black line

@tomleslie 

sir i want the polynomial in the interval [0,1] by using the formula K[i]:=(x-x[i+1])/(x[i]-x[i+1]); by taking 

N:=10;h:=1/N:

and if we add N1+N2.....+N10 we get 1, the sample codes are here
 X:=i->i.h:
X(0);X(1);X(2);X(3);X(4):
N:=(x,i)->piecewise((x<X(i-1),0,x<X(i),(x-X(i-1))/h,x<X(i+1),(X(1+i)-x)/h),0);

plot([N(x, 1), N(x, 2), N(x, 3)], x = 0 .. 1, legend = [N[1], N[2], N[3]])

@Carl Love 

thanks sir, 

I want to solve the second order differential equation using galerkin finite element method and calculate the error also , and plot the graph of direct solution and galerkin FEM solution using maple , code please give some suggestions

@Carl Love thanks

@Kitonum sir how can wedone manually  in legend line + symbol etc...  

@Kitonum 

Thanks a lot but in the legend i am not getting symbol and line i am getting only line , sample graph enclosed

@Kitonum 

Sir i am  getting the error like this

Error, (in plot) expecting option style to be of type identical("point", "patch", "patchnogrid", "line", "hidden", "wireframe", "contour", "patchcontour", "polygonoutline", "polygon", "surface", "surfacecontour", "surfacewireframe", "wireframeopaque", "default") but received pointline
Error, (in plots:-display) expecting plot structure but received: A
 

@Adri van der Meer  sir did we get a linestyle like this 

i.e line diamond symbol, line square, line astrick etc.....

@Adri van der Meer thanks a lot