@AHSAN Even after changing q(y)=j to u(y)=j the problem still makes no sense.  The "bc" and "bc1" in your formulation stand for "boundary condition" but you have not said where the boundaries are.

Go to where you found the problem, read it carefully, then come back with the correct statement.

I looked through your worksheets but was unable to understand them.  It will be very helpful to your readers to have some verbal description of the purpose of the worksheets.  You may want to supply the following information:

1. It appears that you are solving the PDE 
   
   with the initial conditions u=(x,0)=sin(Pi*x), diff(u(x,t),t)=0 on the interval (0.1).  Is that correct?
2. You state the boundary conditions as
   bc1:=u(x,t):        bc2:=u(x,t):
   I don't know that those mean.  Explain.
3. It appears that you discretize the equation through centered finite differencing in the x direction. Is that correct?
4. Then you discretize the equations in the time direction in some way.  That does not seem to be consistent with your stated objective of applying the method of lines, which approaches the problem through solving a system of ODEs, but perhaps I have misunderstood.  Explain.
5. What is the idea behind your non-classical method?
6. You write "I need to compare the result of the method with the exact solution".  What is the exact solution?
7. You must know that you may solve your PDE with Maple's pdsolve() (numerically).  You are aware of that, aren't you?
8. Why have you posted two worksheets?  What is the purpose of each?  Provide a guide.

@janhardo There are neither double integrals nor polar coordinates there.  The area is calculated by simple integration in Cartesian coordinates.

@AHSAN Your bc says q(y) = 0.  Your bc1 says q(y) = j .  Do you want both?

@Fancypants Yes, equation (20) looks good.

@Fancypants Okay, that's better.  Now you have an ODE for S[y](t).  The ODE involves phi(t), psi(t), theta(t), f(t), etc., with various subscripts.  Are those known/given functions?  If so, then your equation (13) is a simple second order linear ODE which is trivial to solve.  On the other hand, if you expect to determine the functions phi(t), psi(t), theta(t), f(t), etc,, then you need to provide additional equations.  In general, you need the number of equations to be equal to the numbers of unknowns.

@Fancypants In the worksheet you say that you expect partial_py1 to have terms involving the derivatives of S[y](t) with respect to py1.  But Maple hasn't been told that S[y](t) depends on py1, therefore it sets all those derivatives to zero.

There may be other errors in the worksheet but I don't quite see the aim of your calculations, and therefore I cannot suggest a remedy.

@spinoza Are you going to pass the solutions that you have obtained here as your own?  Or are you going to tell your teacher that Carl/Kitonum solved them for you?

@Subhan331 I ask again: Why do you need Maple code for an algorithm that produces nonsense?  I wouldn't waste my time on that, and nor should you.

As a student and a beginner, you should spend your time on more productive pursuits.

I have no problem with "with(combinat)" in Maple2017 and Maple2020 on Linux.  I suspect that you have typed something incorrectly.  Post your worksheet so that people can examine the source of the error.

To post a worksheet, reply to this message, and note the big fat green arrow in the toolbar of the window where you compose your reply.  Click on that arrow to upload your worksheet.

@Subhan331 Doing the exercise that I noted takes just a few lines of simple calculations by hand.  What's the point of writing Maple code if the algorithm produces nonsense?

As to fractional order derivatives, isn't is too soon to jump into that?  See if you can make sense of the algorithm applied to traditional PDEs first.

An alternative to dharr's suggestion is to press Ctrl+Shift+G (that is, the Ctrl, Shift, and G keys together) followed by the Greek letter code.  For instance, Ctrl+Shift+G followed by e produces ε, and Ctrl+Shift+G followed by s produces σ.  Most Greek letters map to their corresponding Latin keys.  For the complete table of Greek character, look up the help page 2-D Math Shortcut Keys and Hints.

The Ctrl+Shift+G combination works on Linux and Windows.  On a Mac, type Command+Shift+G instead.

@Preben Alsholm What you have quoted can be an explanation of what we are seeing.  It's not quite clear to me how, but it seems relevant.  Thanks for your input.  Because of that I have converted your Reply to an Answer.

@Pascal4QM You need to say
A,C := test()[[1,3]];