Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Hi, 

I'm trying to calculate the inner product of these two vectors. 

|a> = (1 + i) i, (1) j, (i) k

|b> = (4 - 1) i, (0) j, (2 -2i)k

where bold i, j, k, are the orthonormal bases. 

I can calculate it on paper, the problem I'm having is figuring out how to calculate it using maple. I want to be able to use dirac notation. I've tried using the physics library but none of the example use actual vectors but instead use quantum numbers. 

Thanks, 

Poiseuill's law: the formula is correct if both sides has the same dimension, so it is required that 

     [v]=[(P*Rm*nn)/L]

find n and m that make the above equation dimensionally correct

It might seem like a small thing, but is there a way to separate the text from the numbers? It's only for optical purposes, as i think it looks better. 

 

Thank you,

MO

Hi everyone, I am trying to dsolve a ode but could not get the answer, could anyone please help me or guide me what's wrong with my worksheet


 

restart

u := -3*beta*(2*p^3*sigma^4*(1/3)+(2*p^3*y^2*(1/3)-4*q*(1/3))*sigma^2-4*p^2*y*(k+1)*sigma*(1/3)+p*(k+1)^2)*(y+sigma)*(y-sigma)/(8*sigma^2)+(-p*sigma^3+(p*y^2-k+1)*sigma-(k+1)*y)/(2*sigma)

-(3/8)*beta*((2/3)*p^3*sigma^4+((2/3)*p^3*y^2-(4/3)*q)*sigma^2-(4/3)*p^2*y*(k+1)*sigma+p*(k+1)^2)*(y+sigma)*(y-sigma)/sigma^2+(1/2)*(-p*sigma^3+(p*y^2-k+1)*sigma-(k+1)*y)/sigma

(1)

ode := diff(theta(y), y, y)+G*(diff(u(y), y))^2+G*beta*(diff(u(y), y))^4

