MaplePrimes Questions

I am trying to solve the folowing ODE with initial conditions t0=0,v0=0 and tf =80 with step 0.01 but the matrix that appears is not having the values!please help

f := proc (t, V) options operator, arrow; -9.81+0.563e-3*(0.1832e-2*abs(V)+0.51702e-1*(abs(V)^(3/2))+.4*V^2) end proc

V0 := 0:

t0 := 0:

tf := 80:

n := 1000

h := evalf((tf-t0)/n):

t := t0:

V := V0:

``

rk := proc (x0, tf, V0, n) local t, V, h, i, k1, k2, k3, k4, k, R; t := t0; V := V0; R := t, V; h := evalf((tf-t0)/n); for i to n do k1 := f(t, V); k2 := f(t+(1/2)*h, V+(1/2)*h*k1); k3 := f(t+(1/2)*h, V+(1/2)*h*k2); k4 := f(t+h, V+h*k3); k := (1/6)*k1+(1/3)*k2+(1/3)*k3+(1/6)*k4; V := V+h*k; t := t+h; R := R, t, V end do end proc:

R := rk(0, 80, 0, 300):

0, 0

 

.2666666667, -2.615998110

 

.5333333334, -5.231986045

 

.8000000001, -7.847954569

 

1.066666667, -10.46389463

 

1.333333334, -13.07979725

 

1.600000001, -15.69565351

 

1.866666668, -18.31145452

 

2.133333335, -20.92719143

 

2.400000002, -23.54285541

 

2.666666669, -26.15843764

 

2.933333336, -28.77392932

 

3.200000003, -31.38932166

 

3.466666670, -34.00460588

 

3.733333337, -36.61977321

 

4.000000004, -39.23481490

 

4.266666671, -41.84972219

 

4.533333338, -44.46448634

 

4.800000005, -47.07909861

 

5.066666672, -49.69355027

 

5.333333339, -52.30783259

 

5.600000006, -54.92193685

 

5.866666673, -57.53585433

 

6.133333340, -60.14957632

 

6.400000007, -62.76309410

 

6.666666674, -65.37639897

 

6.933333341, -67.98948223

 

7.200000008, -70.60233518

 

7.466666675, -73.21494912

 

7.733333342, -75.82731536

 

8.000000009, -78.43942520

 

8.266666676, -81.05126996

 

8.533333343, -83.66284095

 

8.800000010, -86.27412949

 

9.066666677, -88.88512689

 

9.333333344, -91.49582448

 

9.600000011, -94.10621358

 

9.866666678, -96.71628551

 

10.13333334, -99.32603160

 

10.40000001, -101.9354432

 

10.66666668, -104.5445116

 

10.93333335, -107.1532281

 

11.20000002, -109.7615841

 

11.46666669, -112.3695709

 

11.73333336, -114.9771799

 

12.00000003, -117.5844024

 

12.26666670, -120.1912297

 

12.53333337, -122.7976532

 

12.80000004, -125.4036642

 

13.06666671, -128.0092541

 

13.33333338, -130.6144142

 

13.60000005, -133.2191358

 

13.86666672, -135.8234103

 

14.13333339, -138.4272290

 

14.40000006, -141.0305833

 

14.66666673, -143.6334645

 

14.93333340, -146.2358640

 

15.20000007, -148.8377732

 

15.46666674, -151.4391834

 

15.73333341, -154.0400859

 

16.26666675, -159.2403334

 

16.53333342, -161.8396611

 

16.80000009, -164.4384465

 

17.06666676, -167.0366810

 

17.33333343, -169.6343560

 

17.60000010, -172.2314629

 

17.86666677, -174.8279930

 

18.13333344, -177.4239376

 

18.40000011, -180.0192882

 

18.66666678, -182.6140360

 

18.93333345, -185.2081725

 

19.20000012, -187.8016890

 

19.46666679, -190.3945769

 

19.73333346, -192.9868275

 

20.00000013, -195.5784322

 

20.26666680, -198.1693824

 

20.53333347, -200.7596694

 

20.80000014, -203.3492846

 

21.06666681, -205.9382194

 

21.33333348, -208.5264652

 

21.60000015, -211.1140133

 

21.86666682, -213.7008550

 

22.13333349, -216.2869818

 

22.40000016, -218.8723850

 

22.66666683, -221.4570560

 

22.93333350, -224.0409862

 

23.20000017, -226.6241669

 

23.46666684, -229.2065895

 

23.73333351, -231.7882454

 

24.00000018, -234.3691260

 

24.26666685, -236.9492226

 

24.53333352, -239.5285266

 

24.80000019, -242.1070294

 

25.06666686, -244.6847224

 

25.33333353, -247.2615969

 

25.60000020, -249.8376443

 

25.86666687, -252.4128560

 

26.13333354, -254.9872234

 

26.40000021, -257.5607378

 

26.66666688, -260.1333906

 

26.93333355, -262.7051732

 

27.20000022, -265.2760770

 

27.46666689, -267.8460934

 

27.73333356, -270.4152137

 

28.00000023, -272.9834293

 

28.26666690, -275.5507316

 

28.53333357, -278.1171120

 

28.80000024, -280.6825619

 

29.06666691, -283.2470726

 

