## How do I solve a complex equation in function of p...

Guys,

I am not very familiar with Maple and have to solve quite a complex equation.

I have an equation which is complex ,containing I . I split this equation up in Re=0 an Im=0 . I have to get an answer in function of other parameters, in order to plot these... Maybe it s easier if you look at the work sheet

 >
 >
 >
 (1)
 >
 (2)
 >
 (3)
 >
 >
 >
 (4)
 >
 (5)
 >
 (6)
 >
 >
 (7)
 >
 (8)
 >
 (9)
 >
 (10)
 >
 (11)
 >
 (12)
 >
 (13)
 >
 (14)
 >
 (15)
 >
 (16)
 >
 >
 (17)
 >
 (18)
 >
 (19)
 >
 >
 >
 >
 >
 >
 (20)
 >
 (21)
 >
 (22)
 >
 (23)
 >
 (24)
 >
 (25)
 >
 (26)
 >
 (27)
 >
 (28)
 >
 (29)
 >
 (30)
 >
 (31)
 >
 (32)
 >
 >
 (33)
 >
 >
 >
 (34)
 >

## How can I get the Transfer Function from an Equati...

Hi, I am new in Maple. If I have an electric network as in the figure, I want to get the Transfer Function V2(s)/Vi(s) from this equation system:

Vi=R1.I1+V1

Vi=R2.I2+V2

I1-I2=C1.dV1/dt

I2=C2.dV2/dt

Which are the commands that I may write to get this?? Before hand, Thanks by your answers!

I have four matrix equations

P1, P2, P3 are known 4x4 matrix.

A1 A2 A3 A4 are known 1x4 matrix.

x1 x2 x3 x4 are 1x1 known matrix.

U is 4x4 unknown matrix.

These equations are

(A1T*U*P1*A1) +( (P2*A1)T*U*P1*A1) + ( (P3*A1)T*U*A1) + ( ( P3*A 2)T*U*P1*A1) + x1 =0;

(A2T*U*P1*A2) +( (P2*A2)T*U*P1*A2) + ( (P3*A2)T*U*A2) + ( ( P3*A2 )T*U*P1*A2) + x2 =0;

(A3T*U*P1*A3) +( (P2*A3)T*U*P1*A3) + ( (P3*A3)T*U*A3) + ( ( P3*A3 )T*U*P1*A3) + x3 =0;

(A4T*U*P1*A4) +( (P2*A4)T*U*P1*A4) + ( (P3*A4)T*U*A4) + ( ( P3*A4 )T*U*P1*A4) + x4 =0;

How can i find 4x4 matrix U by using these above four equations??

Thank you

## Extract number from equation...

E_T := (2/mu-2/r)*exp(-r/mu)*Pi^2;

How do I extract the numbers out of the equation, so it becomes

## Unable to convert to an explicit first order syste...

 >
 >
 (1)
 >
 (2)
 >
 (3)
 >
 (4)
 >
 (5)
 >
 (6)
 >
 (7)
 >
 Error, (in dsolve/numeric/process_input) system must be entered as a set/list of expressions/equations
 >
 (8)
 >

## How to solve numerically a Fredholm integral equat...

Hi guys,

I am trying to solve a Fredholm equation of the second kind using Maple. An analytical expression cannot be in principle found. I was wondering whether Maple does numerical evaluation of such integral equations. Please see the equation in attach. Any help is highly appreciated.

Thanks

F

Question.mw

## Digitizing mathematics: ODEs, Special Functions...

by: Maple

The material below was presented in the "Semantic Representation of Mathematical Knowledge Workshop", February 3-5, 2016 at the Fields Institute, University of Toronto. It shows the approach I used for “digitizing mathematical knowledge" regarding Differential Equations, Special Functions and Solutions to Einstein's equations. While for these areas using databases of information helps (for example textbooks frequently contain these sort of databases), these are areas that, at the same time, are very suitable for using algorithmic mathematical approaches, that result in much richer mathematics than what can be hard-coded into a database. The material also focuses on an interesting cherry-picked collection of Maple functionality, that I think is beautiful, not well know, and seldom focused inter-related as here.

 Digitizing of special functions, differential equations, and solutions to Einstein’s equations within a computer algebra system   Edgardo S. Cheb-Terrab Physics, Differential Equations and Mathematical Functions, Maplesoft Editor, Computer Physics Communications

 • Big amounts of knowledge available to everybody in local machines or through the internet
 • Take advantage of basic computer functionality, like searching and editing

 • By digitizing mathematical knowledge inside appropriate computational contexts that understand about the topics, one can use the digitized knowledge to automatically generate more and higher level knowledge
 Challenges 1) how to identify, test and organize the key blocks of information,   2) how to access it: the interface,   3) how to mathematically process it to automatically obtain more information on demand
 Three examples

