Where u and v are the displacement components in x and y directions respectively.

Hi,

I'm having trouble of converting my y-axis to decibel.  Here is is my script:

restart;
with(plots);
loglogplot(1/(1+f), f = 10^(-4) .. 10^4);

Everything is fine except for the y-axis, which is not in dB.

I have been having problems with using the BodePlot function with units:

R1 := 18.2*10^3*Unit('Omega');

R2 := 10^3*Unit('Omega');

C1 := 470*10^(-12)*Unit('F');

C2 := 4.7*10^(-9)*Unit('F');

# wo is in hertz

wo := 1/sqrt(R1*R2*C1*C2);

# Q is unitless

Q := wo*R1*R2*C2/(R1+R2)

with(DynamicSystems);

sys := TransferFunction(wo^2/(s^2+wo*s/Q+wo^2));

This is the error message I got:

Error, (in Units:-Standard:-+) the units `1` and `Hz` have incompatible dimensions

I think the problem is that the BodePlot function doesn't expect 'wo' to have units.  

So I tried to work around the issue by using the loglogplot but it doesn't seem to like 

complex function even when I used abs to find the magnitude (with or without units).

 Any workaround is appreciated.

I am trying to simplify the eigenvalues of a 2x2 matrix [[a,b],[c,d]] subject to the condition a,b,c, and d are integers such that a+b = c+d. Why do the following commands not achieve this?

with(LinearAlgebra):

A:=Matrix[[a,b],[c,d]]):

Eigenvalues(A) assuming a::integer,b::integer,c::integer,d::integer,a+b=d+c

How might I achieve what I need?

I would be grateful if anyone would comment on the compatibility of Maple 18 and MacOS 10.10 (Yosemite).

The ability of Maple to solve differential equations is unsurpassed, but when the solutions appear in terms of Heun functions that result is disappointing because it is either difficult or impossible to convert those functions to other functions more commonly used and for which plots are readily generated.

Specifically, does any reader have a suggestion what to do with Heun C and Heun G functions?  In principle, they seem to be related to 1F1 and 2F1 hypergeometric functions, but the conversion seems not to succeed, and it is not obvious how to make it succeed.  In both cases of interest, the literature contains hints of solutions in other functions.

It seems that a solution of a differential equation in terms of Heun functions is not a solution at all.

Today science professionals in engineering software used to only work on the desktop and even just looking to download and use mobile apps math; but they are not able to design their own applications.Maplesoft to set the solution to it through its Maple package; software supports desktop and mobile; solves problems of analysis and calculation with Embedded Components. To show this we have taken the area of different mathematical topics; fixed horizontally to a certain range of parameters and not just a constant as it is customary to develop. This paper shows how the Embedded Components allow us to develop mathematics in all areas. Achieving build applications that are interactive in mobile devices such as tablets; which are used at any time. Maple gives us design according to our university or research need, based on contemporary and modern mathematics.With this method we encourage students, teachers and researchers to use graphics algorithms.

CSMP_PUCP_2014.pdf

Coloquio_PUCP.mw

Lenin Araujo Castillo

Physics Pure

Computer Science

Dear all, I have been trying to use Runge-Kutta method to plot an approximate solution with the following code. However, although I can get the numerical approximation the plot would not show.

h := .1;

x[0] := 0;

y[0] := 1;

xf := 3;

n := floor(xf/h)

f:= (x,y)->1/(3 y-x-2)

x := x[0]

y := y[0]

for i to n do

k1 := f(x, y);

k2 := f(x+(1/2)*h, y+(1/2)*h*k1);

k3 := f(x+(1/2)*h, y+(1/2)*h*k2);

k4 := f(x+h, h*k3+y);

k := (k1+2*k2+2*k3+k4)*(1/6);

y := h*k+y;

x := x+h

end do;

y[n]

data := [seq([x[n], y[n]], n = 0 .. 30)];

p[2] := plot(data, style = point, color = blue);