29.33333358, -285.8106356

 

29.60000025, -288.3732422

 

29.86666692, -290.9348838

 

30.13333359, -293.4955517

 

30.40000026, -296.0552374

 

30.66666693, -298.6139322

 

30.93333360, -301.1716276

 

31.20000027, -303.7283149

 

31.46666694, -306.2839855

 

31.73333361, -308.8386307

 

32.00000028, -311.3922420

 

32.26666695, -313.9448108

 

32.53333362, -316.4963284

 

32.80000029, -319.0467862

 

33.06666696, -321.5961756

 

33.33333363, -324.1444880

 

33.60000030, -326.6917148

 

33.86666697, -329.2378474

 

34.13333364, -331.7828772

 

34.40000031, -334.3267955

 

34.66666698, -336.8695937

 

34.93333365, -339.4112633

 

35.20000032, -341.9517956

 

35.46666699, -344.4911820

 

35.73333366, -347.0294139

 

36.00000033, -349.5664827

 

36.26666700, -352.1023797

 

36.53333367, -354.6370964

 

36.80000034, -357.1706242

 

37.06666701, -359.7029544

 

37.33333368, -362.2340784

 

37.60000035, -364.7639877

 

37.86666702, -367.2926736

 

38.13333369, -369.8201275

 

38.40000036, -372.3463408

 

38.66666703, -374.8713049

 

38.93333370, -377.3950112

 

39.20000037, -379.9174511

 

39.46666704, -382.4386159

 

39.73333371, -384.9584971

 

40.00000038, -387.4770861

 

40.26666705, -389.9943742

 

40.53333372, -392.5103529

 

40.80000039, -395.0250135

 

41.06666706, -397.5383474

 

41.33333373, -400.0503461

 

41.60000040, -402.5610009

 

41.86666707, -405.0703032

 

42.13333374, -407.5782444

 

42.40000041, -410.0848159

 

42.66666708, -412.5900091

 

42.93333375, -415.0938154

 

43.20000042, -417.5962262

 

43.46666709, -420.0972329

 

43.73333376, -422.5968269

 

44.00000043, -425.0949996

 

44.26666710, -427.5917424

 

44.53333377, -430.0870466

 

44.80000044, -432.5809037

 

45.06666711, -435.0733051

 

45.33333378, -437.5642421

 

45.60000045, -440.0537062

 

45.86666712, -442.5416888

 

46.13333379, -445.0281812

 

46.40000046, -447.5131749

 

46.66666713, -449.9966612

 

46.93333380, -452.4786316

 

47.20000047, -454.9590775

 

47.46666714, -457.4379902

 

47.73333381, -459.9153612

 

48.00000048, -462.3911818

 

48.26666715, -464.8654435

 

48.53333382, -467.3381376

 

48.80000049, -469.8092556

 

49.06666716, -472.2787888

 

49.33333383, -474.7467287

 

49.60000050, -477.2130666

 

49.86666717, -479.6777940

 

50.13333384, -482.1409022

 

50.40000051, -484.6023827

 

50.66666718, -487.0622268

 

50.93333385, -489.5204260

 

51.20000052, -491.9769716

 

51.46666719, -494.4318551

 

51.73333386, -496.8850678

 

52.00000053, -499.3366012

 

52.26666720, -501.7864466

 

52.53333387, -504.2345955

 

52.80000054, -506.6810393

 

53.06666721, -509.1257693

 

53.33333388, -511.5687770

 

53.60000055, -514.0100537

 

53.86666722, -516.4495909

 

54.13333389, -518.8873800

 

54.40000056, -521.3234124

 

54.66666723, -523.7576794

 

54.93333390, -526.1901725

 

55.20000057, -528.6208831

 

55.46666724, -531.0498026

 

55.73333391, -533.4769224

 

56.00000058, -535.9022339

 

56.26666725, -538.3257285

 

56.53333392, -540.7473976

 

56.80000059, -543.1672326

 

57.06666726, -545.5852249

 

57.33333393, -548.0013659

 

57.60000060, -550.4156470

 

57.86666727, -552.8280597

 

58.13333394, -555.2385953

 

58.40000061, -557.6472452

 

58.66666728, -560.0540009

 

58.93333395, -562.4588537

 

59.20000062, -564.8617951

 

59.46666729, -567.2628164

 

59.73333396, -569.6619091

 

60.00000063, -572.0590645

 

60.26666730, -574.4542741

 

60.53333397, -576.8475293

 

60.80000064, -579.2388215

 

61.06666731, -581.6281420

 

61.33333398, -584.0154823

 

61.60000065, -586.4008338

 

61.86666732, -588.7841879

 

62.13333399, -591.1655361

 

62.40000066, -593.5448697

 

62.66666733, -595.9221801

 

62.93333400, -598.2974587

 

63.20000067, -600.6706970

 

63.46666734, -603.0418863

 

63.73333401, -605.4110181

 

64.00000068, -607.7780837

 

64.26666735, -610.1430746

 

64.53333402, -612.5059822

 

64.80000069, -614.8667979

 

65.06666736, -617.2255130

 

65.33333403, -619.5821190

 

65.60000070, -621.9366073

 

65.86666737, -624.2889694

 