Mathematical Functions

"Mathematical functions, are defined by algebraic expressions. So consider algebraic expressions in general ..."

"Supporting information on definitions, identities, possible simplifications, integral forms, different types of series expansions, and mathematical properties in general"

 Examples
 General description
 References

Differential equation representation for generic nonlinear algebraic expressions - their use

"Compute differential polynomial forms for arbitrary systems of non-polynomial equations ..."

 The Differential Equations representing arbitrary algebraic expresssions
 Deriving knowledge: ODE solving methods
 Extending the mathematical language to include the inverse functions
 Solving non-polynomial algebraic equations by solving polynomial differential equations
 References

Branch Cuts of algebraic expressions

"Algebraically compute, and visualize, the branch cuts of arbitrary mathematical expressions"

 Examples
 References

Algebraic expresssions in terms of specified functions

"A conversion network for arbitrary mathematical expressions, to rewrite them in terms of different functions in flexible ways"

 Examples
 General description
 References

Symbolic differentiation of algebraic expressions

"Perform symbolic differentiation by combining different algebraic techniques, including functions of symbolic sequences and Faà di Bruno's formula"

 Examples
 References

Ordinary Differential Equations

"Beyond the concept of a database, classify an arbitrary ODE and suggest solution methods for it"

 General description
 Examples
 References

Exact Solutions to Einstein's equations

"The authors of "Exact solutions toEinstein's equations" reviewed more than 4,000 papers containing solutions to Einstein’s equations in the general relativity literature, organized the whole material into chapters according to the physical properties of these solutions. These solutions are key in the area of general relativity, are now all digitized and become alive in a worksheet"

The ability to search the database according to the physical properties of the solutions, their classification, or just by parts of keywords (old paradigm) changes the game.

More important, within a computer algebra system this knowledge becomes alive (new paradigm).

 • The solutions are turned active by a simple call to one commend, called the g_  spacetime metric.
 • Everything else gets automatically derived and set on the fly ( Christoffel symbols  , Ricci  and Riemann  tensors orthonormal and null tetrads , etc.)
 • Almost all of the mathematical operations one can perform on these solutions are implemented as commands in the Physics  and DifferentialGeometry  packages.
 • All the mathematics within the Maple library are instantly ready to work with these solutions and derived mathematical objects.

Finally, in the Maple PDEtools package , we have all the mathematical tools to tackle the equivalence problem around these solutions.

 Examples
 References

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

## need to procedure for calculating......

hi....how i can extract Coefficients  (i.e. {f1[2],f2[2],f2[3],f3[2],.....f3[6]}) from every algebric equations and create matrix A ,in form AX=0, (X are f1[2],f2[2],f2[3],f3[2],.....f3[6] ) then the determinant of the matrix of coefficients (A) set to zero for obtaining unknown parameter omega.?

Note that  if m=3 then 6 equations is appeare and if m=4 then 9 equations is appeare.thus i need a procedure that works for every arbitary value of ''m''.

in attached file below m=4 thus we have 9 equations, i.e. 3 for eq1[k_] and 3 for eq2[k_] and so on...

also we should use boundary conditions for some amount of fi[j] (i=1,2,3 and j=2,3,...,7)

be extacting above Coefficients for example from first equation ,