diff(diff(theta(y), y), y)+G*(-(3/8)*(diff(beta(y), y))*((2/3)*p(y)^3*sigma(y)^4+((2/3)*p(y)^3*y(y)^2-(4/3)*q(y))*sigma(y)^2-(4/3)*p(y)^2*y(y)*(k(y)+1)*sigma(y)+p(y)*(k(y)+1)^2)*(y(y)+sigma(y))*(y(y)-sigma(y))/sigma(y)^2-(3/8)*beta(y)*(2*p(y)^2*sigma(y)^4*(diff(p(y), y))+(8/3)*p(y)^3*sigma(y)^3*(diff(sigma(y), y))+(2*p(y)^2*y(y)^2*(diff(p(y), y))+(4/3)*p(y)^3*y(y)*(diff(y(y), y))-(4/3)*(diff(q(y), y)))*sigma(y)^2+2*((2/3)*p(y)^3*y(y)^2-(4/3)*q(y))*sigma(y)*(diff(sigma(y), y))-(8/3)*p(y)*y(y)*(k(y)+1)*sigma(y)*(diff(p(y), y))-(4/3)*p(y)^2*(diff(y(y), y))*(k(y)+1)*sigma(y)-(4/3)*p(y)^2*y(y)*(diff(k(y), y))*sigma(y)-(4/3)*p(y)^2*y(y)*(k(y)+1)*(diff(sigma(y), y))+(diff(p(y), y))*(k(y)+1)^2+2*p(y)*(k(y)+1)*(diff(k(y), y)))*(y(y)+sigma(y))*(y(y)-sigma(y))/sigma(y)^2-(3/8)*beta(y)*((2/3)*p(y)^3*sigma(y)^4+((2/3)*p(y)^3*y(y)^2-(4/3)*q(y))*sigma(y)^2-(4/3)*p(y)^2*y(y)*(k(y)+1)*sigma(y)+p(y)*(k(y)+1)^2)*(diff(y(y), y)+diff(sigma(y), y))*(y(y)-sigma(y))/sigma(y)^2-(3/8)*beta(y)*((2/3)*p(y)^3*sigma(y)^4+((2/3)*p(y)^3*y(y)^2-(4/3)*q(y))*sigma(y)^2-(4/3)*p(y)^2*y(y)*(k(y)+1)*sigma(y)+p(y)*(k(y)+1)^2)*(y(y)+sigma(y))*(diff(y(y), y)-(diff(sigma(y), y)))/sigma(y)^2+(3/4)*beta(y)*((2/3)*p(y)^3*sigma(y)^4+((2/3)*p(y)^3*y(y)^2-(4/3)*q(y))*sigma(y)^2-(4/3)*p(y)^2*y(y)*(k(y)+1)*sigma(y)+p(y)*(k(y)+1)^2)*(y(y)+sigma(y))*(y(y)-sigma(y))*(diff(sigma(y), y))/sigma(y)^3+(1/2)*(-(diff(p(y), y))*sigma(y)^3-3*p(y)*sigma(y)^2*(diff(sigma(y), y))+((diff(p(y), y))*y(y)^2+2*p(y)*y(y)*(diff(y(y), y))-(diff(k(y), y)))*sigma(y)+(p(y)*y(y)^2-k(y)+1)*(diff(sigma(y), y))-(diff(k(y), y))*y(y)-(k(y)+1)*(diff(y(y), y)))/sigma(y)-(1/2)*(-p(y)*sigma(y)^3+(p(y)*y(y)^2-k(y)+1)*sigma(y)-(k(y)+1)*y(y))*(diff(sigma(y), y))/sigma(y)^2)^2+G*beta*(-(3/8)*(diff(beta(y), y))*((2/3)*p(y)^3*sigma(y)^4+((2/3)*p(y)^3*y(y)^2-(4/3)*q(y))*sigma(y)^2-(4/3)*p(y)^2*y(y)*(k(y)+1)*sigma(y)+p(y)*(k(y)+1)^2)*(y(y)+sigma(y))*(y(y)-sigma(y))/sigma(y)^2-(3/8)*beta(y)*(2*p(y)^2*sigma(y)^4*(diff(p(y), y))+(8/3)*p(y)^3*sigma(y)^3*(diff(sigma(y), y))+(2*p(y)^2*y(y)^2*(diff(p(y), y))+(4/3)*p(y)^3*y(y)*(diff(y(y), y))-(4/3)*(diff(q(y), y)))*sigma(y)^2+2*((2/3)*p(y)^3*y(y)^2-(4/3)*q(y))*sigma(y)*(diff(sigma(y), y))-(8/3)*p(y)*y(y)*(k(y)+1)*sigma(y)*(diff(p(y), y))-(4/3)*p(y)^2*(diff(y(y), y))*(k(y)+1)*sigma(y)-(4/3)*p(y)^2*y(y)*(diff(k(y), y))*sigma(y)-(4/3)*p(y)^2*y(y)*(k(y)+1)*(diff(sigma(y), y))+(diff(p(y), y))*(k(y)+1)^2+2*p(y)*(k(y)+1)*(diff(k(y), y)))*(y(y)+sigma(y))*(y(y)-sigma(y))/sigma(y)^2-(3/8)*beta(y)*((2/3)*p(y)^3*sigma(y)^4+((2/3)*p(y)^3*y(y)^2-(4/3)*q(y))*sigma(y)^2-(4/3)*p(y)^2*y(y)*(k(y)+1)*sigma(y)+p(y)*(k(y)+1)^2)*(diff(y(y), y)+diff(sigma(y), y))*(y(y)-sigma(y))/sigma(y)^2-(3/8)*beta(y)*((2/3)*p(y)^3*sigma(y)^4+((2/3)*p(y)^3*y(y)^2-(4/3)*q(y))*sigma(y)^2-(4/3)*p(y)^2*y(y)*(k(y)+1)*sigma(y)+p(y)*(k(y)+1)^2)*(y(y)+sigma(y))*(diff(y(y), y)-(diff(sigma(y), y)))/sigma(y)^2+(3/4)*beta(y)*((2/3)*p(y)^3*sigma(y)^4+((2/3)*p(y)^3*y(y)^2-(4/3)*q(y))*sigma(y)^2-(4/3)*p(y)^2*y(y)*(k(y)+1)*sigma(y)+p(y)*(k(y)+1)^2)*(y(y)+sigma(y))*(y(y)-sigma(y))*(diff(sigma(y), y))/sigma(y)^3+(1/2)*(-(diff(p(y), y))*sigma(y)^3-3*p(y)*sigma(y)^2*(diff(sigma(y), y))+((diff(p(y), y))*y(y)^2+2*p(y)*y(y)*(diff(y(y), y))-(diff(k(y), y)))*sigma(y)+(p(y)*y(y)^2-k(y)+1)*(diff(sigma(y), y))-(diff(k(y), y))*y(y)-(k(y)+1)*(diff(y(y), y)))/sigma(y)-(1/2)*(-p(y)*sigma(y)^3+(p(y)*y(y)^2-k(y)+1)*sigma(y)-(k(y)+1)*y(y))*(diff(sigma(y), y))/sigma(y)^2)^4

