Hi!

This question is related to http://www.mapleprimes.com/questions/204419-Derivatives-Of-Splines-Are-Not-Defined and http://www.mapleprimes.com/questions/42114-Problem-With-Spline-Integrating , however I have not been able to apply the solutions given there to my problem.

I have a set of points given by

and certain function value points given by

where e1 is the function I am approximating.

Using

I come up with my piecewise function.

When I do diff(e4,x), however, the points at the nodes show "float(undefined) x=0.2..."(the node).

As it turns out, the value of the derivative on the left of the node is not equal to that of the right side by a factor of 10^(-7), in other words, numerically unimportant but high enough for maple to realise it is not the same number. How could I tell Maple that I am happy choosing, for example, the value given by the function on right side of the node?

I welcome any suggestions.

Many thanks in advance.

Every time I try to write a procedure I get stuck.

This time is no different:

restart;

global a:=0.081819221, PI:=3.1415926535897932384626433832795;
Ecce:=proc(lt)
lat:=lt*(PI/180);
b:=((1-a*sin(lat))/(1+a*sin(lat)))^(a/2));
t:=ln((tan(PI/4)+lat/2));
d:=3437.7468*t*b;
return d;
end proc;

I digit the following to get a result
Ecce(45.2112);

and this is what I get (in blue)
Ecce(45.2112)

Every single time. I can never have a procedure that works right away. It's getting on my nerves

I want to extract all the symbols and last trade from 

http://finance.yahoo.com/q/cp?s=%5EIXIC

I have tried:

status, data, headers := HTTP:-Get("http://finance.yahoo.com/q/cp?s=%5EIXIC");
data;

but it just gives me a bunch of jiberish! Any ideas?

