## How to create color bar in Maple 13...

I use the following plot in Maple 13. I would like to have the colorbar also along with the plots. Help me to write the commands to get the colorbar.

## Maple 13 also a derivative of the third derivative...

restart:
> y0:=0;
> D[0](y):=1;
> D(D(y))[0]:=A;
> G(x):=diff(y(x),[x\$3])=-1/2*y(x)*diff(y(x),[x\$2]);
> for k from 1 to 12 do
> G(k):=D@@(k)(G)[0];
> value(G(k));
> od;
> for k from 1 to 12 do
> D@@(k)(y):=G[k];
> od;
> for k from 0 to 12 do
> y0:=y0+(G[k]/k!)*x^k;
> od;
> y(x):=y0;

## Evaluate and print values...

how to evaulate the value of R(z) for different values of z=0 to 1 with an interval of 0.1 and print ten values  in one column

R(z):= 1-cos^2*(Pi*z);

## maple take long time...

why i can not evaluate 29 polynomial. maple try to evaluate last 7hr, how many time required too solve it?

## How to plot graphs with different parameters chang...

hai,
greetings,

I require help to plot graphs by changing differernt parameters. i am enclosing my codes and sample codes,

nanofluid_sample.mw

## Getting error in plotting graph...

Hai any one help me to plot the curves and remove the errors  i am attaching the codes and  sample graph. thanks in advanced.

test.mw

## Unable to get output...

I am unable to the get the output  in ans2 , error is comming

restart:
with(DETools):
with(PDEtools):
u[o](r,z):=(-1/4)*diff(p[o](z),z)*(1-r^2):
ode:=gamma1*diff(u[o](r,z),z)+(1/r)*diff(v[1](r)*r,r)=0:
#dsolve( (ode), { v[1](r) } ):
IC1 := {v[1](0) = 0}:
ans2 := combine(dsolve(`union`(ode, IC1),{v[1](r)}));

## display pde in desired form...

How to display the result in the desired form in partial differential equation and also collect the same terms

i am attaching sample codes

restart:
with(DETools):
P[o](z):=C[o]*exp(lambda*z)+D[o]*exp(-lambda*z):
u[o](r,z):=(1-r^2)*diff(P[o](z),z):
v[o](r,z):=(2*r-r^3)*diff(P[o](z),z,z):
A1:=P[o](z)*u[o](r,z):
# we want the result out put  like this
A1 := lambda*(C[o]^2*exp(2*lambda*z)-D[o]^2*exp(-2*lambda*z))*(1-r^2):
A2:=P[o](z)*(1+r^2)+v[o](r,z)+A1:
#we want the  resultant output like this
A2:=lambda^2*(C[o]*exp(lambda*z)+D[o]*exp(-lambda*z))*(1+r^2+2*r-r^3)+A1:

## Differential equations solving ...

Any one can help me to solve the differential equations using maple to get the velocities u ,v and pressure p for the problem mentioned below

## How Do I Save This As My Default Font?...

How do I save these fonts as my default ones?

"Maple Input: Arial 24, Red"

and "2D Output: Arial 18, Blue".

I am using Maple 13.

Thank you!

## How to convert Matlab codes to Maple code...

