@Mac Dude To be clear, I did not originally mean that a custom file would be created anew for each call (although that is an alternate approach which I was considering).

I meant that the single custom file would accept the full set of arguments that you are now passing to ssystem. This file would be written just once. Its first line would be, say,

#!/bin/tcsh

on Unix/Linux (and I don't know what on OSX). It could set its PATH and PYTHONPATH or other environment variables. It would pass along its arguments to a call to python.

Do you have a valid Maple license?

Stop posting duplicates of this question.

If you have additional details, add yhe here as a Reply/Comment on the original Question.

@adel-00 Sorry, I will be travelling for 3 weeks and have no Maple access.

Perhaps next time you could provide the correct expressions up front.

@mmcdara It's great that you found it, Statistics:-Specialize.

You last comment also has a struck-out fragment containing (mu, sigma) -> RandomVariable(Normal(mu, sigma)) . Please note that that is a procedure which generates new RVs from scratch, and not a tool for specializing/instantiating a given, extant RV (which was actually your Question). It might accomplish your goal, but it doesn't do precisely what you asked.

Please add your followup queries on this here, instead of as separate Questions.

@denitsastaicova Now you've changed your request. At first you asked how to use your given z1 to shade via colorscheme with zgradient. And now it seems as if you'd prefer a slightly different scaling and transition for the gradient of values than are in the z1 which you originally supplied.

You've subsequently added the information that what you really want is a way to visually discern the highest z1 values easily. Well, your problem there is in your choice of a zgradient from Black to RoyalBlue. But I didn't suggest that scheme... you did!

But there is another approach which can use any computed z1 in its own "zgradient" colorscheme to generate coloring data. That coloring data can then be extracted and re-used as a coloring replacement in a x1,y1 point-plot. See below for an example.

You can also adjust the scaling in the z1 to get a different pace of color gradient. (I've provided an additional example below.)

Or you can choose different colors for the gradient. Or you can use 1-~z1 to reverse your original color gradient. It's up to you which color gradient pleases you most. What I'm showing is that, yes, a 3rd list like z1 can indeed be used to denote the colorscheme zgradient in a point plot of x1 vs y1. I leave it to you to devise which z1 values and which colors for the gradient you want.

Download color_pointplot2.mw

@Kitonum On the line that assigns to single try using the prefix form for the call to &not, as a function call using brackets.

Error_ac.mw

It works in Maple 2018.2 using

  &not( ... )

and, fwiw,

  NULL &not ...

@mmcdara Unfortunately the approach you've given does not do as well in later versions. You used Maple 2015, it seems.

But with two little tweaks it can still serve. (See my Answer as to why, if you'd like.)

Download int_mmcd_ac.mw

@Carl Love The crash is due to the evalf of these Int instances at higher precision (see my Answer, eg, at Digits=25 or so).

Yes, it's true that it would be better if fsolve didn't try and demand that kind of precision increase. But the bug is in evalf(Int(...)) .

@pauldaas 

It's too bad that you've ignored all my previous suggestions. Your code still uses deprecated linalg and matrix, unnecessary stacking to iteratively build the Matrices, unnecessary globals, and substitution done via unnecessary assignment to the list of key names.

This code is severely flawed. It is unnecessarily fragile and will likely cause future grief.

It's also too bad that you could not properly explain the purpose of the code. Comments about load forces on a beam are irrelevant. What would be useful would be information about why you need to build R row by row, whether you really need C at the end (or just R), whether you know the number of rows in final R, etc.

I would not be surprised if a robust and sensible version of code to accomplish your true goals could be only a few lines long.

It's too bad that you ignored all my suggestions in a previous thread. Your code still uses deprecated linalg and matrix, unnecessary stacking to iteratively build the Matrices, unnecessary globals, and substitution done via unnecessary assignment to the list of key names.

This code is severely flawed. It is unnecessarily fragile and will likely cause future grief.

@dharr Yes, it looks like it may just be reformulated as,

   K1*Int(cos(theta),theta=0..2*Pi) + K2*Int(sin(theta),theta=0..2*Pi)

where K1 and K2 do not depend on theta.

070919_ThetaIntegral_ac.mw

@9009134

For example,

int3_ac.mw

@adel-00 Is this what you are trying to accomplish?

Apart from other problems including misuse of assumed names, your sprintf calls tried to use Lambda as if its value were a float, but it is unassigned.

Download adel00.mw