p[3] := plot(data, style = line, color = blue);
display(seq(p[n], n = 2 .. 3));

Hello,

Is it possible in Maple to select summand with some command?

For example:

eq := a+b = 0;
                           a + b = 0
eq-b;
                             a = -b

How to generalize second line?(in this case I explicitly told to subtract b, but I want to pick b by some command)

And I am not talking about solving, isolation or anything like that in this case. I just want to transfer something from lhs to rhs by selecting it beforehand.

Thank you.

Hello,

How can one simplify following expression:

After applying 'simplify' command I am getting this:

Powers are not distributed between bases.

How to force Maple simplify it further to

Thank you.

I have about a year of experience in Maple, but as we know, we often come accross very simple problems that are difficult to find solutions to. Part of the problem is that I have not used conditional tests within for loops before.

I have a very large matrix, but until I solidify the algorithm, I'm practicing on a small matrix.

I created a 3 by 3 fibonacci matrix A:=<21,5,1|13,3,1|6,2,0>

I am trying to change the 1s to 0s without specific index reference. Here is my psuedocode that I hope to work on after some tips:

• for i from 1 to 9
  if A[k] = 1 then
  A[k]=0
  end if
  end do

I get an error message saying "unterminated for loop," which is confusing me.

I want to integrate the following cumulative distribution function with variable bounds but maple returns the integral.

>with(Statistics);

>N := RandomVariable(Normal(0, 1)):

>a := (ln(98.53*(1/95))+(0.1e-1+.5*x^2)*.5)/(x*sqrt(.5))

>y:= CDF(N, a, inert = true);

How do i make maple return a function of x instead of an integral?

Hi,

I used to use windows 32 bit, but I have 64 bit windows now. I installed the MapleSim 6.4 and Maple 18 and try to run the simulations that I created with the 32 bit windows. I have this error 'Unable to compile (rc=1), please try again, and if that fails verify your Windows compiller installation'. Could you please let me know what should I do to eliminate this problem?

Best

Onder

Hello,

      I would like to solve a system of 9 nonlinear equations, with the constraints on all 9 variables to be that they are nonnegative. How can I do this?

My code is below - I am trying NLPSolve and have tried solve, but am getting stuck.

with(Optimization);

restart; eq1 := 531062-S/(70*365)-(.187*(1/365))*(H+C+C1+C2)*S/N = 0;eq2 := (4/365*(T+C))*S/N-(.187*(1/365))*(H+C+C1+C2)*T/N-(1/(70*365)+1/(5*365))*T = 0; eq3 := (.187*(1/365))*(H+C+C1+C2)*S/N-(4/365)(T+C)*H/N-(1/(70*365)+1/(4*365))*H = 0; eq4 := (.187*(1/365))*(H+C+C1+C2)*T/N+(4/365*(T+C))*H/N-(1/(70*365)+3/(8*365)+.2*(1/365)+.1)*C = 0; eq5 := .1*C-(1/(70*365)+1/(4*365)+1/60+.5)*C1 = 0; eq6 := (1/60)*C1-(1/(70*365)+1/(4*365)+1/210+.5)*C2 = 0; eq7 := .5*C1-(1/(70*365)+1/60+0.1e-2)*CT1 = 0; eq8 := .5*C2-(1/(70*365)+1/210+(1/9)*(0.1e-2*7))*CT2+(1/60)*CT1 = 0; eq9 := N-S-T-H-C-C1-C2-CT1-CT2 = 0; soln := NLPSolve({eq1, eq2, eq3, eq4, eq5, eq6, eq7, eq8, eq9}, {C, C1, C2, CT1, CT2, H, N, S, T}, assume = nonnegative);

The problem is The square root of 16-x^2 over the interval [0,-4]  0 being the upper bound, -4 being the lower bound.  I have solved 3/4s of this problem but I don't understand what they mean by "Solve the definite integral exactly by geometry".