How do people test BIBD design in maple?

what are the different when using different field when test design ?

how to design a game that is suitable for the specified field of the design?

I have assumed lambda <-1 yet signum(lambda) returns lambda~.  It should be -1. I need to show that the expression lambda is used in is negative. 

Download signum_problem.mw

I am trying to add a specific property (radiation length) to the Elements in ScientificConstants.

To do this  I first 

ScientificConstants:-AddProperty(radiationlength); # this seems to work

I then load a table of element name - radiationLength pairs (it is really a bigger table but only these are of interest):

data:=ImportMatrix(cat(libname[1],"/RadiationLengthTable.txt"),source=delimited,\
                       datatype=anything,delimiter="\t"); # works too

(I made the file using downloaded data).

I then want to fill the radiationlength properties for all elements for which I've got the value:

for i from 1 to LinearAlgebra:-RowDimension(data) do
  elemt:=data[i,1]; # the element name as a string
  ScientificConstants:-ModifyElement(cat(elemt),radiationlength=[value=data[i,4],units=g/cm^2]);
end do;

which fails with the error

Error, (in ScientificConstants:-ModifyElement) element ``H `` is not known

H is the first element in the table of radiation lengths so it is getting the right one.

At  issue is that somehow the element name (called descriptor in the Help files) does not seem to get properly evaluated. In fact I tried various ways "by hand", none except just the unevaluated element name (i.e H,He etc. worked).

There must be a way to set these in a loop. Does someone know how?

TIA,

Mac Dude

help me

into

Why I could not unaasign U_hat in my following program?

Download moshkel.mw

Here is a simple little procedure that does not work

et :=proc(x) 

description   "this is a variant of evalf that gets rid of almost 0 nos.(rounding error) and shortens the display to 2 `digits"`;  

if  abs(x) < 10^((-14))  then 0 else  evalf(x,2) end if     end proc;

Error, unable to delimit strings/identifiers

and here is the same procedure that works because I added the # before the description.

et :=proc(x)

# description   "this is a variant of evalf that gets rid of almost 0 nos.(rounding error) and shortens the display to 2 `digits"`;  

if  abs(x) < 10^((-14))  then 0 else  evalf(x,2) end if     end proc;

In the help pages they do not say that the description must be hidden by the # sign (nor do they mention that quotes should be used although they are used in the examples.

What am I missing?

How can I activate my maple 13?

I have no activatIon left but installed maple

once

My purchase Code is BNN7VMT4BPFN7ESG

Thank. You 

Hi,

Is there any one can tell me how to plot a curve(e.g f(x,1)=2x^2-x-1) by the command spacecurve?

Thanks in advance!

I would like to compute the function 

                       f(x) := exp(x) * HeunC(2*x, 1/2, -1/2, -x^2, 1/8, 99/100)

for very large real values of x. However, in the following plot it can be seen that there are already strong numerical instabilities for x>33. What can I do to reach more stability for large x? Or, does somebody know how to get a closed form asymptotic expression for large x?

I am working on a simple set of notes for students. Much of the on-line stuff is either out-of-date or too complex.

Here is the problem I wish to show

A corpse was found in a room with an ambient temperature of 64°F. The body temperature on finding the corpse was 75°F; one hour later it was 73°F. Use Newtons Law of cooling to find the time of death.

Here is my solution. I realise the assignment to variable cool is not used. Please tell me any suggestions for improvements.

many thanks

Hi,

When you do this a := plot(1/x, x=0..1);  the figure contains a "smart" graph which extends roughly from 0 to 30 in the vertical direction.
(no discont=true nor smartview=false used here, numpoints set to its default falue)

If you use plottools:-getdata(a); you find that the vertical range is about  0..1800.
It seems to mean that smartview=true (the default setting) overrides the range determined from the values of discont and numpoints?
 

My question is: Is it possible to retrieve the vertical range that plot uses when it displays the graph?
 

Dear all,

I need your help to compute the multiplicative  group generated by the following two matrices.

A :=  Matrix(2, 2, [[2, 0], [0, 1]]) ;

B := Matrix(2, 2, [[1, 1], [0, 1]]) ;

 such that A.B different to B.A

many thanks

Hello Dr/Pof/Colleague

Please help me on the ODE BVP problem

my target i want RUN double loop

Download 3DAKc2alpha.mw

Here is a simple collection of Python scripts that allows Maple code to be embedded directly within a LaTeX document. The scripts will process the LaTeX document, calling Maple as required, making the Maple output directly accessible in the LaTeX document. The source, including extensive examples, can be found at https://github.com/leo-brewin/hybrid-latex. This  library also supports inclusion of active Mathematica, Matlab, Python and Cadabra code within LaTeX documents.

Here is a simple example (based on maple/example-01 in the GitHub site)

\documentclass[12pt]{mpllatex}

\begin{document}

\section*{Calculus}

\begin{maple}
   ans := diff(x*sin(x),x):                          # mpl (ans.501,ans)
   ans := eval(diff(x*sin(x),x),x=Pi/4):             # mpl (ans.502,ans)
   ans := int(2*sin(x)^2, x=a..b):                   # mpl (ans.503,ans)
   ans := int(2*exp(-x^2),x=0..infinity):            # mpl (ans.504,ans)
   ans := ''int(2*exp(-x^2),x=0..infinity)'':        # mpl (lhs.504,ans)
   ans := int(int(x^2 + y^2,  y=0..x),x=0..1):       # mpl (ans.505,ans)
   ans := ''int(int(x^2 + y^2,  y=0..x),x=0..1)'':   # mpl (lhs.505,ans)
\end{maple}

\begin{align*}
   &\mpl*{ans.501}\\
   &\mpl*{ans.502}\\
   &\mpl*{ans.503}\\
   \mpl{lhs.504}&=\Mpl{ans.504}\\
   \mpl{lhs.505}&=\Mpl{ans.505}
\end{align*}

\end{document}

and here is the corresponding output

example-01.pdf