Frage


hold on option is not working in the present code
M = 0.2; L = 0.01; K = 0.1; G = 2; Pr = 2; s = 0.1; Ec = 0.01; m = 0.5; R = 0.5; fw = 0.5; a = 1; n = 0.5; for M = [0 .1 .2] ...

etwa 3 Jahre vor | 1 Antwort | 0

1

Antwort

Frage


Can the code be modified and matched with the attached pdf
%% Can someone modify the code to run and match figs in the attached pdf AND also want to draw some surf plot % u_t(x,t) = ...

mehr als 3 Jahre vor | 0 Antworten | 0

0

Antworten

Frage


mixed odes and pdes in matlab solver
Can it be possible to solve Eqns. (26) - (30), To draw all the figs of attached pdf. Is there any solver? Help me!

mehr als 3 Jahre vor | 0 Antworten | 0

0

Antworten

Frage


Legend option and numbering in boxes
y1 = [85 83 82 87; 81 79 84 86; 82 81 89 87; 89 85 87 80; 90 91 85 87]; figure(1), bar(y1,'stacked'); labels = arrayfun(@(value...

mehr als 3 Jahre vor | 2 Antworten | 0

2

Antworten

Frage


vpa(Fv,5) can not consider all the terms
Fv = - 4.4574e-25*x^21 + 3.0044e-21*x^20 - 5.6303e-17*x^19 + 1.7379e-13*x^18 - 1.4891e-9*x^17 + 1.3636e-8*x^16 - 2.8806e-9*x^15 ...

mehr als 3 Jahre vor | 1 Antwort | 0

1

Antwort

Frage


Which format can be choosen
syms x f = - 0.00000000000000000000000044574112852848068606808271927672*x^21 + 0.0000000000000000030044127779970749176557242324...

mehr als 3 Jahre vor | 1 Antwort | 0

1

Antwort

Frage


Addition of derivative goes wrong
syms M k Ec Pr y T(M, k, Ec, Pr, y) = (A*Ec*Pr/(Ec*Pr-4*A))*( exp(-2*y*sqrt(A)) - exp(-y*sqrt(Ec*Pr)) ) + (y*sqrt(Ec*Pr)/2 ...

mehr als 3 Jahre vor | 0 Antworten | 0

0

Antworten

Frage


Derivative of a function of three variables
%% Can it be possible to find derivative of 'C' simultaneously with all the variables syms y S K C = exp(-y*sqrt(S*K)); CD =...

mehr als 3 Jahre vor | 1 Antwort | 0

1

Antwort

Frage


loop in dsolve code needed
D = 10; R = 0.5; a = 0.1; b = 0.5; A = 1.5; B = 1; C = 1; P = 7; xb = 5; N = 5; syms x f0(x) g0(x) f(k) g(k) eqn0 = [ diff(f...

mehr als 3 Jahre vor | 0 Antworten | 0

0

Antworten

Frage


modification needed in dsolve code
Da = 10; Rd = 0.5; Tw = 1.5; C1 = 1; a1 = 0.1; a2 = 0.5; A1 = 1.5; B1 = 1; Pr =7; syms x f0(x) g0(x) f1(x) g1(x) f2(x) g2(x) f...

mehr als 3 Jahre vor | 1 Antwort | 0

1

Antwort

Frage


dsolve code is not running
syms x f(x) g(x) h(x) M2 M3 a1 k K L A Rd Br Pr Nb Nt G Le H eqn = [ diff(f,2) - L*diff(f)*diff(g) + A*a1*exp(L*g) == 0,... ...

fast 4 Jahre vor | 3 Antworten | 0

3

Antworten

Frage


Need dsolve code to run for Symbolic solution
syms U(y) T(y) C(y) M Kp Pr phi Gr Gc Sc Kc a % % % Gr = 1; Gc = 1; M =1; Kp = 0.1; Pr = 1; phi = 0.1; Sc = 0.22; Kc = 1; ...

etwa 4 Jahre vor | 0 Antworten | 0

0

Antworten

Frage


Can dsolve is enough to solve these ODEs
syms a s t y u(y,t) g(y,t) h(y,t) M Kp Pr phi Gr Gc Sc Kc U(y,s) T(y,s) C(y,s) % % Gr = 1; Gc = 1; M =1; Kp = 0.1; Pr = 1; p...

etwa 4 Jahre vor | 0 Antworten | 0

0

Antworten

Frage


Nu and Sh value should be non-zero
function main x = linspace(0,1,101); t = x; m = 0; sol = pdepe(m,@pde,@ic,@bc,x,t,[]); U = sol(:,:,1); T = sol(:,:,2); ...

etwa 4 Jahre vor | 0 Antworten | 0

0

Antworten

Beantwortet
help solving differential equations
@ eden meirovich What about bvp 4c?

etwa 4 Jahre vor | 0

Frage


Can this dsolve code find solution recursively
syms x f f0(x) f1(x) f2(x) f3(x) f4(x) f5(x) f6(x) eqn0 = diff(f1,3) == 0; cond0 = [f1(0) == 0, subs(diff(f1),0) =...

etwa 4 Jahre vor | 0 Antworten | 0

0

Antworten

Frage


Unable to solve the collocation equations -- a singular Jacobian encountered
b1 = 0.1; b2 = 0.01; A = 0.1; L = 1; b = 0.3; M = 1; R = 2; E = 1; Pr = 2; N = 2; rho = 1; a1 = 1; Ec = 10; G = 1; BC = @(ya,yb...

etwa 4 Jahre vor | 0 Antworten | 0

0

Antworten

Frage


Unable to find value of T1
syms L m x xx s theta_b h Z y U T0 T1 T2 k0 b Tinf Tb sumsol k=6; sym('u0(x)');('u1(x)');('u2(x)');('u3(x)');('u4(x)');('u5(x)...

etwa 5 Jahre vor | 1 Antwort | 0

1

Antwort

Beantwortet
ERROR came as "Index exceeds matrix dimensions".
u, θand ϕ are coded in u, v and w respectively.

mehr als 5 Jahre vor | 0

Beantwortet
ERROR came as "Index exceeds matrix dimensions".
Dear Walter Please see Equation (12), (Attached pdf) 2nd line is for gl 3rd line is for gr

mehr als 5 Jahre vor | 0

Frage


How to get single curve.
H=10;R=5;Pr=1;Q=H-(R/Pr); xl=0; xr=5; J = 10; dx = (xr-xl) / J; tf = 01; Nt = 100; dt = tf/Nt; mu = dt/(dx)^2; % Evaluate the ...

mehr als 5 Jahre vor | 1 Antwort | 0

1

Antwort

Frage


ERROR came as "Index exceeds matrix dimensions".
%% When I incorporate 3 eqns, ERROR came as "Index exceeds matrix dimensions". %%Here I have initial condition f=[0 0 0 ], BC: ...

mehr als 5 Jahre vor | 3 Antworten | 0

3

Antworten

Beantwortet
How to code in bvp4c
I need the to write fw= - 0.2 and 0.2 in the CODE so that the legend will come automatically. Attached fig. can be follow...

mehr als 5 Jahre vor | 0

Frage


How to code in bvp4c
How to code in bvp4c I have drawn 6 curves in a fig. and want 1st, 3rd and 5th curve as dashed lines others solid lines in...

mehr als 5 Jahre vor | 1 Antwort | 0

1

Antwort

Frage


How to plot for: -T0(5) (Y-axis) vs Nt (ranges from 0.1 to 0.5 in X-axis) where already 'x' varies in the differential
function main format('long'); gg=['r','k','b','g','m','c','y','r.','m.','k.']; Pr=10; Le=10; %Nt=.1; Nb=.1; xa=...

mehr als 6 Jahre vor | 0 Antworten | 0

0

Antworten