Question

Generalization needed in dsolve code

Pr = 1; ODE = @(x,y) [y(2); y(3); -(1/2)*y(3)*y(1); y(5); - (Pr/2)*y(1)*y(5);]; BC = @(ya,yb)[ya(1);ya(2);...

5 months ago | 0 answers | 0

Question

functional values and graphs are not matching with this code

Pr = 1; ODE = @(x,y) [y(2); y(3); -(1/2)*y(3)*y(1); y(5); - (Pr/2)*y(1)*y(5);]; BC = @(ya,yb)[ya(1);ya(2...

5 months ago | 0 answers | 0

Question

Can this be arranged to run

function [u,x,y] = poisson(f,g,bx0,bxf,by0,byf,D,Mx,My,tol,MaxIter) x0 = D(1); xf = D(2); y0 = D(3); yf = D(4); dx = (xf - x0)...

5 months ago | 1 answer | 0

Question

Wrong graph is coming with this code

function main A=0.5; pr=1; a=1; phi=0.1;rhos=997.1;Cps=4179;ks=0.613;rhof=8933;Cpf=385;kf=401; a1=(1-phi)^-2.5*((1-phi)+phi...

7 months ago | 1 answer | 0

Question

How to incorporate v and w

%% u_t = D1 * u_xx, u(0,t) = u(1,t) = 0, u(x,0)= cos(pi*(x-0.5)), D1 = 0.1. % % v_t = D2 * v_xx + a*w , v(0,t) = v(1,t) = 0, ...

8 months ago | 0 answers | 0

Question

Can someone rearrange the code to run

%subdivisions space /time ht=0.01; Tmax=1.2; nx=33; hx=1/(nx-1); x=[0:hx:1]'; %matrices K=stiff(1/pi^2,hx,nx); M=mass(1/ht,hx...

8 months ago | 1 answer | 0

Question

How to retrieve the values and substitute to get all the functions

x = sym('x');f(x) = sym('f(x)');g(x) = sym('g(x)'); h(x) = sym('g(x)');a = sym('a');b = sym('b');c= sym('c'); f = x+(1/2).*a....

8 months ago | 1 answer | 0

Question

can it be possible to use odeVectorField

syms f(x) g(x) h(x) Pr Nb Nt Le [V] = odeToVectorField( diff(f,3) - diff(f)^2 + f * diff(f,2) == 0,diff(g,2) + Pr * ( f * diff(...

9 months ago | 1 answer | 0

Question

symbolic code ends with an error

Tr = 0.5; M = 1; Kp = 0.1; K = 02; L = 0.5; D = 3; Rd = 0.5; Pr = 2; S1 = 0.01; Nb = 02; Nt = 1; Le = 1; Kc = 1; ...

10 months ago | 1 answer | 0

Question

My answer is not matching with attached file

syms k r a=sym('a'); b = sym('b');L=sym('L'); M = sym('M'); b1 = sym('b1'); m=7; F = sym(zeros(m,1)); F(1)=0; F(2)=1; ...

1 year ago | 0 answers | 0

Question

How to retrieve the successive values from the given code

syms k r a b k=7;F = zeros(k,1);F(1)=0;F(2)=1;F(3)=a/2; for r = 1:k F(k+3)=(symsum((r+1)*F(r+1)*(k-r+1)*F(k-r+1),r,0,k) -...

1 year ago | 0 answers | 0

Question

how to extract F(k+3) and G(k+2) from symsum code and applying Initial Condition to find series solution

syms x k r f(x) g(x) a b beta b1 M L syms F(k) G(k) % F(0)=0;F(1)=1;F(2)=a/2;G(0)=0;G(1)=1/2;G(2)=b/2;b1=1/beta; %%%%dnf=dif...

1 year ago | 0 answers | 0

Question

Solution of Recurrence relation to find a series expression

syms x k r f(x) g(x) a b beta b1 M L syms F(k) G(k) F(0)=0;F(1)=1;F(2)=a/2;G(0)=0;G(1)=1/2;G(2)=b/2;b1=1/beta; %%%%dnf=diff(...

1 year ago | 0 answers | 1

Question

taylor series expansion with initial condition

syms x f(x) p f1=taylor(f(x),x,'order',3) (D(D(f))(0)=p; D(f)(0)=1; f(0)=0; %%% I want to put initial conditios (D(...

1 year ago | 1 answer | 0

Question

CONVERSION OF ODE TO RECURRENCE RELATION

syms x k r f(x) g(x) a b beta b1 M L syms F(k) G(k) F(0)=0;F(1)=1;F(2)=a/2;G(0)=0;G(1)=1/2;G(2)=b/2;b1=1/beta; %%%%dnf=diff(...

1 year ago | 0 answers | 0

Question

symbolic code with boundary condition

t = sym('t'); x = sym('x'); a=5;a1=2;xa=0;xb=3; g = zeros(1,3,'sym');ga = zeros(1,3,'sym');g(1)=a*x+(5/2)*x.^2; fo...

1 year ago | 0 answers | 0

Question

symbolic and numerical code merged to solve

function main %%%%%%NUMERICAL CODE A=0.5; pr=1; a=1; phi=0.1;a1=2;a2=1;xa=0;xb=6; solinit=bvpinit(linspace(xa,xb,10),[0 1 ...

1 year ago | 2 answers | 0

Question

The symbolic code is not running

syms t x a p q r a1 a2 A pr f(1)=x+p*x^2/2;g(1)=a*x+q*x^2/2;h(1)=1+r*x; for i=1:5 %(Can I take i=0:5) fa(i) = subs(...

1 year ago | 0 answers | 0

Question

Need symbolic calculation, but unable to do

Please follow the attached pdf to consider.

1 year ago | 0 answers | 0

Question

Figs. 1, 2, 3 are not coming while using the following code

Pr=1;Sc=1; R=0.1; xa=0;xb=6; Lv = [-2.5:0.1:0]; for D=[ 0 0.5 1] for i=1:length(Lv) L = Lv(i); fODE =...

1 year ago | 1 answer | 0

Question

How to draw fig.1 with this code

Pr=1;Sc=1; R=0.1; xa=0;xb=6; Lv = [-2.5:0.1:0]; for D=[ 1 1.5 2] for i=1:length(Lv) L = Lv(i); fODE =...

2 years ago | 0 answers | 0

Question

Subscripted assignment dimension mismatch.

function main D=1; %L=0; Pr=1;R=0.1;Sc=1; xa=0;xb=6; Lv = [-2.5:0.025:0]; p = []; for i=1:length(Lv) L = Lv(i); ...

2 years ago | 2 answers | 0

Question

What changes should have been made to run the code

%% HERE is the code format long wt=pi/2; t=1; % t=0.01:0.01:5; % w=pi/2; w=wt./t; R=10;H=2;M=5;Kp=0.5;Gr=5;Gc=5;Kc=...

2 years ago | 0 answers | 0

Question

Needs derivative w.r.t. ' y ' & then graph w.r.t ' t '

y=0:0.01:5; P2=5; %%% t=1; U=exp(y*sqrt(P2)).*(1-erf(y./(2.*sqrt(t))+ sqrt(P2.*t))); S= - diff(U,y); ...

2 years ago | 1 answer | 1

Question

Graph is not coming with this MATHEMATICA code (new to this)

TypeEQ = 2; NumEQ = 3; f[1, z_, {f_, g_, s_}, sigma_] := D[f, {z, 2}] + Gr*Pr*g - Nr*s + sigma; f[2, z_, {f_, g_, s_}, sigma...

2 years ago | 0 answers | 0

Question

How to draw Fig. 4

function main Pr=1; L=-1;D=1; R=0.1;Sc=1; % D=input('D='); %%D=0.5, 1, 1.5 xa=0;xb=6; solinit=bvpinit(linspace(xa,xb,10...

2 years ago | 1 answer | 0

Question

Why this ERROR is coming

M=1;Kp=100; P2=M+(1/Kp);P3=P2+1i*w; U7=exp(y*sqrt(P3)).*(1-erf(y/(2*sqrt(t))+ sqrt(P3*t))); Error using erf (error functi...

2 years ago | 0 answers | 0

Question

How to draw Fig. 1 from the attached pdf with this code

function main Pr=1; G=0.1; % phi=input('phi='); %%0,.05, .1, .15, .2 phi=0.0; rhof=997.1;Cpf=4179;kf=0.613; %for...

2 years ago | 1 answer | 0

Question

can matlab find ilaplace

F1=ilaplace(exp(-y*sqrt(s+a))/s); pretty(F1)

2 years ago | 2 answers | 0

Question

Using LAPLACE TRANSFORM Can it be possible to draw the attached problem

Using LAPLACE TRANSFORM Can it be possible to draw the attached problem.

2 years ago | 0 answers | 0