Lecture 4 - Problem 2
a)
> a := -2100; b := 2100; N := 8; h := (b-a)/N; C := 4491; x := [seq(h, h = -2100 .. 2100, h)]; y := proc (x) options operator, arrow; sqrt(1+((1/2)*exp(x/C)-(1/2)*exp(-x/C))^2) end proc; evalf(y(x[1])); evalf(y(x[N+1])); L := (1/3)*h*(y(x[1])+4*(y(x[2])+y(x[4])+y(x[6])+y(x[N]))+2*(y(x[3])+y(x[5])+y(x[N-1]))+y(x[N+1])); evalf(L);
                             -2100
                              2100
                               8
                              525
                              4491
     [-2100, -1575, -1050, -525, 0, 525, 1050, 1575, 2100]
         /                           2\
         |    /1    /x\   1    /  x\\ |
x -> sqrt|1 + |- exp|-| - - exp|- -|| |
         \    \2    \C/   2    \  C// /
                          1.111332350
                          1.111332350
                                           (1/2)
         /                               2\    
         |    /1    /-700\   1    /700 \\ |    
     175 |1 + |- exp|----| - - exp|----|| |    
         \    \2    \1497/   2    \1497// /     

                                               (1/2)
              /                              2\    
              |    /1    /-175\   1    /175\\ |    
        + 700 |1 + |- exp|----| - - exp|---|| |    
              \    \2    \499 /   2    \499// /     

                                                (1/2)
              /                               2\    
              |    /1    /-175\   1    /175 \\ |    
        + 700 |1 + |- exp|----| - - exp|----|| |    
              \    \2    \1497/   2    \1497// /     

                                                (1/2)
              /                               2\    
              |    /1    /175 \   1    /-175\\ |    
        + 700 |1 + |- exp|----| - - exp|----|| |    
              \    \2    \1497/   2    \1497// /     

                                               (1/2)
              /                              2\    
              |    /1    /175\   1    /-175\\ |    
        + 700 |1 + |- exp|---| - - exp|----|| |    
              \    \2    \499/   2    \499 // /     

                                                (1/2)     
              /                               2\          
              |    /1    /-350\   1    /350 \\ |          
        + 350 |1 + |- exp|----| - - exp|----|| |      + 350
              \    \2    \1497/   2    \1497// /           

                                                (1/2)
              /                               2\    
              |    /1    /350 \   1    /-350\\ |    
        + 350 |1 + |- exp|----| - - exp|----|| |    
              \    \2    \1497/   2    \1497// /     

                                                (1/2)
              /                               2\    
              |    /1    /700 \   1    /-700\\ |    
        + 175 |1 + |- exp|----| - - exp|----|| |    
              \    \2    \1497/   2    \1497// /    
                          4354.742601
> L := (1/3)*h*(y(x[1])+4*add(y(i), i = y(x[2]) .. y(x[N]), 2)+2*add(y(i), i = y(x[3]) .. y(x[N-1]), 2)+y(x[N+1])); evalf(L);
Error, invalid input: add expects 2 arguments, but received 3
                          4354.742601

We've just decided to jump from Mathcad to Maple, and are doing our very first steps with Maple.

Our company is a small sized engineering company (15 employees), and what I am looking for is a bit help to make a working company setup.

This includes in the first run - how do I create and tell Maple where to find those items?

• Templates with company logo
• Programs and packages
• Are there any configuration variables that can / should be set?
• Are there any other common files that can be made accessible for all on a common network server?

Error

can not raise the datalimit above hard limit

my new notebook has 32GB memory

kernelopts(datalimit=28*unit(gibibtye)):

do i need to assign datalimit with kernelopts?

and

I remember i set a trick to protected this file

now I can read and run 

but return above error

does it mean that I do not need to set limit?

Minimise(eq1, {eq2>=0,eq3>=0,eq4>=0, a1<>0,a2<>0,a3<>0,a4<>0})

can not handle strict inequalities error

I created a user library and added a component that I had created from a shared subsystem and this worked fine.

I edited the component in the library to add some additional functionality and saved it.  I then had to reload the library in Maplesim but now I get this error everytime I try to run the simulation:

Cannot resolve `ControlLogix.PIDE` in model `Main`; there is no `ControlLogix` visible

Yet I dragged and dropped the component PIDE from the ControlLogix library.

Can anyone help me to resolve this?

Thanks

I have specified the problem in the red comments of the worksheet i am uploading:

Download possible_string_tools_bug.mw

Can someone explain the content and detail of this message:

com.maplesoft.mathdoc.model.plot.PlotException: Unrecognized option in COLOR: RGBA

Have quadrature detected time domain complex data arrays.  Trying to fft array to get frequency domain arrays.  Arrays returned are always sinusoidal in nature rotating between + and - y axis.  Can not phase correct these sinusoidal frequency arrays.  Any ideas how to pretreat time domain arrays before fft or posttreat frequency domain arrays after fft in order to determine phase correction angles that will allow combination of real and imaginary frequency domain data points that are only in + y axis.

How can I set the alpha for the plot symbols? I would like to add some alpha for blending purposes(will help with the visual in my case).

Idealy I would like to plot a 2d "guassian" fade.

A warm greeting for all

How to import a figure from Mathematica to Maple.

Amr

I have a dependent and independent variables u(x),v(x),w(x),....diff( u(x),x$n)=U(k),diff( v(x),x$n)=V(k)......

Is it possible to create, diff(     , x$n)  is an operator or any differentiable function?

I'm having trouble connecting from Maple on Windows 10 to MSQL Server. I tried Microsoft recommended drivers such as sqljdbc_6.4.0.0, did (as I thought) all required steps. The only invariable result I get is "Cannot load driver". I was wandering if anyone had implemented such a construction. Driver name & version , connection string and Java version would be greatly appreciated. Another option is to have any driver, which connects to any of standard databases (Oracle, MySQL).  The only limitation is- it must be from Windows 7 or 10.
          Thanks.
           A.B.

The results should be the same, right?

Download test.mw

I am trying to compare time taken in minutes for each iterative model(Jacobi, Gauss-Seidel and SOR) to complete, so as to figure out the iterativre model with a faster time of convergence but i don't know the command to initiate.

Hello,

I want to evaluate the change of temperature and energy loss during the flow through an expansion valve.

But the command fsolve this does not work with CoolProp.

The following command is just repeated, but gives no result.

fsolve({ThermophysicalData:-Property("D", "H2", "temperature" = TTT, "pressure" = ppp) = 31.13, ThermophysicalData:-Property("H", "H2", "temperature" = TTT, "pressure" = ppp) = 4.098640000*10^6}, {TTT, ppp})

Regards,

Andreas