I can't make sense of this question.  It appears to be referring to the equation

diff(h(x,y,z),z) = I*(A*f(x,y,z) + B*g(x,y,z));

then it evaluates the equation at z=R where R is some constant.
In Maple's notation that is written as

D[3](h)(x,y,R) = I*(A*f(x,y,R) + B*g(x,y,R));

But since the left-hand side has been evaluated at z=R, we no longer have a differential equation.  You need to be specific in posing your question.

To make this post (and your other posts) useful, I suggest that you begin with stating the purpose of your code.  And at the end give an example of its use.  As is, it is difficult to understand the significance of what you have done.

From the post's title it appears that your code is related to the Euler-Bernouli beam equation.  But that equation is simply d^4 u / dx^4 = f(x),  the solution to which is a matter of integrating f(x) four times and applying the boundary conditions. This can be done in a single line by calling Maple's dsolve function. Does your code provide anything beyond that?

@MapleUser2017 You are doing 
LibraryTools:-Save(Mylib, LibLocation);

That should be

LibraryTools:-Save(MyMat, LibLocation);

Show what you have done so far, and if you have specific questions, post it here to get help.

You don't want someone else to do your entire homework for you, do you?

That's an odd set of differential equations.  Are you sure you are not confusing a differential equation with its solution?

If we take the equations that you have given literally, then asking for a phase portrait is not a good question.  A phase portrait is meaningful concept for an autonomous system, that is, a system of differential equation in which the time variable t does not appear explicitly.  For instance,
x1'(t) = 3*x1(t) + 2*x2(t),
x2'(t) = 5*x1(t) + 1*x2(t),
is an autonomous system, but yours has sin(t) and cos(t), and therefore it isn't autonomous.

@ecterrab Things work fine now.  Thank's for the quick fix!

@ecterrab 

Edgardo, I installed the latest Physics update.  The problem with the `+` is now resolved but `-` needs your attention.

Download mw.mw

@vv Thanks for the suggested workaround.  I would have used it, but with Edgardo's latest fix that won't be necessary.

Edgardo, thank you very much for your lighteningly fast remedy.  I will download and install the update now.

And while we are on the subject of Physics[Vectors[, I have a question for you.

The standard orthonormal basis vectors are _i, _j, _k.  Is it possible to introduce another triple of vectors, let's say _a, _b, _c, and tell Maple that these are also orthonormal?  That is, we want Maple to automatically apply the rules  _a . _a = 1,   _a ._b = 0,   _a &x _b = _c,  etc, whenever these vectors are encountered.

I recall having read reports on MaplePrimes of instances where Maple produces corrupt files when the file name contains Danish characters. I suggest limiting file names to ASCII characters to avoid the issue.

@amrramadaneg I see no connection between plotting fields and plotting contours.

@tomleslie He is looking for the equilibrium (that is, steady state) solutions of the heat equation, in which case the time derivative is zero and the heat equation reduces to an ODE.  His shows (correctly) how one solves the ODE by hand. He want to see how to do that step-by-step in Maple.  But the steps are so trivial that bringing in Maple will obfuscate rather than illuminate.

@TeyhaNeedHelp Neither the solution nor its derivative match the figure that you have shown. I suspect that you have entered the equation or the boundary conditions incorrectly.  Check!

Download mw.mw

@Kitonum It looks like the graphs are of the derivatve of f, not f itself.

@Christian Wolinski Another way of getting around the bug is to specify the range=... option to dsolve, as in:

dsol := dsolve({de1,de2,eq,ic}, numeric, output=operator, range=0..10);