Any one can help me to convert matlab codes to maple.Matlab_code.mw

 > %Finite element method code for solving bvp nonlinear ODEs%
 > % u''+uu'-u=exp(2x), u(0)= 1, u(1)=e     %
 > %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 > function FEM_Code()
 > clear all; close all; clc
 > n=5;                     % NO of element
 > nn=n+1;                  % No of nodes
 > lgth=1;                  % Domain length
 > he=lgth/n;               % lenth of each elemnet
 > x=[0:he:lgth];           % Data point for independant variable
 > AC=0.00005;              % Accuracy
 > F=zeros(nn,1);           % Initialization
 > F(1)=exp(0); F(nn)=exp(1);  % Boundary conditions
 >
 > %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 > % Direct Iterative process to handle nonlinear problem
 > %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 > c=1.0;
 > count=0;                   % Initializations for count for iterations
 > tic                        % Time start
 > while (c>0)
 > [F1]=assembly(F,n,he);
 > c=0.0;
 > for i=1:nn
 > if (abs(F(i)-F1(i))>AC)
 > c=c+1;
 > break;
 > end
 > end
 > F=F1;
 > count=count+1;
 > end
 > disp('Hence solution=:');
 > %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 > % Output for prinmary and secondary variables %%%
 > %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 > diff=abs(F-exp(x)');
 > fprintf('No of element=%d\n',n)
 > disp('      x       FEM          Exact       Error')
 > disp([x',F,exp(x)',diff])
 > fprintf('No of iterations=%d\n',count)
 > %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 > %%% Ploting of primary variable %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 > %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 >
 > plot(x,F,'--rs','Linewidth',2)
 > xlabel('x')
 > ylabel('u(x)')
 > title('solution plot to given BVP')
 > toc                                % given totlal time
 >
 > end
 > %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 > %%% Derivative of element matrix and Assembly%%%%%%%%%%%%%%%%%%%%
 > %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 >
 > function [F1]=assembly(F,n,he)
 > nn=n+1;
 > k = zeros(nn,nn);               % Initialization of main Matrix
 > R = zeros(nn,1);                % Initialization of RHS Matrix
 > syms x                          % x as symbolic variable
 > s=[1-x/he,x/he];                % linear shape function
 > ds=diff(s,x);                   % Differentiations of shape function
 > lmm =[];
 > for i=1:n
 > lmm=[lmm;[i,i+1]];          % connectvity Matrix
 > end
 > for i=1:n
 > lm=lmm(i,:);
 > %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 > %%% Generation of Element Matrix k11 and RHS Matrix f1%%%%%%%%%%
 > %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 > k11=-int(ds'*ds,x,0,he)+(int(s'*ds*s(1),x,0,he)*F(lm(1))...
 > +int(s'*ds*s(2),x,0,he)*F(lm(2)))-int(s'*s,x,0,he);
 > f1 = int(exp(2*(x+(i-1)*he))*s',x,0,he);
 > %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 > %%% Assembly accroding to connectivity Matrix%%%%%%%%%%%%%%%%%%%%
 > %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 >
 > k(lm,lm) = k(lm,lm) + k11;
 > R(lm) = R(lm) + f1;
 > end
 > %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 > %%% Imposing Boundary Conditions %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 > %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 >
 > k(1,:) = 0.0; k(nn,:) = 0.0;
 > k(1,1) = 1.0; k(nn,nn) = 1.0;
 > R(1,1) = F(1); R(nn,1) = F(nn);
 > %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 > %%% Solution of equations (F1) %%%%%%
 > %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 >
 > d = k\R;           % better than using inverse k*R
 > F1 = d;
 > end
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >

i am attaching the codes and file. i will be thankful to you .

## How to activate maple...

How can I activate my maple 13?

I have no activatIon left but installed maple

once

My purchase Code is BNN7VMT4BPFN7ESG

Thank. You

## Error with Procedure output...

 >
 >
 (1)
 >
 >
 >
 (2)
 >
 (3)
 >

Way i get the wrong output with procedure S1(uu)

## Arrow plot color error...

arrow([u1, u2, u3], shape = cylindrical_arrow, length = sqrt(LinearAlgebra:-Transpose(f0).f0), color = [red, black, green]):

Why i get this error

Error: the number of colours in the COLOUR structure is not sufficient.Error: the number of colours in the COLOUR structure is not sufficient.Error: the number of colours in the COLOUR structure is not sufficient.Error: the number of colours in the COLOUR structure is not sufficient.Error: the number of colours in the COLOUR structure is not sufficient.Error: the number of colours in the COLOUR structure is not sufficient.