Okay, here goes I am have some questions. I am solving the Mathieu differential equation whit a damping term, so the equation I am solving is: u''(x)+b*u'(x)+(a-q*cos(\Omega t))=0 Where I have a=0 and b=1,9, and whit the initial conditions u(0)=0 and u'(x)=1, this gives me the following: u(x) = exp(-9/10*x)*MathieuS(-81/100, 1/2*q, x) I would like to see for which values of q, that the function u(x) is in the interval [-20;20]. I can solve the problem by drawing some plots, but I still haven’t been able to solve it without. I know this is a long short, since you guys have so little info on the problem at hand, but any help will be much appreciates.

So I have a vector X that contains my least squares coefficients. How do I incorporate these coefficients back into a polynomial function automatically using loops or the like? Right now I have to define the funciton manually f:=x-> X[1]*x^3 + X[2]*x^2..... but I want to be able to quickly vary the degree of the polynomial without having to edit the function definition by hand each time. Any geniuses out there that can help me? Thanks, Matt

If f(x)=-4x^(2) and g(x)=2/x, find [g*f](x).

Implement a computer programme which, for a given function f(x), interval [a,b] and real numbers epsilon1>epsilon2>0, finds an approximation x~ from the inerval (a,b) with an error less than epsilon1 by bisection and, with x~ as the stating value, finds an approximation of this root by the second modification of tj=he Newton method with an error less than epsilon2.

Download 285_quantmechprob1.14.mws
View file details I'm trying to perform an integration on a function involving a couple constants in order to normalize (setting result of integral = 1) and solve for the constant, A. The closest I'm able to reach with respect to completing the integral is a result that involves the error function, erf. How do I go about getting a neater closed form solution to the this integral ? Am I illegal integration limits ?

In a worksheet I have the following subs({tan(a)=sin(a)/cos(a),tan(b)=sin(b)/cos(b)}, eq1) The following doesn't work, but I was wondering if there was any wild card character that would enable me to do something like this subs(tan(?)=sin(?)/cos(?), eq1)

I have been trying to come up w/ a fast way to get identity Matrices which can later accept non-zero non-diagonal entries. If I use the shape=identity parameter then the Matrix can't be manipulated. This is what I've come up with KronDelt :=proc(i,j) => if i=j then 1 else 0 end if end proc: Matrix(1..4,1..4,(i,j) -> KronDelt(i,j)) Was it necessary to seperately define KronDelt ? I tried putting a mapping directly into the Matrix constructor but always got a syntax error, e. g., Matrix(1..4,1..4,(i,j)-> if i=j then 1 else 0 end if ) This doesn't work. Is there any way to d

I am having the hardest time plotting the qualitative behavior of the solutions of these differntial equations. I keep getting Error, (in plots/animate) no non-zero vectors found. If someone could walk me through plotting these equations, it would be greatly appreciated. Equation1: dx/dt=x^2, x(0)=1, 0≤t<>

My goal is to obtain a formula F:= x -> fa(x)*c1_1(x)+ fb(x)*c1_2(x)+fc(x)*c2_1(x)+... where every ci_j(x) is a function of x. fa, fb, fc are known functions of x. I do not need to have it printed, I just need it to return a numerical value for every value of x I throw in. I have obtained the solutions of ci_j's in the form Sol[1]:=[c1_1=f11(x), c1_2=f12(x), c1_3=f13(x), ...], N1 terms Sol[2]:=[c2_1=f21(x,c1_j’), c2_2=f22(x,c1_j’), c2_3=f23(x,c1_j’), ...], N2 terms Sol[3]:=[c3_1=f31(x,c1_j’,c2_j’), c3_2=f32(x,c1_j’,c2_j’), c3_3=f33(x,c1_j’,c2_j’), ...], N3 terms Sol[4]:=[c4_1=f41(x,c1_j’,c2_j’,c3_j’), c4_2=f42(x,c1_j’,c2_j’,c3_j’), c4_3=f43(x,c1_j’,c2_j’,c3_j’), ...], N4 terms

For testing an external routine against Maple I have something like

  Gamma_mpl:=define_external(
    'Gamma_mpl',
    'C',
    'Z'::ARRAY(1..n,datatype=float[8]),
    'result'::ARRAY(1..n,datatype=float[8]),
    LIB=theDLL);

where the function in the DLL (compiled with Dev-C++) is like

  extern "C" __declspec( dllexport ) __stdcall 
  void Gamma_mpl(double *Z, double *R){ ... some code  ...}

It takes the input in Z and updates values in R.

The following works as expected:

  z:= -1.1 + 2.1*I*0:
  # provide memory
  Z:=Array(1..2, order= Fortran_order, datatype= float[8]):

How do I instruct a worksheet to automatically store a sequence of plots as .gif images? In my worksheet, I generate a sequence of say 50 plots and then look at them within the worksheet. I really need to store these someplace else on my computer as .gif files and do not want to export each plot individually.

I need to get the slope of this plotted line that I did, and I have no clue how. Here is the points I used: points := [[22, 100.19], [40, 200.19], [50, 300.19], [62, 400.19], [75, 500.19], [95, 600.19]] I'm not sure what do to, I tried slope[[22,100.19],[95,600.19]] so try and get the slope of the line between those two points, but nothing came up, it just went to the next line and said "slope".

I want to write a maple program/script involving matrix computations. The program will consist of approximently 50 statements and will require user supplied input data. I have Maple 9 and would like to call the script from Matlab. Matlab would provide the input data to the Maple script and then I would like to have Maple pass the output data to Matlab. Is there a way to do this using Maple 9 or do I need the new Maple Matlab toolbox? Thank you, Richard

Hello together, I've just started to do some calculations with maple 9.5 and still in trouble. I would like to modify the step function below to a n-periodic function. In maple version 10 there's a function called "rept" which is not available in my version unfortunately. Does anyone have a ease of use solution? f:=(MTBF,MTTR,t,n) -> piecewise(t<>

Hi, I have a small problem. I want to use the Schurform computation in the LinearAlgebra-package. The main objective is to compute the stabilizing solution to a Algebraic Matrix Riccati equation. (perhaps someone knows some code for this?) I call the SchurForm function as (T,Z) := SchurForm(A,output=['T','Z']) where A has the correct type (ie complex(sfloat)). This works well in a Maple sheet, but in a procedure (where T and Z are defined local) the kernel says that the output chosen is invalid. If I don't specify any output, it also works good (but then I get no Z). The Maple version is 10 and I am using Linux.