Hello all,
I have looked at the int2d and int3d commands for double and triple
integrals, but was wondering if there was a simple method of defining the
input triangles or tetrahedrons (int2d and int3d respectively)?
It would be nice to just specify the function, limits and integration
variable(s
Hello all,
This has been queried previously without an obvious solution.
For example: [output, bOK]=dos('grdinfo Test_10-FAA.grd')
Fails with bOK = F
Adding 'start grdinfo Test_10-FAA.grd gives bOK=T but now dos window
console with the result
The only way to get the out put shown below is to pu
ommand or a program.
>
> On the other hand the expression "single text" is not clear. What is a
> "text" in Scilab? Text is meaningless, in Scilab we have strings, and
> "single" is somewhat vague. Does it mean a one-component string?
>
> Your worka
Hello,
Basic query. I have a simple code that applies the Collatz conjecture
equation (3n+1) by running a function and then runs a loop over the values
stored in prime (the first 8 Prime numbers):
clear
exec('collatz.sci',-1);
prime = primes(20);
for i = 1:length(prime)
[ns, seq]=collatz(pri
2021 at 09:24, Samuel Gougeon wrote:
> Le 15/08/2021 à 09:00, Lester Anderson a écrit :
>
> Hello,
>
> Basic query. I have a simple code that applies the Collatz conjecture
> equation (3n+1) by running a function and then runs a loop over the values
> stored in prime (the
Hi Stefan,
I did try that before, but got an error - "Submatrix incorrectly defined"
Lester
On Sun, 15 Aug 2021 at 10:56, Stefan Du Rietz wrote:
>
>
> On 2021-08-15 09:00, Lester Anderson wrote:
> > Hello,
> >
> > Basic query. I have a simple cod
Hi Stefan,
Thank you for clarifying the meaning of the error message.
Lester
On Sun, 15 Aug 2021 at 11:59, Stefan Du Rietz wrote:
> Hello Lester,
> the problem is that seq in each loop is a vector of increasing length!
>
> Stefan
>
>
> On 2021-08-15 12:05, Lester Anderson
Thanks for all the guidance guys!
Lester
On Sun, 15 Aug 2021 at 15:33, Stefan Du Rietz wrote:
>
>
> On 2021-08-15 15:54, Samuel Gougeon wrote:
> > Le 15/08/2021 à 11:28, Lester Anderson a écrit :
> >> Hello Samuel,
> >>
> >> The size of ns (number
Gougeon wrote:
> Le 15/08/2021 à 11:28, Lester Anderson a écrit :
>
> Hello Samuel,
>
> The size of ns (number of steps) and seq (sequence of values) are variable
> depending on the integer input, and this seems to be one issue.
>
> For this reason, seq must be a list, leading
This is the simplest option, based on primes(35):
[p1,p2,p3,p4,p5,p6,p7,p8,p9,p10,p11]=seq([1:$]);
Extract values based on the indices in the seq list.
On Mon, 16 Aug 2021 at 10:49, Lester Anderson wrote:
> Hi Samuel,
>
> Further to your method, is there a way to vectorise the li
Hello all,
I am interested in the Reimann Zeta function, but it does not appear as an
option directly with Scilab (v6.1.1).
Has anyone already prepared a function for this?
Lester
___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/ma
51, Heinz Nabielek wrote:
> Riemann zeta function !
>
> Georg Friedrich Bernhard Riemann (17 September 1826 – 20 July 1866) was a
> German mathematician who made contributions to analysis, number theory, and
> differential geometry.
>
> Heinz
>
>
>
> > On 29.11.2021, at 0
e :
>
> http://numbers.computation.free.fr/Constants/Miscellaneous/zeta.html
>
> S.
>
> Le 29 nov. 2021 à 09:39, Lester Anderson a écrit :
>
>
> Hello all,
>
> I am interested in the Reimann Zeta function, but it does not appear as an
> option directly with Scilab (v6.1.1).
>
Hello,
A quick query. Have adapted existing Matlab code for Scilan to calculate
Bernoulli numbers using an explicit formula (www.bernoulli.org). The code
works fine for numbers from 0 to 10, but 11 onwards the values are
incorrect:
Not the most efficient code but do need to understand where the e
, Samuel Gougeon wrote:
> Hello Lester,
>
>
> Le 29/12/2021 à 09:00, Lester Anderson a écrit :
>
> Hello,
>
> A quick query. Have adapted existing Matlab code for Scilan to calculate
> Bernoulli numbers using an explicit formula (www.bernoulli.org).
>
>
> Not so ex
Hello,
I am trying to understand how to work the Cauchy integral inputs and
replicate the results of a published example:
.e.g. Compute the integral of e^(z^2) / (z-2) assumes C is closed
(anticlockwise) and z=2 is inside C (a simple circle). The solution should
be 2*pi*i*f(2) = 2*pi*i*e^4
In Sc
he expected result:
>
> --> function y=f(z)
> > y = exp((z.^2))./(z-2)
> > endfunction
>
> --> fz=intl(0, 2*%pi, 2, 1, f,1e-10) // gives round-off error
> fz =
>4.199D-13 + 343.05029i
>
> --> 2*%pi*%i*%e^4
> ans =
>0. + 343.05029i
>
>
the default one,
> > we get the expected result:
> >
> > --> function y=f(z)
> >> y = exp((z.^2))./(z-2)
> >> endfunction
> >
> > --> fz=intl(0, 2*%pi, 2, 1, f,1e-10) // gives round-off error
> > fz =
> > 4.199D-13 + 343.05029i
&
Hello all,
Very basic query but how do you print all the values from a result via
mprintf?
e.g.
fac = factor(1729)
7 13 19
mprintf('Factors of 1729 are: %d\n', fac)
just prints 'Factors of 1729 are: 7
Sorry for the basic question!
Cheers
Lester
___
u
= msprintf("%d\n", fac')'
> > fac =
> >"7" "13" "19"
> > --> mprintf("%s%s", s, strcat(fac, " "))
> > Factors of 1729 are: 7 13 19
> >
> > Regards
> > Stefan
>
Hello all,
I have prepared a code to generate Collatz sequences based on the even
(n/2) and odd ((3n+1)/2) rules.
clear
function [ns, seq] = collatz(n)
// Reference: https://en.wikipedia.org/wiki/Collatz_conjecture
//ns =number of steps; seq=Collatz sequence
seq(1) = n;
// Position an index on t
Hi Stephane,
Many thanks for the code pointers - graphics look fab!
Cheers
Lester
On Thu, 14 Apr 2022 at 15:29, Stéphane Mottelet
wrote:
> You can also enjoy the 3d version:
>
> function X=anglePath3(r, th, ph)
> cumth = cumsum(th);
> cumph = cumsum(ph);
> X = cumsum([0 r.*sin(cump
Hello all,
I have a simple code which computes applies the Chinese Remainder theorem
(for a single variable) given three input congruence equations to solve for
x; not very elegant but it works:
My query, is it possible to make this more generalised and vary the number
of input equations and give
'+string(m2)' ...
'x ≡ '+string(a3)'+' mod '+string(m3)' ...
'The solution for (x): '+string(x)+' mod '+string(M)'], ...
'Solution of 3 congruence equations
Sorted the code with multiple equations. Not too sure how to deal with the
messagebox list of equations.
Still have to do that manually - unless there is a way to automatically
build this?
Can this multiple text string be done via a for loop to build a list ?
Thanks
__
input
values of a and m, it is not generating M_inv(5), or M_inv(2) (2
equations).
I have verified what the results should be in Maxima CAS,
[image: image.png]
Any ideas what is going wrong?
Thanks
On Sun, 1 May 2022 at 05:32, Lester Anderson wrote:
> Many thanks Samuel
>
> On Sat, 30 Ap
As a test, I increased the value n (iterations) and that forced it to
converge on the value 150999.
Sounds like I need a test of convergence.
On Sun, 1 May 2022 at 08:37, Lester Anderson wrote:
> Having an issue with the code:
>
> Tried: x ≡ 2 mod 11; 3 mod 12; 4 mod 13; 5 mod 17;
Hi all,
After a lot of trial and error, I have managed to get a set of functions to
compute the approximations of Riemann's Zeta for negative and positive real
values; values of n > 1e6 seem to give better results:
function zs=zeta_s(z, n)
// Summation loop
zs=1;
if z == 0
zs =
quot; uses vectorisation to improve
speed.
Lester
On Sun, 22 May 2022 at 07:31, Lester Anderson wrote:
> Hi all,
>
> After a lot of trial and error, I have managed to get a set of functions
> to compute the approximations of Riemann's Zeta for negative and positive
> real values;
Hi all,
I just found a useful way of getting netCDF files into Scilab via HDF5
(just found this!):
a=h5open('Test_10-FAA.grd');
h5ls(a)
// "lat" "dataset" // "lon"
"dataset" // "z" "dataset" // "Conventions"
"attribute"// "GMT_version" "at
lon',lon);
h5write(b,'data',data);
h5close(b)
save('Test_output.sod','lon','lat','data')
Saving the output creates a file composed of the three datasets (lon, lat,
data). Plotting the data in Scilab via Sgrayplot(x,y,z) it looks correct,
oriented
Hi all,
I previously queried an issue with the dos function in Scilab 6.1.1 not
working as it should; i.e. dos('grdinfo file.grd') would fail and return F.
However, the workaround is the following: [result, bOK]=dos('exit | grdinfo
file.grd')
This will write the result and return T, but not leav
Hi all,
Is it possible to use one vector to find only those closest values in
another vector with a tolerance limit?
A=[2 3 4 5 6]; // source array
B=[0.25 1.2 2.01 5.2 7.4 9.3 10.6 ]; // target array to find closest values
to A
So in the example above, the closest values would be 2.01 and 5.2 i
not too big you can proceed as follow:
> E=abs(ones(B')*A-B'*ones(A))
>
> [i,j]=ind2sub(find(E<=0.3),size(E));
> Serge
>
> Le 02/06/2022 à 10:17, Lester Anderson a écrit :
>
> Hi all,
>
> Is it possible to use one vector to find only those closest val
Hello all,
I await with interest the new version of Scilab. I only use the 64-bit version
so have no issues with the 32-bit version being cut. Most people will be
working on 64-bit systems anyway, and so makes sense to reduce or cut support
for the 32-bit architecture.
With regard to versionin
Hello,
I am revisiting the Collatz plots and wondered how one can (a) randomise the
polyline colours and (b) colour the polyline segments from first to last
segment given a range of start-to-end colours?
Thanks for any pointers
Lester
This email and any attachments are intended solely for th
Hello,
A quick query, how can I make the 's' in (Rs) a subscript?
Thanks
Lester
xlabel 'Relative distance in Schwartzschild Radii (Rs)' fontname 4 fontsize 4
This email and any attachments are intended solely for the use of the
individual or entity to whom it is addressed and may be confiden
Hello,
I have a small FITS file that is 4D and was wondering if this can be loaded
into Scilab and saved out as a csv file?
The header looks like this:
SIMPLE =T /file does conform to FITS standard
BITPIX = -32 /bits per data value
NAXIS =
Hello,
I have made a function to approximate the Riemann Zeta Function (see attached),
but complex input only works by directly calling zeta_3(s,n).
Calling zeta((0.5 + 3*%i),1e6) gives an error message: "Complex comparison not
supported. Please define %s_2_s_custom() or check your code."
For e
precise.
function zz=zeta(s, n)
if isreal(s)==%F
zz = zeta_3(s,n)
elseif s>0 & s<1
zz = zeta_3(s,n)
elseif s==0
zz=-0.5
elseif s==1
zz=%inf
elseif s < 0
zz = zeta_2(s,n)
else s > 1
zz = zeta_1(s,n)
end
endfunction
Lester
On Mon, 17 Jul 2023 at 11:41,
2011.
On Sat, 9 Sept 2023 at 15:03, Lester Anderson
mailto:arctica1...@gmail.com>> wrote:
Hello all,
I have a plot of an Equal-area stereonet with small circles and great circles
plotted. However, all the polylines are stored as compound objects (in this
case 38); is there a way to
Hello all,
I have a plot of an Equal-area stereonet with small circles and great circles
plotted. However, all the polylines are stored as compound objects (in this
case 38); is there a way to unglue all compound objects in one go? Unglue seems
to just do one at a time.
Also, any pointers on h
Hello,
I am trying to loop through a vector of strings but having issues with the
strsplit function.
The vector consists of HEX values for colours. I created a simple function to
get RGB values.
function rgb=hex2rgb(x)
// x=string input
hexrgb=strsplit(x,[2 4])';
rgb=hex2dec(hexrgb
));
blue(i)=hex2dec(rgb(i,3));
RGB_colours=cat(2,red,green,blue);
end
On Tue, 3 Oct 2023 at 12:42, Lester Anderson
mailto:arctica1...@gmail.com>> wrote:
Hello
Thanks for the suggestions. The 'part' option worked but the function fails at
hexrgb=strsplit(x(i),[2 4])';
strsplit:
d
endfunction
You can also use the part function which works with vectors
rgb=[part(x,1:2),part(x,3:4),part(x,5:6)];
Le 02/10/2023 à 09:52, Lester Anderson a écrit :
Hello,
I am trying to loop through a vector of strings but having issues with the
strsplit function.
The vector consists of HEX va
Hello,
I have code for the periodic table of elements (all 118) and selection is via
x_choices and a list. Is there any way to control the strat-up height of the
selection box as at the moment, it stretches top to bottom of the screen and
you have to manually shrink it to get to the ok/cancel b
ons", "labels"], 1:1:10, string(1:1:10));
// no sub-tick graduations
f.sub_ticks=[0,1]
On Sat, 11 Nov 2023 at 13:22, Lester Anderson
mailto:arctica1...@gmail.com>> wrote:
As a quick update:
I have tried the following:
f.auto_ticks(1)='off'
f.x_ticks.locations=[1;2;
Hello
A quick query. I have an axis with a range from 1-10 but the plot has ticks and
labels at 0.5 intervals (1.0, 1.5, 2.0, 2.5 etc). How can I set the interval
and ticks at 1, 2 3 4, etc ?
In the figure properties there is an option to set subticks =0 and step by=1
which does the job. Howev
cks=[0,1]
But get an error "labels should be a string vector." referenced at
x_ticks.locations. Labels are a string vector?
Where is the "Step by" option available in figure edit?
On Fri, 10 Nov 2023 at 09:14, Lester Anderson
mailto:arctica1...@gmail.com>> wrote:
Hello
A
Hello,
I have tried to convert the following Matlab code as a whole and even a very
basic function, but it fails to generate anything!
function [x0,y0,iout,jout] = intersections(x1,y1,x2,y2,robust)
%INTERSECTIONS Intersections of curves.
% Computes the (x,y) locations where two curves intersec
Hello,
Is it possible to make a polygon fill say 50% transparent in a plot? It is easy
to just leave it as an outline poly.The code plots a stereographic net and two
planes with a vector line marking the line of intersection.
Any suggestions would be helpful.
Thanks
Lester
This email and
As Stephane suggested I exported as SVG and used Inkscape to change the opacity.
Lester
On Mon, 26 Feb 2024 at 15:20, Lester Anderson
mailto:arctica1...@gmail.com>> wrote:
Hello,
Is it possible to make a polygon fill say 50% transparent in a plot? It is easy
to just leave it as an o
Hello,
I am looking to create a hillshade function, but issues arise defining the the
X and Y gradients.
The test code so far:
function [grad_x, grad_y] = compute_gradients(dem)
grad_x = diff(dem, 1, 'c');
grad_y = diff(dem, 1, 'r');
endfunction
// Create a test DEM grid
dem=[1.23 1.4
Hello,
I have a plot which forms the basis of an earthquake focal mechanism solution
(attached). The code is adapted from the Matlab BB.m code.
What I am looking to do is clip the curves along segments (A-B), (B-C) and
(A-C), concatenate the x and y values and fill the polygon using xfpoly..
I
Following on from Stephane's suggestion I tried the code suggested but only
partly solved the problem:
Code is attached with the program data. The discourse site did not allow CSV
uploads.
Lester
On Wed, 15 May 2024 at 11:42, Lester Anderson
mailto:arctica1...@gmail.com>> wrote
Data upload as a zip file
On Fri, 17 May 2024 at 11:57, Lester Anderson
mailto:arctica1...@gmail.com>> wrote:
Following on from Stephane's suggestion I tried the code suggested but only
partly solved the problem:
Code is attached with the program data. The discourse site did no
Hello,
Given the attached plot example, is it possible to contour or use Sgrayplot on
such data if it is not on a regular grid?
The plot is a Kaverina earthquake analysis, and I am looking to generate a
density plot within the Kaverina plot region only.All arrays have the same
size, X,Y,data1,
Hello,
I have a simple contour2d plot that contours the zero level in a dataset, and
that is fine. However, it is not easy to extract the data for the two curves as
the figure properties suggest that the data are duplicated.
Checking the graphic entities:
f=gca().children.children
f = [12x1
A workaround finds the unique segments from the plot:
contour2d(x, y, uz, [0 0]);
f = gca().children.children
polyline_idx = find(f.type == 'Polyline');
curve1 = f(polyline_idx(1)).data;
curve2 = f(polyline_idx(2)).data;
curve3 = f(polyline_idx(3)).data;
On Tue, 4 Jun 2024 at 08:
Hello,
I have converted the Matlab code bb.m into Scilab and the calculations are all
fine. Verified the results in Octave. However, the beachball fill plot is not
working:
The Octave/Matlab plot is fine, but the Scilab plot is not filled properly
suggesting that xfpoly is reading the vertex o
Hello,
Just a quick query. Not been using Scilab for a while, so maybe I have
missed something, but is there a built-in precision setting for working
with real and imaginary functions?
Tried the Euler formula (special case) in Scilab 5.4.1:
-->x=cos(%pi)+(%i*sin(%pi))
x =
- 1. + 1.225D-16i
Thanks for the pointer - knew it was a simple one!
Cheers
On 19 February 2014 14:25, Dang, Christophe wrote:
> Hello,
>
> > De la part de Lester Anderson
> > Envoyé : mercredi 19 février 2014 15:16
> >
> > Is there a way to tell Scilab if this gets beyond a limi
It would be helpful to see the script that you are running to understand
the problem, and a screen-capture to show the plot issue.
Lester
On 9 March 2014 22:52, Yavuz Can wrote:
> Hallo together,
>
> I´ve installed the last Version of Scilap. I tried to plot a simple graph,
> but in the opene
Just checked on my Scilab (5.4.1) Windows 7 64-bit
Same issues as you quoted:
plot2d("nl", [1:2], [%nan;1]);// causes crash (or infinite loop)
plot2d("nl", [1:2], [1; %nan]); // no problem - plot drawn
Lester
On 12 March 2014 10:39, Antoine Monmayrant wrote:
> On 03/12/2014 11:03 AM, Iai Masaf
A quick query.
Been a while since I have used Scilab now, but I need to read in a column
data set and process this, and then write out the results as the input
column plus the computed value.
Have an array of the form:
1.0
2.0
3.0
.
.
.
.
.
value(n)
e.g.
read (input.txt) --> array(z) i.e. read
[mailto:users-boun...@lists.scilab.org] *De la part de*Lester
> Anderson
> *Envoyé :* vendredi 25 avril 2014 18:01
> *À :* International users mailing list for Scilab.
> *Objet :* [Scilab-users] Reading single column data for processing
>
>
>
> A quick query.
>
> Bee
Hello,
Is there a simple way to install modules for use in Scilab 5 when no
direct Internet access. Having dowloaded a module zip file - what is
the procedure to install?
Thanks
Lester
___
users mailing list
users@lists.scilab.org
http://lists.scilab.or
/toolboxes/swt - says
it is being "packaged". The zip file contains all the source code and
the builder.sce etc. (I asume that is the reason.)
Lester
On 30 November 2015 at 16:25, Antoine Monmayrant
wrote:
> Le 11/30/2015 04:41 PM, Lester Anderson a écrit :
>>
>> Hello,
>>
Hello,
A basic query, but how can I make sure that the default starting
directory is set to say c:\programs\scilab_work ? Every time Scilab
starts it is deafulted to c:\windows\system32\ - there is no
preference option to set this in 5.4.1 (windows 64-bit).
I can change to the working directory
Hello,
I have the windows version of Scilab (currently 5.4.1), but would like
to run an external program (e.g. gmt software) and return or use the
output in the Scilab code:
As a rough algorithm example (.sce):
//run the gmt program grdinfo.exe - extract data as single column z-matrix)
grd2xyz
first
>
> complex_command = msprintf(...)
>
> and then issue it using the string variable
>
> dos(complex_command,'-echo')
>
> The '-echo' is optional but nice to see in the command window.
>
> The resulting 'faa.txt' file from GMT can
Thanks Samuel for the clarification. Adding the double quotes worked
fine. Will this work with single quotes in v6 as in the original
syntax?
On 3 December 2015 at 18:14, Samuel Gougeon wrote:
> Le 03/12/2015 17:47, Dang Ngoc Chan, Christophe a écrit :
>>
>> Hello,
>>
>
Hello,
A basic query, but I cannot seem to get the Scilab code to run as it
should and plot the data. Not sure where it is falling over. Attached
the code with the original Matlab commands at the front for reference.
Oddly, the straight Matlab code will run and make the plot after first
startingS
Hello,
I am having issues getting the font size to change for just the
legend. I have managed to alter the size for the axis and title. Also,
should both legend and legends work for the same task? I tried to use
legends and set the opt=4 (lower right position) and font_size=4, but
it did not seem
t methods of doing this in
Scilab, but the code does what it should. It does seem that the
formatting is not consistent between axis labels and legends.
Thanks
On 11 January 2016 at 21:24, Samuel Gougeon wrote:
> Hello,
>
> Le 11/01/2016 17:31, Lester Anderson a écrit :
>>
>>
Hello,
I currently have v5.4.1 of Scilab, and have to install modules manually.
I am trying to install apifun (as quite a few modules use it), but get
this error
-->atomsInstall('c:\ss20-work\scilab\apifun_0.4.2-3.bin.zip')
!--error 21
Invalid index.
at line 170 of function atomsInstall cal
Just installed v5.5.2 and get the same issue when trying to install a
module manually, error 21 at line 170 of atomsInstall
Scilab is installed in Program Files on Windows, and so I am guessing
that it might be a permissions issue for a work machine?
On 13 January 2016 at 08:51, Lester Anderson
January 2016 at 06:37, Lester Anderson wrote:
> Just installed v5.5.2 and get the same issue when trying to install a
> module manually, error 21 at line 170 of atomsInstall
>
> Scilab is installed in Program Files on Windows, and so I am guessing
> that it might be a permissions i
;proxyUser","myUserName")
> atomsSetConfig("proxyPassword","myPassword")
>
>
> Read the more documentation below:
> https://wiki.scilab.org/ATOMS
>
> Reagrds,
>
> Lavita Pereira
>
>
>
> On Fri, Jan 15, 2016 at 3:59 PM, Lester Anderson
ly after one has cached module list.
Thanks
On 15 January 2016 at 13:13, wrote:
> Le 01/15/2016 01:45 PM, Lester Anderson a écrit :
>>
>> Hello,
>>
>> Tried setting the proxy with the settings our IT guy gave in the
>> Scilab preferences and it still fail
tried atomsSetConfig('Offline','True') with no change in result
As you say Antoine, room for improvement in terms of installing
offline for sure.
On 15 January 2016 at 13:29, Stéphane Mottelet wrote:
> Le 15/01/2016 14:25, Lester Anderson a écrit :
>>
>> Hell
Still not working.
On 15 January 2016 at 14:26, Antoine Monmayrant
wrote:
> Le 01/15/2016 03:16 PM, Lester Anderson a écrit :
>>
>> tried atomsSetConfig('Offline','True') with no change in result
>
>
> That's weird because I found the rel
?
>
> Last resort:
> edit atomsInstall
> change the file/function names to "MyatomsInstall"
> comment out the lines 221->229
> save
> F5
> MyatomsInstall('mylocalpackageIcannotinstall.zip');
>
>
> Antoine
>
>>
>> On 15 January 2016 at 14:26, Antoine Mon
the end :) I will also let our
IT guys know.
On 15 January 2016 at 15:49, Antoine Monmayrant
wrote:
> Le 01/15/2016 04:26 PM, Lester Anderson a écrit :
>>
>> Hi Antoine,
>>
>> Yes all is set fine:
>>
>> -->atomsSetConfig('offline','True'
Hello,
I am trying to create a loop to plot three curves from an input variable (Te):
clear()
rho_m = 3330;
rho_fill = 2400;
g = 9.81;
Pb = -1.5e12;
v = 0.25;
E = 1e11;
x=linspace(0,200,41);
for Te = 5000:25000:1 // start:end:step to make 5000, 15000, 25000
D = E*Te .^3/(12*(1-v^2));
Thanks for that pointer.
The plot works but only does Te=25 in this case
for Te = 5000:1:25000 // start: step: end
D = E*Te .^3/(12*(1-v^2));
lamda = ((rho_m-rho_fill)*g ./(4*D)).^0.25;
lamda_k = lamda*1000;
flex = 2*Pb*lamda/((rho_m-rho_fill)*g)*exp(-lamda_k*x).*cos(lamda_k*x
> flex =
> [flex,2*Pb*lamda/((rho_m-rho_fill)*g)*exp(-lamda_k*x).*cos(lamda_k*x)];
> end
> plot(x,flex)
>
>
>
> Le 18/01/2016 11:50, Lester Anderson a écrit :
>>
>> Thanks for that pointer.
>>
>> The plot works but only does Te=25 in this case
Hello,
I have a quick query realted to the change from 5.4.1 to 5.5.2 (and
6), related to a vector^scalar operation. The basic code runs and
plots the curves as needed but generates warnings on running:
clear()
clc()
functions_library = lib('C:\geophys\Lester\SCILAB_work\Functions');
//Finds the
Hello all
I am looking to use the sciNetCDF module to work with GMT grid files
but have some issues.
Reading a GMT grid is not a problem as far as I can tell, so the
following does read the file (and z variable), storing the values in
the grav array:
source = fullfile('Q:\geophys\Gz.grd')
grav =
Hello
A quick query. How would one define the Hilbert transform of a grid
for X and Y directions; looking for two solutions Hx and Hy (for the
real values).
Thanks
Lester
___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/list
ncwrite by typing help ncwrite in the console.
>
> Regards
> Harpreet
>
> On Fri, Feb 5, 2016 at 3:18 PM, Lester Anderson
> wrote:
>>
>> Hello all
>>
>> I am looking to use the sciNetCDF module to work with GMT grid files
>> but have some issues.
>&g
;Hilb_Hy', 'Dimensions', list('r', no_row, 'c', no_column); after that you
> can pass the ncwrite command and you have to pass like this
> ncwrite(filename, 'Hilb_Hy',Hy) as Hy is a matrix.
>
> Regards
> Harpreet Singh
>
> On Fri,
Hi Serge,
I am working with grid data, so looking for the 2D Hilbert, the
results I think should appear similar to doing a directional
derivative, where X highlights features with N-S trends and E-W for Y.
Lester
On 5 February 2016 at 12:24, Serge Steer wrote:
> Le 05/02/2016 10:56, Les
not
> create it twice.
>
> Delete the file, and try again.
>
>
> An example of grid writing, taken from the ncwrite() help page:
>
> nccreate(filename, 'mat', 'Dimensions', list('r', 3, 'c', 4));
> ncwrite(filename, 'mat',
lumn and row-by-row. That SHOULD
> work, and if you do it as matrix operations it should be pretty fast in
> Scilab.
>
> On Fri, 2016-02-05 at 12:52 +, Lester Anderson wrote:
>> Hi Serge,
>>
>> I am working with grid data, so looking for the 2D Hilbert, the
>&g
Hi all,
A quick one. I am looking for the Scilab equivalent to padarray in Matlab.
In the Matlab code the line is:
datapad = padarray(data, [rdiff cdiff], 'replicate') where cdiff/rdiff
are column/row, with the replicate function padding with the values at
the borders of "data"
The closest I can
Thanks for the ideas, it kind of solves the need, but pads to left and
bottom only:
adding 2 cols/rows say to [1 2 3; 4 5 6; 7 8 9]
1 2 3 3 3
4 5 6 6 6
7 8 9 9 9
7 8 9 9 9
7 8 8 9 9
Is it possible to replicate on the top and right edges? (1 2 3) and (1 4 7)
Also looked at resize_matrix(data, nr,
Thanks for the update Samuel. Got the padding working as you suggested:
clear()
data = [1 2 3; 4 5 6; 7 8 9]
addedRows = 3; addedCols = 3;
tmp = [data ; data($, :) .*. ones(addedRows, 1)];
tmp = [tmp tmp(:, $) .*. ones(1,addedCols)];
tmp = [repmat(tmp(1, : ), addedRows, 1) ; tmp];
paddedMat = [r
Hello
At the moment the only way I can install additional modules is in
"offline" mode.
Tried all the usual routes for proxy settings but the program will not
find the default proxy in place, even though other applications will
do (without manually setting the values). Tested with the free GIS
pr
Hello,
Quick query. If I have a data file containg three columns e.g.
Long Lat gravity
12 2 34.4
20 4 60.3
30 6 23.7
I want to read column 3 (gravity) to process and then write the
results back out as:
gravity_processed = f(gravity)
Long Lat gravity_processed
12 2 g1
20 4 g2
30 6 g3 etc.
Is th
1 - 100 of 134 matches
Mail list logo