66.13333404, -626.6391966

 

66.40000071, -628.9872803

 

66.66666738, -631.3332119

 

66.93333405, -633.6769829

 

67.20000072, -636.0185846

 

67.46666739, -638.3580085

 

67.73333406, -640.6952460

 

68.00000073, -643.0302885

 

68.26666740, -645.3631274

 

68.53333407, -647.6937540

 

68.80000074, -650.0221599

 

69.06666741, -652.3483364

 

69.33333408, -654.6722749

 

69.60000075, -656.9939668

 

69.86666742, -659.3134036

 

70.13333409, -661.6305766

 

70.40000076, -663.9454773

 

70.66666743, -666.2580971

 

70.93333410, -668.5684273

 

71.20000077, -670.8764594

 

71.46666744, -673.1821848

 

71.73333411, -675.4855950

 

72.00000078, -677.7866813

 

72.26666745, -680.0854351

 

72.53333412, -682.3818479

 

72.80000079, -684.6759110

 

73.06666746, -686.9676159

 

73.33333413, -689.2569539

 

73.60000080, -691.5439165

 

73.86666747, -693.8284951

 

74.13333414, -696.1106811

 

74.40000081, -698.3904660

 

74.66666748, -700.6678411

 

74.93333415, -702.9427978

 

75.20000082, -705.2153276

 

75.46666749, -707.4854218

 

75.73333416, -709.7530719

 

76.00000083, -712.0182693

 

76.26666750, -714.2810054

 

76.53333417, -716.5412716

 

76.80000084, -718.7990593

 

77.06666751, -721.0543600

 

77.33333418, -723.3071650

 

77.60000085, -725.5574658

 

77.86666752, -727.8052537

 

78.13333419, -730.0505202

 

78.40000086, -732.2932567

 

78.66666753, -734.5334546

 

78.93333420, -736.7711053

 

79.73333421, -743.4686884

 

``

 

 

Download runge_kutta4.mw

What's the best way to get the eigenvector associated with a certain eigenvalue?

Specifically, given the nature of my matrix A, I know that there ALWAYS exists an eigenvector with eigenvalue 1.  Is there a quick way to extract this without looping through the output of Eigenvectors() and checking each one?

Hello,

 

I'm modeling the simple DC motor system in Maple.
The equations describing the system;

eq1:=J*diff(theta(t),t,t)+b*diff(theta(t),t)=K*i(t):
eq2:=L*diff(i(t),t)+R*i(t)=V(t)-K*diff(theta(t),t):
DCMotor:=[eq1,eq2];

First, I create the system using DiffEquation:

Sys:=DiffEquation(DCMotor,[V(t)],[theta(t)]);

And now I have problem. The input var is V(t) (input voltage) and the output var is theta(t) (position of the rotor).

But I wont to have in output var not position of the rotor but speed of the fotor - diff(theta(t),t)

How to set output var for diff(theta(t),t) (the speed of the motor)?

 

Best

Rariusz

 

Hi,

I'd like to create a graded drag and drop exercise so that my students can click on an item and place it into various places (one of which is the correct place). I don't think any of the Question types support this though .. is there a way to build such a quesiton in Maple TA ?

Thank you for your help!

Elisabeth

Dear All

I want to use subs command to replace products in expression with entries from matrix. When I use "subs" command it returns same result as original without any substitution. But when I use "algsubs", it works fine, but is complicated to apply if there are large number of substitutions.

Please see followings:

 

M := Matrix([[U[1], U[2], -epsilon*U[1]+U[3], U[4], U[5], U[6]], [U[1], U[2], -epsilon*U[2]+U[3], epsilon*U[1]+U[4], epsilon*U[2]+U[5], U[6]], [e^epsilon*U[1], e^epsilon*U[2], U[3], U[4], U[5], U[6]], [U[1], -epsilon*U[1]+U[2], U[3], U[4], -epsilon*U[4]+U[5], U[6]-epsilon*(U[3]+2*U[5])+epsilon^2*U[4], U[1], e^(-epsilon)*U[2], U[3], e^epsilon*U[4], U[5], e^(-epsilon)*U[6]], [U[1], e^(-epsilon)*U[2], U[3], e^epsilon*U[4], U[5], e^(-epsilon)*U[6]], [epsilon*U[2]+U[1], U[2], U[3], U[4]+epsilon*(U[3]+2*U[5])+epsilon^2*U[6], epsilon*U[6]+U[5], U[6]]])

M := Vector(4, {(1) = ` 6 x 12 `*Matrix, (2) = `Data Type: `*anything, (3) = `Storage: `*rectangular, (4) = `Order: `*Fortran_order})

(1)

M[5, 1]; 1; M[5, 3]

U[1]

 

U[3]

(2)

expand((w*T[5]+z*T[4]+T[6])*(U[1]*a[1]+U[2]*a[2]+U[3]*a[3]+U[4]*a[4]+U[5]*a[5]+U[6]*a[6]))

