MaplePrimes Questions

I am solving a complicated ODE and I would like to know if there is a way for Maple to output the ODEs without doing any numerical substitutions for known parameters. Say one of my parameters, call it P, is initialized (there are many more but Ill just simplify and consider one here) to a value of 10. In order to chekc that I have coded the ODEs correctly it would help me if Maple does not substitue with the numerical value 10 for P when displayng teh ODEs, but rather keeps P, as a parameter. Is there a way to achieve this?

hi.please help me for solve algebric equations below?

solve or fsolve dos not any answer?

i think use Newton  Raphson, because solve or fsolve not work.

thankscomparision.mw

restart; mu1 := .1; mu2 := .1; sigma1 := -40; F := 25; upsilon1 := (1/2)*sigma2; upsilon2 := (3/2)*sigma2-sigma1; sigma2 := 100

100

(1)

gamma11 := 2.686901; -1; gamma12 := 7.175339; -1; gamma21 := 2.436735; -1; gamma22 := 12.94855; -1; gamma22 := 12.94855; -1; delta1 := .928207; -1; delta2 := .105073; -1; s11 := .629894; -1; s12 := .217601; -1; s21 := 0.73897e-1; -1; s22 := .805815

.805815

(2)

 

 

 

 
Q1 := -mu1*p1-upsilon1*q1+gamma11*q1*(p1^2+q1^2)+gamma12*q1*(p2^2+q2^2)-delta1*(2*p1*q1*p2-q2*(p1^2-q1^2))-s11*F*q1+s12*F*q2 = 0

-.1*p1-65.747350*q1+2.686901*q1*(p1^2+q1^2)+7.175339*q1*(p2^2+q2^2)-1.856414*p1*q1*p2+.928207*q2*(p1^2-q1^2)+5.440025*q2 = 0

(3)

Q2 := -mu1*q1+upsilon1*p1-gamma11*p1*(p1^2+q1^2)-gamma12*p1*(p2^2+q2^2)-delta1*(2*p1*q1*q2+p2*(p1^2-q1^2))-s11*F*p1-s12*F*p2 = 0

-.1*q1+34.252650*p1-2.686901*p1*(p1^2+q1^2)-7.175339*p1*(p2^2+q2^2)-1.856414*p1*q1*q2-.928207*p2*(p1^2-q1^2)-5.440025*p2 = 0

(4)

Q3 := -mu2*p2-upsilon2*q2+gamma21*q2*(p1^2+q1^2)+gamma22*q2*(p2^2+q2^2)+delta2*(3*p1^2*q1-q1^3)+s21*F*q1 = 0

-.1*p2-190*q2+2.436735*q2*(p1^2+q1^2)+12.94855*q2*(p2^2+q2^2)+.315219*p1^2*q1-.105073*q1^3+1.847425*q1 = 0

(5)

Q4 := -mu2*q2+upsilon2*p2-gamma21*p2*(p1^2+q1^2)-gamma22*p2*(p2^2+q2^2)+delta2*(-p1^3+3*p1*q1^2)-s21*F*p1 = 0

-.1*q2+190*p2-2.436735*p2*(p1^2+q1^2)-12.94855*p2*(p2^2+q2^2)-.105073*p1^3+.315219*q1^2*p1-1.847425*p1 = 0

(6)

 

-.1*q2+((3/2)*sigma2-40)*p2-2.436735*p2*(p1^2+q1^2)-12.94855*p2*(p2^2+q2^2)-.105073*p1^3+.315219*q1^2*p1-1.847425*p1

(7)

NULL

fsolve({Q1, Q2, Q3, Q4}, {p1, p2, q1, q2})

fsolve({-.1*p1-65.747350*q1+2.686901*q1*(p1^2+q1^2)+7.175339*q1*(p2^2+q2^2)-1.856414*p1*q1*p2+.928207*q2*(p1^2-q1^2)+5.440025*q2 = 0, -.1*p2-190*q2+2.436735*q2*(p1^2+q1^2)+12.94855*q2*(p2^2+q2^2)+.315219*p1^2*q1-.105073*q1^3+1.847425*q1 = 0, -.1*q1+34.252650*p1-2.686901*p1*(p1^2+q1^2)-7.175339*p1*(p2^2+q2^2)-1.856414*p1*q1*q2-.928207*p2*(p1^2-q1^2)-5.440025*p2 = 0, -.1*q2+190*p2-2.436735*p2*(p1^2+q1^2)-12.94855*p2*(p2^2+q2^2)-.105073*p1^3+.315219*q1^2*p1-1.847425*p1 = 0}, {p1, p2, q1, q2})

(8)

solve(Q1, Q2, Q3, Q4)