(2)

dsolve(ode)

Error, (in dsolve) Required a specification of the indeterminate function

 

bc := theta(-sigma) = 0, theta(sigma) = 1

theta(-sigma) = 0, theta(sigma) = 1

(3)

sol := dsolve({bc, ode})

(4)

``


 

Download help_ode.mw

When I type the command

sum(2^n, n=0..N)

Maple calculates it correctly. However, if I then go back to that expression and click on it to change the value of N, Maple tries to execute something before I hit enter. It won't stop unless I press the Interrupt button, and after that I have to insert some other calculation, e.g. 1+1, for it to work once again.

This happens even on completely new files and for any explicit values of N.

Why is this happening?

Hello! I am trying to plot a parabola given by the equation x=y^2. I used the following code

plot(y^2, y = -1 .. 1, color = blue, thickness = 3, title = "T1.1 x=y^2 [ECE]")

but my plot looks like this

I want it to have the x-axis horizontally displayed and the y-axis to be vertically displayed. Something like this

File:X=y^2.svg - Wikimedia Commons

Any advice? 

Thank you!

I think I am doing this question right. Is there anything to improve the code? 

How can I get exactly graph plotting as the sample show? All blue area under the curve. 

 


 

" ref: https://www*maplesoft*com/applications/view*aspx"?SID=153480&view=html"  Data :=  [ [0,0.0,0.0],      [5,0.2,0.1],      [10,0.25,0.2],      [18,0.3,0.3],      [25,0.4,0.4],      [32,0.6,0.45],      [38,0.72,0.50],      [43,0.6,0.40],      [47,0.3,0.30],      [50,0.0,0.0] ]: xdata := [0,5,10,18,25,32,38,43,47,50];  ddata := [0.0,0.2,0.25,0.3,0.4,0.6,0.72,0.6,0.3,0.0];  vdata := [0.0,0.1,0.2,0.3,0.4,0.45,0.5,0.4,0.3,0.0];  "

[0, 5, 10, 18, 25, 32, 38, 43, 47, 50]

 

[0., .2, .25, .3, .4, .6, .72, .6, .3, 0.]

 

[0., .1, .2, .3, .4, .45, .5, .4, .3, 0.]

(1)

with(plots); with(plottools); n := nops(xdata)-2; riverbed := plot([seq([xdata[i], -ddata[i]], i = 1 .. n+2)], style = line, color = blue); readings := seq(plot([[xdata[i], 0], [xdata[i], -ddata[i]]], style = line), i = 1 .. n+2); for i to n do x := xdata[i+1]; d := ddata[i+1]; v := vdata[i+1]; s := sprintf("%4.2f m/s", v); xplot[i] := textplot([x, 0, cat("x", i)], align = above); vplot[i] := textplot([x, -d-0.5e-1, s], font = [HELVETICA, 14]) end do; xtickmarks := seq(xplot[i], i = 1 .. n); velocities := seq(vplot[i], i = 1 .. n); display([readings, riverbed, velocities, xtickmarks], view = [0 .. 51, -1 .. 0], labels = ["river cross section x (m)", "depth (m)"]); n := nops(xdata); for i from 2 to n-1 do x := xdata[i]; d := ddata[i]; xL := xdata[i-1]; xR := xdata[i+1]; poly[i] := polygon([[(xL+x)*(1/2), 0], [(xL+x)*(1/2), -d], [(xR+x)*(1/2), -d], [(xR+x)*(1/2), 0]], color = cyan) end do; polys := seq(poly[i], i = 2 .. n-1); plots[display]([readings, riverbed, polys, velocities, xtickmarks], view = [0 .. 51, -1 .. 0], labels = ["river cross section x (m)", "depth (m)"])

``


 

Download Untitled.mw

When Maple converts sin(x)^n to Latex, the result remain  sin(x)^n.  But in Mathematical typesetting, this is normally written as sin^n(x).   Ofcourse this is only for Latex. In Maple code this not valid.