w*T[5]*U[1]*a[1]+w*T[5]*U[2]*a[2]+w*T[5]*U[3]*a[3]+w*T[5]*U[4]*a[4]+w*T[5]*U[5]*a[5]+w*T[5]*U[6]*a[6]+z*T[4]*U[1]*a[1]+z*T[4]*U[2]*a[2]+z*T[4]*U[3]*a[3]+z*T[4]*U[4]*a[4]+z*T[4]*U[5]*a[5]+z*T[4]*U[6]*a[6]+T[6]*U[1]*a[1]+T[6]*U[2]*a[2]+T[6]*U[3]*a[3]+T[6]*U[4]*a[4]+T[6]*U[5]*a[5]+T[6]*U[6]*a[6]

(3)

subs(T[5]*U[1] = M[5, 1], T[5]*U[2] = M[5, 2], T[5]*U[3] = M[5, 3], w*T[5]*U[1]*a[1]+w*T[5]*U[2]*a[2]+w*T[5]*U[3]*a[3]+w*T[5]*U[4]*a[4]+w*T[5]*U[5]*a[5]+w*T[5]*U[6]*a[6]+z*T[4]*U[1]*a[1]+z*T[4]*U[2]*a[2]+z*T[4]*U[3]*a[3]+z*T[4]*U[4]*a[4]+z*T[4]*U[5]*a[5]+z*T[4]*U[6]*a[6]+T[6]*U[1]*a[1]+T[6]*U[2]*a[2]+T[6]*U[3]*a[3]+T[6]*U[4]*a[4]+T[6]*U[5]*a[5]+T[6]*U[6]*a[6])

w*T[5]*U[1]*a[1]+w*T[5]*U[2]*a[2]+w*T[5]*U[3]*a[3]+w*T[5]*U[4]*a[4]+w*T[5]*U[5]*a[5]+w*T[5]*U[6]*a[6]+z*T[4]*U[1]*a[1]+z*T[4]*U[2]*a[2]+z*T[4]*U[3]*a[3]+z*T[4]*U[4]*a[4]+z*T[4]*U[5]*a[5]+z*T[4]*U[6]*a[6]+T[6]*U[1]*a[1]+T[6]*U[2]*a[2]+T[6]*U[3]*a[3]+T[6]*U[4]*a[4]+T[6]*U[5]*a[5]+T[6]*U[6]*a[6]

(4)

Why this command subs in not replacing  products  T[5]*U[1], T[5]*U[2], T[5]*U[3] from matrix (1)NULL?????

algsubs(T[5]*U[1] = M[5, 1], w*T[5]*U[1]*a[1]+w*T[5]*U[2]*a[2]+w*T[5]*U[3]*a[3]+w*T[5]*U[4]*a[4]+w*T[5]*U[5]*a[5]+w*T[5]*U[6]*a[6]+z*T[4]*U[1]*a[1]+z*T[4]*U[2]*a[2]+z*T[4]*U[3]*a[3]+z*T[4]*U[4]*a[4]+z*T[4]*U[5]*a[5]+z*T[4]*U[6]*a[6]+T[6]*U[1]*a[1]+T[6]*U[2]*a[2]+T[6]*U[3]*a[3]+T[6]*U[4]*a[4]+T[6]*U[5]*a[5]+T[6]*U[6]*a[6])

z*T[4]*U[2]*a[2]+z*T[4]*U[3]*a[3]+z*T[4]*U[4]*a[4]+z*T[4]*U[5]*a[5]+z*T[4]*U[6]*a[6]+T[6]*U[2]*a[2]+T[6]*U[3]*a[3]+T[6]*U[4]*a[4]+T[6]*U[5]*a[5]+T[6]*U[6]*a[6]+(w*U[2]*a[2]+w*U[3]*a[3]+w*U[4]*a[4]+w*U[5]*a[5]+w*U[6]*a[6])*T[5]+(z*T[4]*a[1]+w*a[1]+T[6]*a[1])*U[1]

(5)

It worked !!!! But If I try "algsubs" like :

algsubs(T[5]*U[1] = M[5, 1], T[5]*U[2] = M[5, 2], T[5]*U[3] = M[5, 3], w*T[5]*U[1]*a[1]+w*T[5]*U[2]*a[2]+w*T[5]*U[3]*a[3]+w*T[5]*U[4]*a[4]+w*T[5]*U[5]*a[5]+w*T[5]*U[6]*a[6]+z*T[4]*U[1]*a[1]+z*T[4]*U[2]*a[2]+z*T[4]*U[3]*a[3]+z*T[4]*U[4]*a[4]+z*T[4]*U[5]*a[5]+z*T[4]*U[6]*a[6]+T[6]*U[1]*a[1]+T[6]*U[2]*a[2]+T[6]*U[3]*a[3]+T[6]*U[4]*a[4]+T[6]*U[5]*a[5]+T[6]*U[6]*a[6])

Error, invalid input: algsubs expects its 3rd argument, x, to be of type {name, function, list({name, function}), set({name, function})}, but received T[5]*U[3] = U[3]

 

See it is giving error !!! but still I can use "algsubs" as follow:

 