Error, invalid input: too many and/or wrong type of arguments passed to solve; first unused argument is -.1*q1+34.252650*p1-2.686901*p1*(p1^2+q1^2)-7.175339*p1*(p2^2+q2^2)-1.856414*p1*q1*q2-.928207*p2*(p1^2-q1^2)-5.440025*p2 = 0

 

Sol := [fsolve(ZZ, omega)]



Download comparision.mw

if three particles set into 3 energy states ( each values of energy states 0 , e, 2e)

 

  

 

 

how to calculate P in each configuration by a particle that follows boson and f-b

 

 

 

thank you

I'm having trouble plotting a couple things. I have

eq := diff(y(x), x$3)+3*diff(y(x),x$2)+12*y(x);
soln := dsolve(eq, y(x));
soln := evalf(soln);
PartSoln1 := dsolve({eq, y(0) = a,y'(0) = 0,y''(0) = 0}, y(x));
curves := {seq(PartSoln1, a = -3 .. 3)};
Then when I try plot(curves, x = -1..5, y = -5..5); I get Warning, expecting only range variable x in expression PartSoln1 to be plotted but found name PartSoln1.

Also,
charEq := r^3+r+1 = 0;
soln := solve(charEq, r);
soln := [evalf(soln, 5)];
soln := map(Re, soln);

I tried a few things, but can't figure out how to plot charEq, including the real roots.

Thanks for any help,

Heather

Here is a code

Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/EQ.mw .

Download EQ.mw

A fragment of code

for b in extra_bcs do try print(b = 10^(-2)); res[b] := dsolve(dsys4 union {b = 10^(-2)}, numeric, initmesh = 2024, output = listprocedure, approxsoln = [omega2 = 0.1e-2, s(x) = cosh(upsilon*x)-cos(upsilon*x)-(cosh(upsilon)+cos(upsilon))*(sinh(upsilon*x)-sin(upsilon*x))/(sinh(upsilon)+sin(upsilon)), g(x) = sin(((2*n+1)*(1/2))*Pi)], abserr = 0.1e-1) catch: print(lasterror) end try end do; indx := indices(res, nolist); nops([indx]); res[indx[i]]; seq(subs(res[indx[i]](1), omega2(1)), i = 1 .. nops([indx]))

Hello everybody,

I have an interesting math problem, here it is:

 

Add the following operators: +, -, * in the middle of the following string 123456789 to get the final result: 2016

The operators can be added randomly at any place, it isn't necessarily placed in between single numbers.

 

To convince your result, one should provide the code or explain explicitly your method.

Have fun!!!!

 

 

hi.

how i can select or chose proper polynomials or another functions that attached boundary conditions at points zero and one , weakly or strongly satisfy??polynomial.mw

s(0) = 0, ((D@@1)(s))(0) = 0, g(0) = 0, ((D@@2)(s))(1) = 0, ((D@@3)(s))(1) = 0, ((D@@1)(g))(1)+(1/2)*((D@@1)(s))(1)^2 = 0

s(0) = 0, (D(s))(0) = 0, g(0) = 0, ((D@@2)(s))(1) = 0, ((D@@3)(s))(1) = 0, (D(g))(1)+(1/2)*(D(s))(1)^2 = 0

(1)

``

 

Download polynomial.mw

thanks...

hi.i calculate this integral numerically but after use this result in my differential equation which is involve f3(x), i enconter with error

i found that this problem is due to answer of this integtal.may i use another way to calculte it?as shown power in result are very large!!!

``

restart; Digits := 15; g3 := theta^2*(theta-1)^2; beta := 100; chi := 5; kappa := 5; a := 0

0

(1)

with(Student[Calculus1]); -1; a31 := evalf(ApproximateInt(-beta^2*g3/((1-g3*f3(x))*ln(2*kappa*(1-g3*f3(x)))^2), theta = a .. 1, method = simpson)); -1; a32 := evalf(ApproximateInt(-chi*g3/(1-g3*f3(x))^4, theta = a .. 1, method = simpson))

-1.33333333333333*(-0.191597582462488e152*f3(x)-0.863011826221205e146*f3(x)^9-0.706378636532818e148*f3(x)^7+0.846944805095211e147*f3(x)^8+0.130151241709538e151*f3(x)^4-0.283379399043279e150*f3(x)^5+0.493554199542613e149*f3(x)^6+0.117104095781373e152*f3(x)^2-0.459581479474515e151*f3(x)^3+0.151119505345015e152+0.755561753647024e145*f3(x)^10+0.175123383079432e97*f3(x)^36-0.244416009548583e100*f3(x)^35+0.163803382818292e103*f3(x)^34-0.704973843312161e105*f3(x)^33+0.219681545650952e108*f3(x)^32-0.529211072065593e110*f3(x)^31+0.102639576866076e113*f3(x)^30-0.164658103849427e115*f3(x)^29+0.222629381398692e117*f3(x)^28-0.257124576224432e119*f3(x)^27+0.256176429040311e121*f3(x)^26-0.221804182616290e123*f3(x)^25+0.167834727658308e125*f3(x)^24-0.111474607475097e127*f3(x)^23+0.652137994527267e128*f3(x)^22-0.336921770775584e130*f3(x)^21+0.154034904897255e132*f3(x)^20-0.624063631993155e133*f3(x)^19+0.224251041670249e135*f3(x)^18-0.714943093769717e136*f3(x)^17+0.202175772822691e138*f3(x)^16-0.506683948932158e139*f3(x)^15+0.112368506172683e141*f3(x)^14-0.220033360070981e142*f3(x)^13+0.379271084508473e143*f3(x)^12-0.573179869691401e144*f3(x)^11)/((-160000.+9801.*f3(x))^4*(-16.+f3(x))^4*(-256.+9.*f3(x))^4*(-10000.+441.*f3(x))^4*(-160000.+8281.*f3(x))^4*(-625.+36.*f3(x))^4*(-160000.+361.*f3(x))^4*(-10000.+81.*f3(x))^4*(-160000.+2601.*f3(x))^4*(-625.+16.*f3(x))^4)

(2)

``

 

Download integral.mwintegral.mw

thanks..

I have written a code which generates a Gaussian like curve as a set of points and have written some basic commands to find several important quantities from the plot. The three of importance here are the peak position, peak height and Full Width at Half Maximum.

See the attached minimal working example maple worksheet.

It can find the peak position and the peak height no problem, but another quantity I want it to find is the Full Width at Half Maximum (FWHM). I use the aforementioned calculated values to find the left and right components of the FWHM but the fsolve command just hangs.

In the worksheet this is not evident but in the much longer full code which I run through the terminal and prints the values to a file, it begins the fsolve command and reaches a particular value and it just hangs. No amount of waiting time makes it continue the calculation.

However if I use an approximate integration scheme (such as Riemann sums) then the fsolve part executed perfectly. fslove seems to get stuck during the numerical integration step.

There is probably an easy answer to this problem which I am missing. Any help is appreciated.

fsolve_numerical_integration.mw

 

- Yeti

is there any command to plot ternarry plots that show the relationsh between three sets of data? For example check the link

https://en.wikipedia.org/wiki/Ternary_plot

Hi,

I am solving a system of differential equations with 2 variables. I can find numerical solution using dsolve with numeric option. I can also plot curves. But I also would like to extract data at time t0 and store into a file. I don't know how to extract numerical values from dsolve and store into a text file at various time, t0, t1, t2 etc. I would like to store data in tabular format with first, sceond and third columns are time, x and y respectively.

 

Just as an example, we can use

x'(t) = y(t), 

y'(t) = 6x(t) - y(t)

with x(0)=1 and y(0) = 2

 

Thanks.

I need to calculate Weyl scalar for a metric using null tetrad using debever package in maple 13. However, I am stuck at the defination of h (representing the covariant complex null tetrad). Is it the product of covariant null tetrad? I have worked it out by using covariant, contravariant and both covariant and contravariant null tetrad (like l_a*l_a, l^a*l^a, l_a*l^a), however, i am not getting the right result not even for the example given in the maple help for plane wave. please help me out that how should i define this h.

thanks,
suresh

I have a polynomial with parameter. I hope to solve the equation that every coefficient of the polynomial is zero.

For example

f(t):=(a+b+1)t+(a-b-1)

then I use the function solve(f(t),{a,b}),

I hope to get the solution to a+b+1=0 and a-b-1=0. However the output of Maple is {a = -(b*t-b+t-1)/(t+1), b = b}. How could regard t as a variable?

 

Thank you very much.

Dear All

Please see following:

 

Suppose we have

M := proc (i) options operator, arrow; i/(i+1) end proc

proc (i) options operator, arrow; i/(i+1) end proc

(1)

and

F := proc (i, j) options operator, arrow; M(i)/(M(i)+j) end proc

proc (i, j) options operator, arrow; M(i)/(M(i)+j) end proc

(2)

 For i from 1 to 5 and j from 1 to 5, how one can display all values of F(i,j) using "loop" ???

 

Download loop.mw

Regards

Dear Maple users

I wanted to use the zip-command to muliply two lists with units. It works without units, but not with units:

 

 

If I add units to each entrance it works too, but not very convenient. How is it done with units in the best way?

 

Regards,

Erik

First 1145 1146 1147 1148 1149 1150 1151 Last Page 1147 of 2434