evalf(10.2^20, 50);

evalf((10+1/5)^20, 50);

Where are 50 digits of the first result?

I have a matrix (1002,2) where first column isthe time and the second is numerical solution of an angle.My angle should be oscillating periodically but I am getting a noise at large times. I need to check if this noise contains certain frequencies in it so I need to plot a power spectrum of my data. How can I do that 

> restart;
> with(LinearAlgebra);
> q := a*mu^4+b*mu^3+d*mu^2+e*mu+f = 0;
> sol := solve(a*mu^4+b*mu^3+d*mu^2+e*mu+f = 0, mu);
> S[1] := allvalues(sol);
> PARAM := [a = -54/c^2-1269*A[1]/(8*c^2), b = 108/c^2+5013*A[1]/(8*c^2), d = 27-693/(2*c^2)+117*A[1]-7113*A[1]/(4*c^2), e = -27+585/(2*c^2)-111*A[1]+20439*A[1]/(16*c^2), f = 1-3*A[1]-18/c^2-8*A[1]/c^2];
> S2 := eval(S[1], PARAM);
Error, invalid input: eval received PARAM, which is not valid for its 2nd argument, eqns
> series(%, c = infinity, 3, A[1] = 0, 2]);
>
Please i mean by last statemene the solutions to be written i series as:

m1=.........+O(1/c^3)+O(A1^2)

m2=.......+O(1/c^3)+O(A1^2)

etf.

A1=0,2 means O(A1^2),therefore,i am not evaluating at .2 .

1. how to do optimization with partial differential equation as constraints in maple

2. how to do optimization with  partial differential equation as objective function in order to make output obey this model

how gr operators work?

I tried to run example given there but it is not working,

where can I get more examples to understand working of Gr operators work?

specially for raychaudhuri equations.

Hi MaplePrime-ers,

I'm using Maple17 in Matlab 2012b to evaluate a symbolic function over a grid of values.  The number of values is generally ~500k, and therefore I have a loop that dumps the solutions into MATLAB where the values are parsed and stored more efficiently.

I put this proceess in an optimization routine, and I keep getting crashes.  The crashes are NOT repeatable.  They generally happen after 4-10 times the routine has been called.  

This is one of the stack dumps.  Anybody have any ideas?  I talked to MATLAB support, but they weren't very useful (they started pointing fingers).

------------------------------------------------------------------------
Segmentation violation detected at Wed Jun 4 17:38:11 2014
------------------------------------------------------------------------

Configuration:
Crash Decoding : Disabled
Current Visual : 0x24 (class 4, depth 24)
Default Encoding: UTF-8
GNU C Library : 2.15 stable
MATLAB Root : /opt/Matlab/R2012b
MATLAB Version : 8.0.0.783 (R2012b)
Operating System: Linux 3.2.0-37-generic #58-Ubuntu SMP Thu Jan 24 15:28:10 UTC 2013 x86_64
Processor ID : x86 Family 6 Model 23 Stepping 6, GenuineIntel
Virtual Machine : Java 1.6.0_17-b04 with Sun Microsystems Inc. Java HotSpot(TM) 64-Bit Server VM mixed mode
Window System : The XFree86 Project, Inc (40300000), display wildnode0:15.0

Fault Count: 1

Abnormal termination:
Segmentation violation

Register State (from fault):
RAX = 00007f3bcc0a80c0 RBX = 00007f3b97fe6000
RCX = 0000000000000000 RDX = 00007f3bcc0a84a0
RSP = 00007f3b8cb29d60 RBP = 0000000000000003
RSI = 5851f42d4c957f2d RDI = 00007f3bde026360

R8 = 65a1566174cc9e28 R9 = 0000000000000002
R10 = 00007f3b8cb29d90 R11 = 0000000000002e88
R12 = 0000000000000001 R13 = 00007f3b8cb29df0
R14 = 0000000000000000 R15 = 0000000000000003

RIP = 00007f3b977d0604 EFL = 0000000000010283

CS = 0033 FS = 0000 GS = 0000

