> restart; with(plots); with(LinearAlgebra); with(IntegrationTools); pde := rho(x)*(Diff(u(t, x), t, t))+lambda(x)*(Diff(u(t, x), t))-(Diff(P(x)*(Diff(u(t, x), x)), x))-Q(x)*u(t, x)+R(x);

> eq1 := int(pde*phi[i](x), x = -1 .. 1);

eq2 := expand(eq1);

> eq3 := applyop(u -> Parts(u, phi[i](x)) end proc, 3, eq2);

i get an error

Error, (in applyop) improper op or subscript selector

how to correct it??

Good morning sir.

I request your kind support to the above cited question.

With thanks & Regards

M.Anand

Assistant Professor in Mathematics

SR International Institute of Technology,

Hyderabad, Andhra Pradesh, INDIA.

Just came across this

Press esc twice then enter.  Maple becomes unresponsive. 

Hi -

I'm trying to use Maple to integrate an ODE numerically, but to output the integration steps as I go (in an array, or similar). None of the dsolve/numeric output parameters appear to offer this (ie

x=0.00 f = whatever,

x = 0.05 f = whatever2,

x = 0.10 f = whatever3)

Is there a standard feature in Maple that I'm missing, or would I have to force the output from the dsolve procedure in some way, please?

i have a piecewise like " g "

g := piecewise(0 <= t and t <= 0.2e-1, 122.63, 0.2e-1 < t and t < .287, 0, .287 <= t and t <= .307, 122.63, .307 < t and t < .335, 0); plot(g, t = 0 .. .335)

this function is semi periodic just for 10 times consecutive . i want to define it

example     k=[g g g g g g g g g g  g ] or somthing like that.

1.mw 

i want to plot fft form of final piecewise function

and then filter it . and plot them.

Hi there..

I have a question on how to do a pointplot.

before plotting, I need to know the value of lambda[j] and  all the values of lambda already have.

so now I need to plot a graph with the values of lambda with different range and different colour,

Let say I have

> for j from 17 to 32 do k[j] := j+1;

x[j] := add(P[j, 1], j = j-1 .. j+2);

X[j] := add(P[j, 1]^2, j = j-1 .. j+2);

y[j] := add(P[j, 2], j = j-1 .. j+2);

Y[j] := add(P[j, 2]^2, j = j-1 .. j+2);

xy[j] := add(P[j, 1]*P[j, 2], j = j-1 .. j+2);

cx[j] := evalf(x[j]/k[j]);

cy[j] := evalf(y[j]/k[j]);

c11[j] := evalf(X[j]/k[j]-cx[j]^2);

c22[j] := evalf(Y[j]/k[j]-cy[j]^2);

c12[j] := evalf(xy[j]/k[j]-cx[j]*cy[j]);

C[j] := evalf(Matrix(2, 2, [[c11[j], c12[j]], [c12[j], c22[j]]]));

E[j] := simplify(fnormal(LinearAlgebra[Eigenvalues](C[j])));

if E[j][1] > E[j][2] then lambda[j] := E[j][2]/(E[j][1]+E[j][2]) else lambda[j] := E[j][1]/(E[j][1]+E[j][2])  end if;

lambda[j];

 end do;

the range of lambda [j] are as follows:

 0.02< lambda [j]<0.06

 0.06< lambda [j]<0.12

 0.12< lambda [j]<0.18

for i from 17 to 32, do if   0.02< lambda [j]<0.06 then green[i]:=P[i,j]; j:=i+1 elif

 0.06< lambda [j]<0.12 then red[i]:=P[i,j];j:=i+1 ; elif 0.12< lambda [j]<0.18 then blue[i]:=P[i,j];j:=i+1 end if;end do

how to do a point plot with the above situation so that in my plotting all the information are on the same graph.

All help is greatly appreciated.

Thanks

Hello.
I'm new to Maple, so my question could appear idiotic. But, anyway. :-)

I want to understand the restrictions of code generation in Maple.
What functions could I convert to code?
Only listed in "Translation of Functions" block in this link?
http://www.maplesoft.com/support/help/MapleSim/view.aspx?path=CodeGeneration/General/TranslationDetails

For example, I've found this post - http://www.mapleprimes.com/posts/36616-Dynamic-Programming-In-Maple 

Could functions "Maximize" and "diff" be converted to Java, VB or something or not? If yes, how?

Thanks anyway. :-)

Im trying to plot the different values of shanks as points in this ] loop.

pp:=0.0:sumn:=0:sumnm1:=0:sumnp1:=0:

for i from 1 to 60 do

if(n>2) then sumnm1:=sumn fi;

if (n>1) then sumn:=sumnp1 fi; p:=1/evalf(2*beta*BesselJ(2,zeros[i])+epsilon*zeros[i]*BesselJ(0,zeros[i])+zeros[i]*BesselJ(1,zeros[i])); pp:=pp+p;

sumnp1:=pp;

if (n>2) then shanks:=(sumnp1*sumnm1-sumn*sumn)/(sumnp1-2*sumn+sumnm1) fi;

shanks;

od;

sumnp1: -2*sumn+sumnm1:

I have tried plotting it but it only plots the last value. 

Many thanks James

I recently installed Maple 17 from my university and greek letters do not work. For instance if I type alpha or any other greek letter, Maple just outputs a capital A in a box. I still have Maple 15 installed and it works on there and all the preferences seem the same.

Thanks

I have 2nd order nonlinear ode I try to solve with Runge Kutta 4th order method in maple but all I get from the outcome was 1 and 0.This is the equation:theta_ode1.mw . How do I do it Or how do I write the code to solve it with maple using  Runge Kutta 4th order method?

Dear Mapleprimes members,

I am trying to understand some known classical theorems of geometry (thales intersection thm, ceva, menelaus ) from the perspective of affine geometry and projective geometry.(my paper/pencil work is ok.. But I would like to use  Maple to visualize -if possible - geometrically some theorems and properties in a graph)

How could I do it in maple? Where may I get some examples?

thank you,

JeanJ

I have 2nd order nonlinear ode I try to solve with Runge Kutta 4th order method in maple but all I get from the out is 1 and 0.This is the equation: theta_ode.mw . How do I do it Or how do I write the code to solve it with maple using  Runge Kutta 4th order method?

Need help with starting this question. Thanks!

If a sequence is defined by X₀=0, X₁=1, X₂=2 and X_n=n(X_n-3+X_n-2+X_n-1), n>=3

How many digits will X₂₀₁₃ have?

I have no idea how to start to answer this question. Please help Thanks!

Let G(x) denote the no. of ways of representing n as the sum of two prime numbers.

For eg, G(10)=2: where 10=3+7=5+5

           G(20)=2: where 20=3+17=7+13

           G(30)=3: where 30=7+23=11+19=13+17

How do you find G(10 000 000)?

I assign a function like so: f:=x->0.2*x^2*(x-3)^3

Then using d/dx in the expression palette I differentiate: d/dx f(x)

Get the following: 0.6x^2 (x-3)^2 D(x^2)(x-3)

What is the "D" in the last expression? What hasn't the whole differentiation been completed?