please i would love make a 3d plot of this Bessel function expression against the parameter p and T[0] using maple (p^(((2+deltaT[0])/(2)))*(C[1]*BesselJ((delta)/(-2)*sqrt(T[0]-4 T[g]),(-((T[g])/(T^(2)[0])*gamma*G*delta*sqrt(alpha*k))/(p))+C[2]*BesselY((delta)/(-2)*sqrt(T[0]-4 T[g]),(-((T[g])/(T^(2)[0])*gamma*G*delta*sqrt(alpha*k))/(p)) where [delta = 0.2e-2, T[g] = (1.1, .615, .48, .2962), k = (1.2, 1.3, 1.4, 1.5), G = .2, gamma = .2, alpha = 5.36, C[1] = 500, C[2] = 100,p = (.22, .23, .24, .25), T[0] = (3.666667, 2.307692, 1.714286, 1.377778)

we can get the continued fraction of cot(x) with the command

> convert(cot(x), confrac, x, 6)

the result is

(1+x^2/(-3+x^2/(5-(1/7)*x^2)))/x

how could I convert this result to a form like

1/x-(1/10)*x+49*x/(20*x^2-210)

In this work the theme of vector analysis shown from a computational point of view; this being a very important role in the engineering component; in civil and mechanical special it is why, using the scientific software Maple develops interactive solutions for long processes through MapleCloud calculations. At present the majority of professors / researchers perform static classes open source leaves; so that our students learn and memorize commands, thus generating more time learning in the area. Loading Bookseller VectorCalculus develop topics: vector algebra, differential operators, conservative fields, etc. Maplesoft making processes provide immediate calculations long operation Embedded Components displayed in line with MapleNet integrations. Today our future engineers to design solutions and will be launched in the cloud thus being a process with global qualification in the specialty. Significantly Maple is a scientific software which allows the researcher to design their own innovations and not use themes for their manufacturers.

III_CRF_2015.pdf

CRF_2015.mw

L.AraujoC.

I have some excel data which I need to fit in the formula:

0.5a*erfc(0.5*2^0.5*(-x+m1/s1)) + (0.5-0.5a)*erfc(0.5*2^0.5*(-x+m2/s2))

In this formula the coefficient m ans s are mean and standard deviation and a is the weigth of 2 peak in a cummulative gaussian distribution.

I fitted this (I will insert the maple file) and found:

-48736.43 erfc(-0.707x + 0.707) + 48736.43 erfc(0.707x + 0.707)

Can this be correct?

And how do I find the values of m1, s1, m2, s2 and a?

data.xlsx

data_fit.mw

I want to solve numerically the nonlinear pde:

u_x+u_t - (u_{xt})^2 = u(x,t)

which method do you propose me to use with maple? (I don't mine about which boundary conditions to be used here).

Hello.
Have a question for use Physics[TransformCoordinates] . For example, in a Cartesian coordinate system is an arbitrary tensor. As you know, in this case the covariant, contravariant and mixed components are the same. We have a coordinate transformation such as transforming our Cartesian coordinate system in the oblique coordinate system. In this simple example, correctly calculated the covariant, contravariant and mixed components (structure: covariant and contravariant), but other mixed components (structure: contravariant and covariant) are calculated is not correct. I checked by hand.

Download Transformation_tensor_components.mw

I recently got myself a 64 bit computer and have noticed that I cannot use the option compile=true in dsolve/numeric. Take the following simple example:

dsolve({diff(x(t),t)=x(t),x(0)=1},numeric,compile=true);
Error, (in dsolve/numeric/SC/preproc) unable to compile (rc=1), please try again, and if that fails verify your Windows compiler installation

I'm using Windows 10, but had the same problem with Windows 8.1 on the same machine.
The Compiler:-Compile examples in the help page all work.

What do I have to do to make the option compile=true work in dsolve/numeric?

You may safely assume that I don't know any technicalities about these things.

hi, I am tyying to solve this equation but there is arising an error, plz help me,
VIM.3rd_order.mw

How to change the font style, globally, the inputs and outputs of computations in Maple 2015? That is, every time a new document is started, the chosen style is automatically applied.

Oliveira.

This Question involves using dsolve(..., numeric) for an IVP specifed by a procedure. This is based on a Question asked earlier today. In this Question, I have no interest in how to solve this IVP or in why this solution technique fails. In the worksheet below, the odeplot command seems to get stuck in an infinite loop (I am not interested in why that happens), and I press the stop button (in the Standard GUI). Then, instead of the usual Warning, computation interupted message followed by a return to the command prompt, I get an informative message and the plot that has been computed so far. This seems like a very useful feature: to return the results computed so far after an interuption. Furthermore, those results are programmatically accessible. My Question is How is this done? How do you trap the stop button and return the results?

Download trap_stop_button.mw

restart;
odes:=diff(x1(t),t)*diff(x2(t),t$2)*sin(x1(t)*x2(t))+5*diff(x1(t),t$2)*diff(x2(t),t)*cos(x1(t)^2)+t^2*x1(t)*x2(t)^2=exp(-x2(t)^2),diff(x1(t),t$2)*x2(t)+diff(x2(t),t$2)*diff(x1(t),t)*sin(x1(t)^2)+cos(diff(x2(t),t$2)*x2(t))=sin(t);
ics:=x1(0)=1,D(x1)(0)=1,x2(0)=2,D(x2)(0)=2;
subs(diff(x1(t),t$2)=yp2,diff(x2(t),t$2)=yp4,diff(x1(t),t)=Y[2],diff(x2(t),t)=Y[4],x1(t)=Y[1],x2(t)=Y[3],{odes});
p:=proc(N,t,Y,YP)
local eqs,yp2,yp4;
YP[1]=Y[2];
YP[3]=Y[4];
eqs:=[yp2*Y[3]+yp4*Y[2]*sin(Y[1]^2)+cos(yp4*Y[3]) = sin(t), Y[2]*yp4*sin(Y[1]*Y[3])+5*yp2*Y[4]*cos(Y[1]^2)+t^2*Y[1]*Y[3]^2 = exp(-Y[3]^2)];
YP[2],YP[4]:=op(subs(fsolve(eqs,{yp2=1,yp4=2}),[yp2,yp4]));
end proc:
res:=dsolve(numeric,procedure=p,initial=array([1,1,2,2]),number=4,procvars=[x1(t),diff(x1(t),t),x2(t),diff(x2(t),t)],start=0,maxfun=0);
plots:-odeplot(res,[t,x1(t)],0..5,gridlines=true);

It's a long time to wait for the odeplot.

Any advice is appreciated.