Is it possible to change Maple's Latex to make it do this automatically? Mathematica does this automatically. Here is an example

restart;
expr:=sin(x)^3+cos(3*x)^5;              
Physics:-Latex(expr)
 
            \sin \left(x \right)^{3}+\cos \left(3 x \right)^{5}

Which when compiled gives

Compare to Latex generated by Mathematica

Which compiles to 

Which is more standard in books and papers, than Maple's version.

Both Maple's Physics:-Latex and latex() command do the same thing.

Is there a way to make it generate the improved version for latex?

Maple 2020.1

 

 

sometimes I get intermediate expressions generated from other operations that contain terms such as exp(x)^n in them. As an example, exp(x)^3.  In Mathematica, it automatically replaces these by exp(3*x). But in Maple I need to force this change.

For purposes of Latex only, I like to change these terms to exp(3*x) before converting the whole expression to Latex, as it is looks much better that way.

expr:=exp(x)^3;
Physics:-Latex(expr)

                \left({\rm e}^{x}\right)^{3}

expr:=exp(3*x);
Physics:-Latex(expr)

                {\rm e}^{3 x}


I found that doing simplify(expr,exp)  does the trick. It changes exp(x)^n to exp(n*x). But I am worried about applying this whole simplification command to the whole expression, which can be very large, and do not want to change it all yet.

I just want to change any occurance of exp() there, and nothing more.

I tried using subsindent to do that, but it does not work on terms in denominator

restart;
expr:=exp(x)^3*sin(x)+3/(exp(x)^n);
subsindets(expr,'exp(anything)^anything',f->simplify(f,exp))

I tried

subsindets(expr,'1/exp(anything)^anything',f->simplify(f,exp))

and it did not work.

I am still not good at subsindent. How to make it change all exp(x)^n to exp(n*x) everywhere?

Hello

I need to detect if different expressions contain a radical. For example

aaa := X3*(alpha[2, 8]*(sqrt(X2/alpha[1, 7])*alpha[3, 8] + X1*alpha[3, 6])*X2 + X3*(alpha[2, 8] + alpha[3, 9]/2))/(X2*alpha[2, 8])

As can be noticed there is a square root in the expression.   

I have tried type(expr,sqrt) and has but to no avail. (I am not sure if I use them as they should though).   

Many thanks

 

Ed

 

Hi:

Every once-in-a-while Maple crashes due to a bug in numerical integration that I reported ages ago, but was never fixed (or if it was, I was never told). Anyway, this time, in addition to crashing, it wiped out my worksheet, and froze the machine. When I recovered, Maple had lost all memory of my preferences, recent files, etc - it was as if I had a clean installation.

I restored the link to the usual style file I use, but it still refuses to give me back the style I want - that is - the classical worksheet. Every command I enter now requires F5 or I will get 2D input which I hate. Does anyone know where the preferences, probably an .INI file, are kept, so I can restore it from an older version of Maple.

Thank you

Hi, i am using solve and solve command to find the root but when i used fsolve command to separate only real root, could separate all roots, can anyone correct me, please


 

restart

f := 9.765625000*10^(-6)*(-6671.221362*(x^2+2)^5*sqrt(2)*arctan((1/2)*x*sqrt(2))*x-555.9351135*(x^2+2)^6/((1/2)*x^2+1)-10479.13001*(x^2+2)^5*sqrt(2)*x-(374220*(0.297116730e-1*x^9+.269385824*x^7+.99643086*x^5+5.18951288*x^3+4.42867382*x))*x-1111.870227*x^10-12601.19538*x^8-62147.39274*x^6-485504.8775*x^4-828649.1585*x^2-788850.2769)/(x^2+2)^6-(0.1171875000e-3*(-555.9351135*(x^2+2)^6*sqrt(2)*arctan((1/2)*x*sqrt(2))-873.2608343*(x^2+2)^6*sqrt(2)-(374220*(0.29711673e-2*x^10+0.33673228e-1*x^8+.16607181*x^6+1.29737822*x^4+2.21433691*x^2+2.107985348))*x))*x/(x^2+2)^7+(3.484800000*sqrt(2)*(x^2+2)*arctan((1/2)*x*sqrt(2))*x+.8712000000*(x^2+2)^2/((1/2)*x^2+1)+(5.473911040*(x^2+2))*sqrt(2)*x+5.227200000*x^2-22.99200001)/(16*(x^2+2)^2)-(.8712000000*sqrt(2)*(x^2+2)^2*arctan((1/2)*x*sqrt(2))+1.368477760*sqrt(2)*(x^2+2)^2-36*x*(-0.484000000e-1*x^2+.638666667))*x/(4*(x^2+2)^3)