Stack Trace (from fault):
[ 0] 0x00007f3c962b01de /opt/Matlab/R2012b/bin/glnxa64/libmwfl.so+00516574 _ZN2fl4diag15stacktrace_base7captureERKNS0_14thread_contextEm+000158
[ 1] 0x00007f3c962b14b2 /opt/Matlab/R2012b/bin/glnxa64/libmwfl.so+00521394
[ 2] 0x00007f3c962b2ffe /opt/Matlab/R2012b/bin/glnxa64/libmwfl.so+00528382 _ZN2fl4diag13terminate_logEPKcRKNS0_14thread_contextE+000174
[ 3] 0x00007f3c9559f093 /opt/Matlab/R2012b/bin/glnxa64/libmwmcr.so+00557203 _ZN2fl4diag13terminate_logEPKcPK8ucontext+000067
[ 4] 0x00007f3c9559bb9d /opt/Matlab/R2012b/bin/glnxa64/libmwmcr.so+00543645
[ 5] 0x00007f3c9559d835 /opt/Matlab/R2012b/bin/glnxa64/libmwmcr.so+00550965
[ 6] 0x00007f3c9559da55 /opt/Matlab/R2012b/bin/glnxa64/libmwmcr.so+00551509
[ 7] 0x00007f3c9559e0fe /opt/Matlab/R2012b/bin/glnxa64/libmwmcr.so+00553214
[ 8] 0x00007f3c9559e295 /opt/Matlab/R2012b/bin/glnxa64/libmwmcr.so+00553621
[ 9] 0x00007f3c93a86cb0 /lib/x86_64-linux-gnu/libpthread.so.0+00064688
[ 10] 0x00007f3b977d0604 /opt/maple17/bin.X86_64_LINUX/libmaple.so+01824260
[ 11] 0x00007f3b977cfc48 /opt/maple17/bin.X86_64_LINUX/libmaple.so+01821768
[ 12] 0x00007f3b977d2f88 /opt/maple17/bin.X86_64_LINUX/libmaple.so+01834888
[ 13] 0x00007f3b9763d63f /opt/maple17/bin.X86_64_LINUX/libmaple.so+00173631
[ 14] 0x00007f3c93a7ee9a /lib/x86_64-linux-gnu/libpthread.so.0+00032410
[ 15] 0x00007f3c937ab3fd /lib/x86_64-linux-gnu/libc.so.6+01000445 clone+000109

If this problem is reproducible, please submit a Service Request via:
http://www.mathworks.com/support/contact_us/

A technical support engineer might contact you with further information.

Thank you for your help.

I am trying to design a question to test for divisibility

$n3=range(0,9,1); and students have to choose the next digit to make the number divisible by 4

if $n3 is even the students must input 0,4,8 to be true

if $n3 is odd the students must input 2 or 6 to be correct.

I am not sure how to us the if(a,b,c) in this case or if I should be using indexof....or both?  Any help would be appreciated. 

> with(Statistics);
> X := Vector([1, 2, 3, 4, 5, 6, 7], datatype = float);
> Y := Vector([155625, 172472, 179589, 186579, 205421, 214989, 237937], datatype = float);
> ExponentialFit(X, Y, x, output = residualsumofsquares);
   
Result by Maple: 0.00216641200893470318
But a real result is 8.159611742*10^7

Please, can anybody help me? Thank you.

> restart;
> with(LinearAlgebra);
> q := a*mu^4+b*mu^3+d*mu^2+e*mu+f = 0;
> sol := solve(a*mu^4+b*mu^3+d*mu^2+e*mu+f = 0, mu);
> S[1] := allvalues(sol);
> PARAM := [a = -54/c^2-1269*A[1]/(8*c^2), b = 108/c^2+5013*A[1]/(8*c^2), d = 27-693/(2*c^2)+117*A[1]-7113*A[1]/(4*c^2), e = -27+585/(2*c^2)-111*A[1]+20439*A[1]/(16*c^2), f = 1-3*A[1]-18/c^2-8*A[1]/c^2];
> S2 := eval(S[1], PARAM);
> series(%, c = infinity, 3, A[1] = 0, 2);
> simplify(%)

I was given that solitary initial state to see how it will deform as time goes on. am struggling to get my code so that I can get video frames. please help on how I can generate my code

u_o(x)= a₀x²(1-x)² for x (less than or equal to) x (less than or equal) 1

u₀(x)    = 0 for x > 1

the video clips will be representing the function u(*,t) :x to u(x,t)

for a sequence of choices of t such as t=0; t=0,5...t=3

 “Error, too many levels of recursion” ,but there’s just a single expression.

How can this happen?

I am using Maple 15 to numerically solve a system of differential algebraic euqations (DAE) with given initial conditions, and I've tried rfk45_dae and rosenbrock_dae solver, but both solver responded in error like this

Error, (in dsolve/numeric) cannot numerically solve complex DAE initial value problems, the system must be converted to a real system

I don't understand what is a real system, and how could i convert it to a real system.

I'm running a long loop with a lot of variables. I have noticed that Maple does not seem to clear its memory after each loop iteration, which results in very large memory used after all the iterations. I can't use restart because I need a few variables to be remembered throughout the loop for analysis afterwards. So the question is how to clear the memory for selective variables and procedures at the end of each iteration.

How to make a usual book view of a system of equations? I want that each equations is placed on new line with curly brace on the left.

I need to plot a function:



f(1.0) gives a correct result but plot doesn't. Could you help me with it?