''**:= (1/128)*f1[2]*omega^2-(1/4)*f2[2]-(1/2)*f2[3]+(1/4)*f2[4]+(1/4)*f3[2]-(1/2)*f3[3]+(1/4)*f3[4]+140*f1[2]-80*f1[3]+20*f1[4]'''

must compute

coeff(**, f1[2]); coeff(**, f2[2]) and so on...

fdm-maple.mw

############################Define some parameters

############################Define some equation

######################################  APPLY BOUNDARY CONDITIONS

 (1)

## how convert 3 couple equations to 1 equation...

hi...how i can convert 3 couple equations to 1 equation with Placement each other?

thanks...

 (1)

 (2)

 (3)

## Solving a nonlinear system of 4 equations using sh...

Dear All,

I have a problem solving the attached nonlinear system of equations using shooting method.
I will be grateful if you could help me finding the solutions out.

restart; Shootlib := "C:/Shoot9"; libname := Shootlib, libname; with(Shoot);
with(plots);
N1 := 1.0; N2 := 2.0; N3 := .5; Bt := 6; Re_m := N1*Bt; gamma1 := 1;
FNS := {f(eta), fp(eta), fpp(eta), g(eta), gp(eta), m(eta), mp(eta), n(eta), np(eta), fppp(eta)};
ODE := {diff(f(eta), eta) = fp(eta), diff(fp(eta), eta) = fpp(eta), diff(fpp(eta), eta) = fppp(eta), diff(g(eta), eta) = gp(eta), diff(gp(eta), eta) = N1*(2.*g(eta)+(eta-2.*f(eta)).gp(eta)+2.*g(eta)*fp(eta)+2.*N2.N3.(m(eta).np(eta)-n(eta).mp(eta))), diff(m(eta), eta) = mp(eta), diff(mp(eta), eta) = Re_m.(m(eta)+(eta-2.*f(eta)).mp(eta)+2.*m(eta)*fp(eta)), diff(n(eta), eta) = np(eta), diff(np(eta), eta) = Re_m.(2.*n(eta)+(eta-2.*f(eta)).np(eta)+2.*N2/N3.m(eta).gp(eta)), diff(fppp(eta), eta) = N1*(3.*fpp(eta)+(eta-2.*f(eta)).fppp(eta)-2.*N2.N2.m(eta).(diff(mp(eta), eta)))};
blt := 1.0; IC := {f(0) = 0, fp(0) = 0, fpp(0) = alpha1, g(0) = 1, gp(0) = beta1, m(0) = 0, mp(0) = beta2, n(0) = 0, np(0) = beta3, fppp(0) = alpha2};
BC := {f(blt) = .5, fp(blt) = 0, g(blt) = 0, m(blt) = 1, n(blt) = 1};
infolevel[shoot] := 1;
S := shoot(ODE, IC, BC, FNS, [alpha1 = 1.425, alpha2 = .425, beta1 = -1.31, beta2 = 1.00, beta3 = 1.29]);
Error, (in isolate) cannot isolate for a function when it appears with different arguments
p := odeplot(S, [eta, fp(eta)], 0 .. 15);
Error, (in plots/odeplot) input is not a valid dsolve/numeric solution
display(p);
Error, (in plots:-display) expecting plot structure but received: p
p2 := odeplot(S, [eta, theta(eta)], 0 .. 10);
Error, (in plots/odeplot) input is not a valid dsolve/numeric solution
display(p2);
Error, (in plots:-display) expecting plot structure but received: p2

by: Maple 15

General description of the method of solving underdetermined systems of equations. As a particular application of the idea proposed a universal method  kinematic analysis for all kinds of linkage (lever) mechanisms. With the description and examples.

## solving system of non linear equations...

hello every one.please help me with solving this equations.i can not solve this and i need it.thanks

eq1 := (cos(beta2)-1)*w11-sin(beta2)*w12+(cos(alpha2)-1)*z11-sin(alpha2)*z12-cos(delta2) = 0; eq2 := (cos(beta2)-1)*w12+sin(beta2)*w11+(cos(alpha2)-1)*z12+sin(alpha2)*z11-2*sin(delta2) = 0; eq3 := (cos(beta3)-1)*w11-sin(beta3)*w12+(cos(alpha3)-1)*z11-sin(alpha3)*z12-3*cos(delta3) = 0; eq4 := (cos(beta3)-1)*w12+sin(beta3)*w11+(cos(alpha3)-1)*z12+sin(alpha3)*z11-4*sin(delta3) = 0; eq5 := (cos(beta4)-1)*w11-sin(beta4)*w12+(cos(alpha4)-1)*z11-sin(alpha4)*z12-5*cos(delta4) = 0; eq6 := (cos(beta4)-1)*w12+sin(beta4)*w11+(cos(alpha4)-1)*z12+sin(alpha4)*z11-6*sin(delta4) = 0; eq7 := (cos(beta5)-1)*w11-sin(beta5)*w12+(cos(alpha5)-1)*z11-sin(alpha5)*z12-7*cos(delta5) = 0; eq8 := (cos(beta5)-1)*w12+sin(beta5)*w11+(cos(alpha5)-1)*z12+sin(alpha5)*z11-8*sin(delta5) = 0; alpha2 := -20; alpha3 := -45; alpha4 := -75; alpha5 := -90; delta2 := 15.5; delta3 := -15.9829; delta4 := -13.6018; delta5 := -16.7388; P21 = .5217; P31 = 1.3421; P41 = 2.3116; P51 = 3.1780;

## simultaneous equation...

Hi,

I have attached a Maple file. My problem is that the solve for the simultaneous equation does not give me understandable results. I even simplified my equations by saying some parameters are zero although my final goal is to find an expression for a and varphi. Any idea how to solve this analytically? I know how to do it numerically. I need an analytical expression.

Thanks,

Baharm31

## Roots of the characteristic polynomial ...

`how i can calculate roots of the characteristic polynomial equations {dsys and dsys2}and dsolve them with arbitrary initial condition for differennt amont of m and n?thanksKr.mw  `

 (1)

 (2)

 (3)

` `

` `

## Difficulty in Solving a Problem...

I have been working on a general solution to motion analysis and seem to be going backwards.  I have an numerical solution in Octave I use for comparison.  I have reduced the problem to a small example that exhibits the problem.

I posted a question similar to this, but, without a set of known values.

I am doing something wrong, but, what?

Tom Dean

## bearing.mpl, solve the target motion problem with bearings only.
##
## Consider a sensor platform moving through points (x,y) at times
## t[1..4] with the target bearings, Brg[1..4] taken at times t[1..4]
## with the target proceeding along a constant course and speed.
##
## time t, bearing line slope m, sensor position (x,y) are known
## values.
##
## Since this is a generated problem the target position at time t is
## provided to compare with the results.
##
#########################################################################
##
restart;
##
genKnownValues := proc()
description "set the known values",
"t - relative time",
"x - sensor x location at time t[i]",
"y - sensor y location at time t[i]",
"m - slope of the bearing lines at time t[i]",
"tgtPosit - target position at time t[i]";
global t, m, x, y, tgtPosit;
local dt, Cse, Spd, Brg, A, B, C, R, X;
local tgtX, tgtY, tgtRange, tgtCse, tgtSpd;
## relative and delta time
t := [0, 1+1/2, 3, 3+1/2];
dt := [0, seq(t[idx]-t[idx-1],idx=2..4)];
## sensor motion
Cse := [90, 90, 90, 50] *~ Pi/180; ## true heading
Spd := [15, 15, 15, 22];  ## knots
## bearings to the target at time t
Brg := [10, 358, 340, 330] *~ (Pi/180);
## slope of the bearing lines
m:=map(tan,Brg);
## calculate the sensor position vs time
x := ListTools[PartialSums](dt *~ Spd *~ map(cos, Cse));
y := ListTools[PartialSums](dt *~ Spd *~ map(sin, Cse));
## target values  start the target at a known (x,y) position at a
## constant course and speed
tgtRange := 95+25/32; ## miles at t1, match octave value...
tgtCse := 170 * Pi/180; ## course
tgtSpd := 10; ## knots
tgtX := tgtRange*cos(Brg[1]);
tgtX := tgtX +~ ListTools[PartialSums](dt *~ tgtSpd *~ cos(tgtCse));
tgtY := tgtRange*sin(Brg[1]);
tgtY := tgtY +~ ListTools[PartialSums](dt *~ tgtSpd *~ sin(tgtCse));
## return target position vs time as a matrix
tgtPosit:=Matrix(4,2,[seq([tgtX[idx],tgtY[idx]],idx=1..4)]);
end proc:
##
#########################################################################
## t[], m[], x[], and y[] are known values
##
## equation of the bearing lines
eq1 := tgtY[1] - y[1]    = m[1]*(tgtX[1]-x[1]):
eq2 := tgtY[2] - y[2]    = m[2]*(tgtX[2]-x[2]):
eq3 := tgtY[3] - y[3]    = m[3]*(tgtX[3]-x[3]):
eq4 := tgtY[4] - y[4]    = m[4]*(tgtX[4]-x[4]):
## target X motion along the target line
eq5 := tgtX[2] - tgtX[1] = tgtVx*(t[2]-t[1]):
eq6 := tgtX[3] - tgtX[2] = tgtVx*(t[3]-t[2]):
eq7 := tgtX[4] - tgtX[3] = tgtVx*(t[4]-t[3]):
## target Y motion along the target line
eq8 := tgtY[2] - tgtY[1] = tgtVy*(t[2]-t[1]):
eq9 := tgtY[3] - tgtY[2] = tgtVy*(t[3]-t[2]):
eq10:= tgtY[4] - tgtY[3] = tgtVy*(t[4]-t[3]):
##
#########################################################################
##
## solve the equations
eqs  := {eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9,eq10}:

Sol:= solve(eqs, {tgtVx, tgtVy, seq([tgtX[k], tgtY[k]][], k= 1..4)}):
##

genKnownValues():
## these values are very close to Octave
evalf(t);evalf(m);evalf(x);evalf(y);evalf(tgtPosit);
## The value of tgtX[] and tgtY[] should equal the respective tgtPosit values
seq(evalf(eval([tgtX[idx],tgtY[idx]], Sol)),idx=1..4);

 1 2 3 4 5 6 7 Last Page 1 of 58
﻿