0.9765625000e-5*(-6671.221362*(x^2+2)^5*2^(1/2)*arctan((1/2)*x*2^(1/2))*x-555.9351135*(x^2+2)^6/((1/2)*x^2+1)-10479.13001*(x^2+2)^5*2^(1/2)*x-374220*(0.297116730e-1*x^9+.269385824*x^7+.99643086*x^5+5.18951288*x^3+4.42867382*x)*x-1111.870227*x^10-12601.19538*x^8-62147.39274*x^6-485504.8775*x^4-828649.1585*x^2-788850.2769)/(x^2+2)^6-0.1171875000e-3*(-555.9351135*(x^2+2)^6*2^(1/2)*arctan((1/2)*x*2^(1/2))-873.2608343*(x^2+2)^6*2^(1/2)-374220*(0.29711673e-2*x^10+0.33673228e-1*x^8+.16607181*x^6+1.29737822*x^4+2.21433691*x^2+2.107985348)*x)*x/(x^2+2)^7+(1/16)*(3.484800000*2^(1/2)*(x^2+2)*arctan((1/2)*x*2^(1/2))*x+.8712000000*(x^2+2)^2/((1/2)*x^2+1)+5.473911040*(x^2+2)*2^(1/2)*x+5.227200000*x^2-22.99200001)/(x^2+2)^2-(1/4)*(.8712000000*2^(1/2)*(x^2+2)^2*arctan((1/2)*x*2^(1/2))+1.368477760*2^(1/2)*(x^2+2)^2-36*x*(-0.484000000e-1*x^2+.638666667))*x/(x^2+2)^3

(1)

ip := solve(f = 0, x)

.6540411301, 3126.002498+5414.398621*I, .4137989369+1.038962897*I, .6364817315+1.870977651*I, -.6364817315+1.870977651*I, -.4137989369+1.038962897*I, -.6540411301, -6252.010299, -.4137989369-1.038962897*I, -.6364817315-1.870977651*I, .6364817315-1.870977651*I, .4137989369-1.038962897*I, 3126.002498-5414.398621*I

(2)

cp := fsolve(numer(f) = 0, x)

.6540411302

(3)

``


 

Download help_fsolve_real_root.mw

Hi there.

As we all know if we multiply two polynomials f(x) and g(x) of degrees m and n respectively we get polynomial h(x)= f(x)*g(x) of degree m+n and with m+n+1 coefficients in general. Function modp1(('Multiply')(...)) doing this very well. But sometimes we don't need full resulting h(x) - just subset of monomials and subset of coefficients of h(x) - so we don't need to calculate all m+n+1 coefficients of h(x) and waste time and resources for that.

I would request some additional rework of modp1 package: by adding to modp1(('Multiply')(...)) two optional parameters - degrees of first and last calculating coefficients of h(x).

For example:

h:=modp1(Multiply(f, g,n-1,n+1), p) could calculate only monomials with n-1, n and n+1 degrees and set other monomials to zero.

Or maybe it should be new function:

h:=modp1(Multiply_Truncate(f, g,n-1,n+1), p)

 

Is it possible?

It would be great and very efficient in many tasks.

Thank you.


 

restart;

M__h := 0.352e-1;

0.352e-1

 

0.34e-1

 

0.8354e-1

 

0.96e-2

 

.123

 

0.7258e-1

 

0.214e-1

 

0.219e-1

 

.123

 

.7902

 

.11

 

0.136e-3

 

0.5e-1

 

0.8910e-1

 

0.45e-1

 

.7

 

.7214

 

1.354

 

0.235e-1

(1)