expand(algsubs(T[6]*U[6] = M[6, 6], expand(algsubs(T[6]*U[5] = M[6, 5], expand(algsubs(T[6]*U[4] = M[6, 4], expand(algsubs(T[6]*U[3] = M[6, 3], expand(algsubs(T[6]*U[2] = M[6, 2], expand(algsubs(T[6]*U[1] = M[6, 1], expand(algsubs(T[5]*U[6] = M[5, 6], expand(algsubs(T[5]*U[5] = M[5, 5], expand(algsubs(T[5]*U[4] = M[5, 4], expand(algsubs(T[5]*U[3] = M[5, 3], expand(algsubs(T[5]*U[2] = M[5, 2], expand(algsubs(T[5]*U[1] = M[5, 1], expand(algsubs(T[4]*U[6] = M[4, 6], expand(algsubs(T[4]*U[5] = M[4, 5], expand(algsubs(T[4]*U[4] = M[4, 4], expand(algsubs(T[4]*U[3] = M[4, 3], expand(algsubs(T[4]*U[2] = M[4, 2], expand(algsubs(T[4]*U[1] = M[4, 1], expand(algsubs(T[3]*U[6] = M[3, 6], expand(algsubs(T[3]*U[5] = M[3, 5], expand(algsubs(T[3]*U[4] = M[3, 4], expand(algsubs(T[3]*U[3] = M[3, 3], expand(algsubs(T[3]*U[2] = M[3, 2], expand(algsubs(T[3]*U[1] = M[3, 1], expand(algsubs(T[2]*U[6] = M[2, 6], expand(algsubs(T[2]*U[5] = M[2, 5], expand(algsubs(T[2]*U[4] = M[2, 4], expand(algsubs(T[2]*U[3] = M[2, 3], expand(algsubs(T[2]*U[3] = M[2, 3], expand(algsubs(T[2]*U[2] = M[2, 2], expand(algsubs(T[2]*U[1] = M[2, 1], expand(algsubs(T[1]*U[6] = M[1, 6], expand(algsubs(T[1]*U[5] = M[1, 5], expand(algsubs(T[1]*U[4] = M[1, 4], expand(algsubs(T[1]*U[3] = M[1, 3], expand(algsubs(T[1]*U[2] = M[1, 2], expand(algsubs(T[1]*U[1] = M[1, 1], w*T[5]*U[1]*a[1]+w*T[5]*U[2]*a[2]+w*T[5]*U[3]*a[3]+w*T[5]*U[4]*a[4]+w*T[5]*U[5]*a[5]+w*T[5]*U[6]*a[6]+z*T[4]*U[1]*a[1]+z*T[4]*U[2]*a[2]+z*T[4]*U[3]*a[3]+z*T[4]*U[4]*a[4]+z*T[4]*U[5]*a[5]+z*T[4]*U[6]*a[6]+T[6]*U[1]*a[1]+T[6]*U[2]*a[2]+T[6]*U[3]*a[3]+T[6]*U[4]*a[4]+T[6]*U[5]*a[5]+T[6]*U[6]*a[6]))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))

a[5]*epsilon*U[6]+a[1]*epsilon*U[2]+U[6]*a[6]+U[2]*a[2]+U[3]*a[3]+U[4]*a[4]+U[5]*a[5]+z*a[1]*U[1]+epsilon^2*U[6]*a[4]+epsilon*U[3]*a[4]+2*epsilon*U[5]*a[4]+U[6]*w*a[6]/e^epsilon+U[4]*e^epsilon*w*a[4]+w*a[1]*U[1]+U[2]*w*a[2]/e^epsilon+epsilon^2*z*U[4]*a[6]-epsilon*z*U[3]*a[6]-2*epsilon*z*U[5]*a[6]+w*U[3]*a[3]+w*U[5]*a[5]-z*a[5]*epsilon*U[4]-z*a[2]*epsilon*U[1]+z*U[2]*a[2]+z*U[3]*a[3]+z*U[4]*a[4]+z*U[5]*a[5]+z*U[6]*a[6]+a[1]*U[1]

(6)

But this procedure is complex, I want to use "subs" only. I know it complexity of terms in (3) might me creating problem.

collect(a[5]*epsilon*U[6]+a[1]*epsilon*U[2]+U[6]*a[6]+U[2]*a[2]+U[3]*a[3]+U[4]*a[4]+U[5]*a[5]+z*a[1]*U[1]+epsilon^2*U[6]*a[4]+epsilon*U[3]*a[4]+2*epsilon*U[5]*a[4]+U[6]*w*a[6]/e^epsilon+U[4]*e^epsilon*w*a[4]+w*a[1]*U[1]+U[2]*w*a[2]/e^epsilon+epsilon^2*z*U[4]*a[6]-epsilon*z*U[3]*a[6]-2*epsilon*z*U[5]*a[6]+w*U[3]*a[3]+w*U[5]*a[5]-z*a[5]*epsilon*U[4]-z*a[2]*epsilon*U[1]+z*U[2]*a[2]+z*U[3]*a[3]+z*U[4]*a[4]+z*U[5]*a[5]+z*U[6]*a[6]+a[1]*U[1], [U[1], U[2], U[3], U[4], U[5], U[6]])

(-epsilon*z*a[2]+w*a[1]+z*a[1]+a[1])*U[1]+(a[1]*epsilon+a[2]+w*a[2]/e^epsilon+z*a[2])*U[2]+(-epsilon*z*a[6]+epsilon*a[4]+w*a[3]+z*a[3]+a[3])*U[3]+(epsilon^2*z*a[6]+e^epsilon*w*a[4]-z*a[5]*epsilon+z*a[4]+a[4])*U[4]+(-2*epsilon*z*a[6]+2*epsilon*a[4]+w*a[5]+z*a[5]+a[5])*U[5]+(a[5]*epsilon+a[6]+epsilon^2*a[4]+w*a[6]/e^epsilon+z*a[6])*U[6]

(7)

``

 

Download Subs_not_giving_result.mwSubs_not_giving_result.mw

Regards

plot3d(x^2+y^2, x = -1 .. 1, y = -1 .. 1);

i searched this 

http://www.maplesoft.com/support/helpJP/Maple/view.aspx?path=NAG/d06abc

and

DelaunayTriangulation(points) in help file

is there any example to show how to use these?

 

Hi all.

1. When plotting the function y=(x-1)/(x-2) on Maple 2016, the asymptote along with the actual curve is drawn as "solid lines/curves".  The asymptote should have been drawn as a dotted line to indicate that x=2 is an asymptote so as to distinguish itself from the actual curve.   This confuses my students.

2. The other question is if I modify the equation to y=1+ (x-1)/(x-2), the plot only shows the asymptote at x=2 (again, a solid line).  There should be another asymptote at y=1, which is not shown.

Can someone please advise if Maple would fix this in the next update?    If not, how should one "fix" this when graphing such functions?

Many thanks.

Hello, I have the system of equations in many vars as below, I want to make an implicit plot in Maple with the projection on 3 vars, for example, in this case (x,y,t1). The range is x[-10,10], y[-10,10], t1[-Pi,Pi] and the rest of the vars are [-Pi,Pi]. Does anyone know how to do it? We have also the inequalities in the system.

 

f1:= cos(t1)+1.35*cos(p1)-x;

f2:= sin(t1)+1.35*sin(p1)-y;

f3:= cos(t2)+1.35*cos(p2)-x +1.15;

f4:= sin(t2)+1.35*sin(p2)-y;

f5:= cos(t3)+1.35*cos(p3)-x+0.575;

f6:= sin(t3)+1.35*sin(p3)-y +0.995;

f7:= cos(t1)*cos(t2)*cos(t3)*sin(p1)*sin(p2)*sin(p3)-cos(t1)*cos(t2)*sin(t3)*sin(p1)*sin(p2)*cos(p3)-cos(t1)*sin(t2)*cos(t3)*sin(p1)*cos(p2)*sin(p3)+cos(t1)*sin(t2)*sin(t3)*sin(p1)*cos(p2)*cos(p3)-sin(t1)*cos(t2)*cos(t3)*cos(p1)*sin(p2)*sin(p3)+sin(t1)*cos(t2)*sin(t3)*cos(p1)*sin(p2)*cos(p3)+sin(t1)*sin(t2)*cos(t3)*cos(p1)*cos(p2)*sin(p3)-sin(t1)*sin(t2)*sin(t3)*cos(p1)*cos(p2)*cos(p3);

f8:= cos(t1)*sin(p1)-sin(t1)*cos(p1) >= 0;

f9:= cos(t2)*sin(p2)-sin(t2)*cos(p2) >= 0;

f10:= cos(t3)*sin(p3)-sin(t3)*cos(p3) >= 0;

Consider you have an expression like:

f := (125*x^3+525*x^2*y+735*x*y^2+343*y^3-36*z^2)/(-216*z^3+25*x^2+70*x*y+49*y^2)

By hand and spending much time maybe we can find out so far that:

f := ((5*x+7*y)^3-(6*z)^2)/((5*x+7*y)^2-(6*z)^3)

I thought that the code `combine(f)` would give us the second one, but it didn’t. Is there anything obvious I cannot see?

Thanks for the time.

I am fairly new to Maple and cannot figure this out.

I'm having some trouble in computing formally the derivatives of
some composite functions

h:= x -> f(g(x)),

imposing conditions on the unknown functions f and g.

For example, I cannot find a way to impose g'(0)=0,
such that the computation of

h'(0)

simply gives 0.

 

Thank you in advance!

 

restart

with(LinearAlgebra):

with(StringTools):

FormatTime("%m-%d-%Y, %H:%M")

NULL

Pile spacing of each pile from a refrence pile:

S__p := Matrix(9, 9, {(1, 1) = 0., (1, 2) = 3.30, (1, 3) = 6.60, (1, 4) = 3.30, (1, 5) = 4.67, (1, 6) = 7.38, (1, 7) = 6.60, (1, 8) = 7.38, (1, 9) = 9.33, (2, 1) = 3.30, (2, 2) = 0., (2, 3) = 3.30, (2, 4) = 4.67, (2, 5) = 3.30, (2, 6) = 4.67, (2, 7) = 7.38, (2, 8) = 6.60, (2, 9) = 7.38, (3, 1) = 6.60, (3, 2) = 3.30, (3, 3) = 0., (3, 4) = 7.38, (3, 5) = 4.67, (3, 6) = 3.30, (3, 7) = 9.33, (3, 8) = 7.38, (3, 9) = 6.60, (4, 1) = 3.30, (4, 2) = 4.67, (4, 3) = 7.38, (4, 4) = 0., (4, 5) = 3.30, (4, 6) = 6.60, (4, 7) = 3.30, (4, 8) = 4.67, (4, 9) = 7.38, (5, 1) = 4.67, (5, 2) = 3.30, (5, 3) = 4.67, (5, 4) = 3.30, (5, 5) = 0., (5, 6) = 3.30, (5, 7) = 4.67, (5, 8) = 3.30, (5, 9) = 4.67, (6, 1) = 7.38, (6, 2) = 4.67, (6, 3) = 3.30, (6, 4) = 6.60, (6, 5) = 3.30, (6, 6) = 0., (6, 7) = 7.38, (6, 8) = 4.67, (6, 9) = 3.30, (7, 1) = 6.60, (7, 2) = 7.38, (7, 3) = 9.33, (7, 4) = 3.30, (7, 5) = 4.67, (7, 6) = 7.38, (7, 7) = 0., (7, 8) = 3.30, (7, 9) = 6.60, (8, 1) = 7.38, (8, 2) = 6.60, (8, 3) = 7.38, (8, 4) = 4.67, (8, 5) = 3.30, (8, 6) = 4.67, (8, 7) = 3.30, (8, 8) = 0., (8, 9) = 3.30, (9, 1) = 9.33, (9, 2) = 7.38, (9, 3) = 6.60, (9, 4) = 7.38, (9, 5) = 4.67, (9, 6) = 3.30, (9, 7) = 6.60, (9, 8) = 3.30, (9, 9) = 0.})

Angle of each pile from a reference pile in radian:

xi := Matrix(9, 9, {(1, 1) = 0., (1, 2) = 0., (1, 3) = 0., (1, 4) = 1.57, (1, 5) = .785, (1, 6) = .464, (1, 7) = 1.57, (1, 8) = 1.11, (1, 9) = .785, (2, 1) = 3.14, (2, 2) = 0., (2, 3) = 0., (2, 4) = 2.36, (2, 5) = 1.57, (2, 6) = .785, (2, 7) = 2.03, (2, 8) = 1.57, (2, 9) = 1.11, (3, 1) = 3.14, (3, 2) = 3.14, (3, 3) = 0., (3, 4) = 2.68, (3, 5) = 2.36, (3, 6) = 1.57, (3, 7) = 2.36, (3, 8) = 2.03, (3, 9) = 1.57, (4, 1) = 1.57, (4, 2) = .785, (4, 3) = .464, (4, 4) = 0., (4, 5) = 0., (4, 6) = 0., (4, 7) = 1.57, (4, 8) = .785, (4, 9) = .464, (5, 1) = 2.36, (5, 2) = 1.57, (5, 3) = .785, (5, 4) = 3.14, (5, 5) = 0., (5, 6) = 0., (5, 7) = 2.36, (5, 8) = 1.57, (5, 9) = .785, (6, 1) = 2.68, (6, 2) = 2.36, (6, 3) = 1.57, (6, 4) = 3.14, (6, 5) = 3.14, (6, 6) = 0., (6, 7) = 2.68, (6, 8) = 2.36, (6, 9) = 1.57, (7, 1) = 1.57, (7, 2) = 1.11, (7, 3) = .785, (7, 4) = 1.57, (7, 5) = .785, (7, 6) = .464, (7, 7) = 0., (7, 8) = 0., (7, 9) = 0., (8, 1) = 2.03, (8, 2) = 1.57, (8, 3) = 1.11, (8, 4) = 2.36, (8, 5) = 1.57, (8, 6) = .785, (8, 7) = 3.14, (8, 8) = 0., (8, 9) = 0., (9, 1) = 2.36, (9, 2) = 2.03, (9, 3) = 1.57, (9, 4) = 2.68, (9, 5) = 2.36, (9, 6) = 1.57, (9, 7) = 3.14, (9, 8) = 3.14, (9, 9) = 0.})

NULL

Angle of each pile from a reference pile in degrees:

theta := Matrix(9, 9, {(1, 1) = 0., (1, 2) = 0., (1, 3) = 0., (1, 4) = 90.0, (1, 5) = 44.8, (1, 6) = 26.6, (1, 7) = 90.0, (1, 8) = 63.6, (1, 9) = 44.8, (2, 1) = 180., (2, 2) = 0., (2, 3) = 0., (2, 4) = 135., (2, 5) = 90.0, (2, 6) = 44.8, (2, 7) = 116., (2, 8) = 90.0, (2, 9) = 63.6, (3, 1) = 180., (3, 2) = 180., (3, 3) = 0., (3, 4) = 153., (3, 5) = 135., (3, 6) = 90.0, (3, 7) = 135., (3, 8) = 116., (3, 9) = 90.0, (4, 1) = 90.0, (4, 2) = 44.8, (4, 3) = 26.6, (4, 4) = 0., (4, 5) = 0., (4, 6) = 0., (4, 7) = 90.0, (4, 8) = 44.8, (4, 9) = 26.6, (5, 1) = 135., (5, 2) = 90.0, (5, 3) = 44.8, (5, 4) = 180., (5, 5) = 0., (5, 6) = 0., (5, 7) = 135., (5, 8) = 90.0, (5, 9) = 44.8, (6, 1) = 153., (6, 2) = 135., (6, 3) = 90.0, (6, 4) = 180., (6, 5) = 180., (6, 6) = 0., (6, 7) = 153., (6, 8) = 135., (6, 9) = 90.0, (7, 1) = 90.0, (7, 2) = 63.6, (7, 3) = 44.8, (7, 4) = 90.0, (7, 5) = 44.8, (7, 6) = 26.6, (7, 7) = 0., (7, 8) = 0., (7, 9) = 0., (8, 1) = 116., (8, 2) = 90.0, (8, 3) = 63.6, (8, 4) = 135., (8, 5) = 90.0, (8, 6) = 44.8, (8, 7) = 180., (8, 8) = 0., (8, 9) = 0., (9, 1) = 135., (9, 2) = 116., (9, 3) = 90.0, (9, 4) = 153., (9, 5) = 135., (9, 6) = 90.0, (9, 7) = 180., (9, 8) = 180., (9, 9) = 0.})

The total number of piles in a pile group:

n__pile := 9

The length of each pile in the pile group in (meters):

L__p := 12.182   

 

``

Diameter of pile (meter):

`φ__pile` := .457

The radius of each pile (meter):

r__0 := evalf((1/2)*`φ__pile`, 2)

NULL

`ρ__A` := 1.0

NULL

D__s := 0.5e-1

NULL

Omega := 46``

``

`ν__soil` := .5

NULL

E__p := 3.2*10^10``

NULL

G__s := 5.7*10^7``

NULL

 

V__s := 182.88

NULL

``

NULL

NULL

V__La := 395.845``

NULL

`&alpha;__&nu;` := proc (S__p, xi, n__pile) local `&alpha;__tmp`, flist, i, j, `&alpha;__1`, `&alpha;__2`, `&alpha;__3`; global D__s, V__La, V__s, Omega; description "Calculate the horizontal interaction factor from a refernce pile in a pile group"; flist := []; for i to `n__pile ` do for j to n__pile do `&alpha;__1`[i, j] := evalf[4](.5*(S__p[i, j]/`&phi;__pile`)^.5*exp(-D__s*Omega*S__p[i, j]/V__La)*exp(-I*Omega*S__p[i, j]/V__La)); `&alpha;__2`[i, j] := evalf[4]((3/4)*(S__p[i, j]/`&phi;__pile`)^.5*exp(-D__s*Omega*S__p[i, j]/V__s)*exp(-I*Omega*S__p[i, j]/V__s)/sqrt(2)); `&alpha;__3`[i, j] := evalf[4](`&alpha;__1`[i, j]*cos(xi[i, j])^2+`&alpha;__2`[i, j]*sin(xi[i, j])^2); if i = j then `&alpha;__tmp`[i, i] := 1.0 elif xi[i, j] = 0. then `&alpha;__tmp`[i, j] := `&alpha;__1`[i, j] elif xi[i, j] = 1.57 then `&alpha;__tmp`[i, j] := `&alpha;__2`[i, j] elif xi[i, j] <> 1.57 and xi[i, j] <> 0. then `&alpha;__tmp`[i, j] := `&alpha;__3`[i, j] end if; flist := [op(flist), `&alpha;__tmp`[i, j]] end do end do end proc:

NULL

NULL

NULL

NULLNULL

NULLNULLNULL

`&alpha;__v` := eval(Matrix(n__pile, n__pile, `&alpha;__&nu;`(S__p, xi, n__pile)))``

NULL

``

``

``

NULL

NULL

``

NULL

 

Download Test_pile.mw

Good Evening Maple Community,

Can anybody please help in the attached file. Highlighted in yellow are the three equations I'm trying to program in Maple procedure, but could not get any output due to an error message. Please Help.

 

Boyer 

Is this what is supposed to happen when linking Maple to Matlab?  I was hoping it would open the actual Matlab program and write variables to the current worksheet.  I'm I missing something?

Nothing even pops up on the Matlab Command window. What good is this?

 

I want to define a proc in Maple which gets an integer n and gives a set of arrays like {1,2,3}^n, meaning {[a[1],...,a[n]]|for all 1<=i<=n a[i] in {1,2,3}}. For a fix n, obviously one can use n-for, for example when n is 2

A:=Array(1..3,1..3):
for a[1] from 1 by 1 to 3 do
 for a[2] from 1 by 1 to 3 do
  A(a[1],a[2]):=[a[1],a[2]];
 end for:
end for:

But when n is not fixed, how can I tell Maple to consider n for?! Also how can ask him to make a n-dimensional array?

Hi 

I have 3 equations

eq1 := vs = (Rs+Z)*i1+Z*i3

eq2 := A*vi = Z*i1+(Z+Rf+ro)*i3

eq3:= vo = (Rf+Z)*i3+Z*i1

 

and I want to solve for vo/vs . How to do that ?

the expected solution is 

Hi;

Could anyone help me to compute the lie derivative of the function h:R^3-->R with respect to the vector-valued function f:R^5-->R^3 below?

f(x,y,L,u,v) = [x + u;
y + v;
L]

and

h(x,y,L) = sqrt((y-x)^2 + (L)^2)
 

Thank you in advance.

 

 

 

First 1144 1145 1146 1147 1148 1149 1150 Last Page 1146 of 2428