pdes := [diff(B(t, x), t) = M__h-beta__1*B(t, x)*G(t, x)/N__h+beta__2*B(t, x)*G(t, x)/N__h-mu__h*B(t, x)+sigma__h*E(t, x)*(diff(B(t, x), x, x)), diff(C(t, x), t) = beta__1*B(t, x)*G(t, x)/N__h-u[1]*C(t, x)/(1+C(t, x))-mu__h*C(t, x)*(diff(C(t, x), x, x)), diff(DD(t, x), t) = beta__2*DD(t, x)*G(t, x)/N__h-u[1]*DD(t, x)/(1+DD(t, x))-mu__h*DD(t, x)-delta__1*DD(t, x)*(diff(DD(t, x), x, x)), diff(E(t, x), t) = u[1]*C(t, x)/(1+C(t, x))+u[1]*DD(t, x)/(1+DD(t, x))-(mu__h+sigma__h)*E(t, x)*(diff(E(t, x), x, x)), diff(F(t, x), t) = M__b-beta__3*F(t, x)*C(t, x)/N__b+beta__4*F(t, x)*DD(t, x)/N__b-mu__b*F(t, x)*(diff(F(t, x), x, x)), diff(G(t, x), t) = beta__3*F(t, x)*C(t, x)/N__b+beta__4*F(t, x)*DD(t, x)/N__b-mu__b*G(t, x)*(diff(G(t, x), x, x))];

[diff(B(t, x), t) = 0.352e-1-0.891056911e-1*B(t, x)*G(t, x)-0.96e-2*B(t, x)+0.8910e-1*E(t, x)*(diff(diff(B(t, x), x), x)), diff(C(t, x), t) = .6791869919*B(t, x)*G(t, x)-0.45e-1*C(t, x)/(1+C(t, x))-0.96e-2*C(t, x)*(diff(diff(C(t, x), x), x)), diff(DD(t, x), t) = .5900813008*DD(t, x)*G(t, x)-0.45e-1*DD(t, x)/(1+DD(t, x))-0.96e-2*DD(t, x)-0.235e-1*DD(t, x)*(diff(diff(DD(t, x), x), x)), diff(E(t, x), t) = 0.45e-1*C(t, x)/(1+C(t, x))+0.45e-1*DD(t, x)/(1+DD(t, x))-0.9870e-1*E(t, x)*(diff(diff(E(t, x), x), x)), diff(F(t, x), t) = .7214-.1739837398*F(t, x)*C(t, x)+.1780487805*F(t, x)*DD(t, x)-1.354*F(t, x)*(diff(diff(F(t, x), x), x)), diff(G(t, x), t) = .1739837398*F(t, x)*C(t, x)+.1780487805*F(t, x)*DD(t, x)-1.354*G(t, x)*(diff(diff(G(t, x), x), x))]

(2)

bcs := [(D[2](B))(t, 0) = 0, (D[2](B))(t, 1) = 0, (D[2](C))(t, 0) = 0, (D[2](C))(t, 1) = 0, (D[2](DD))(t, 0) = 0, (D[2](DD))(t, 1) = 0, (D[2](E))(t, 0) = 0, (D[2](E))(t, 1) = 0, (D[2](F))(t, 0) = 0, (D[2](F))(t, 1) = 0, (D[2](G))(t, 0) = 0, (D[2](G))(t, 1) = 0, B(0, x) = 100, C(0, x) = 70, DD(0, x) = 50, E(0, x) = 70, F(0, x) = 100, G(0, x) = 70]

[(D[2](B))(t, 0) = 0, (D[2](B))(t, 1) = 0, (D[2](C))(t, 0) = 0, (D[2](C))(t, 1) = 0, (D[2](DD))(t, 0) = 0, (D[2](DD))(t, 1) = 0, (D[2](E))(t, 0) = 0, (D[2](E))(t, 1) = 0, (D[2](F))(t, 0) = 0, (D[2](F))(t, 1) = 0, (D[2](G))(t, 0) = 0, (D[2](G))(t, 1) = 0, B(0, x) = .100, C(0, x) = .70, DD(0, x) = .50, E(0, x) = .70, F(0, x) = .100, G(0, x) = .70]

(3)

sol := pdsolve(pdes, bcs, numeric);

module () local INFO; export plot, plot3d, animate, value, settings; option `Copyright (c) 2001 by Waterloo Maple Inc. All rights reserved.`; end module

(4)

sol:-plot3d([B(t, x), C(t, x)], t = 0 .. 20, x = 0 .. 20)

Error, (in pdsolve/numeric/plot3d) unable to compute solution for t>HFloat(0.25):
Newton iteration is not converging

 

``


 

Download spatial_1.mw

First 17 18 19 20 21 22 23 Last Page 19 of 1757