Feeds
질문
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);...
대략 4년 전 | 답변 수: 0 | 0
0
답변질문
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...
대략 4년 전 | 답변 수: 0 | 0
0
답변질문
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)...
대략 4년 전 | 답변 수: 1 | 0
1
답변질문
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...
4년 초과 전 | 답변 수: 1 | 0
1
답변질문
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, ...
4년 초과 전 | 답변 수: 0 | 0
0
답변질문
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...
4년 초과 전 | 답변 수: 1 | 0
1
답변질문
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....
4년 초과 전 | 답변 수: 1 | 0
1
답변질문
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(...
4년 초과 전 | 답변 수: 1 | 0
1
답변질문
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; ...
4년 초과 전 | 답변 수: 1 | 0
1
답변질문
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; ...
거의 5년 전 | 답변 수: 0 | 0
0
답변질문
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) -...
거의 5년 전 | 답변 수: 0 | 0
0
답변질문
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...
거의 5년 전 | 답변 수: 0 | 0
0
답변질문
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(...
거의 5년 전 | 답변 수: 0 | 1
0
답변질문
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(...
거의 5년 전 | 답변 수: 1 | 0
1
답변질문
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(...
거의 5년 전 | 답변 수: 0 | 0
0
답변질문
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...
거의 5년 전 | 답변 수: 0 | 0
0
답변질문
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 ...
거의 5년 전 | 답변 수: 2 | 0
2
답변질문
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(...
대략 5년 전 | 답변 수: 0 | 0
0
답변질문
Need symbolic calculation, but unable to do
Please follow the attached pdf to consider.
대략 5년 전 | 답변 수: 0 | 0
0
답변질문
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 =...
5년 초과 전 | 답변 수: 1 | 0
1
답변질문
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 =...
5년 초과 전 | 답변 수: 0 | 0
0
답변질문
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); ...
5년 초과 전 | 답변 수: 2 | 0
2
답변질문
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=...
5년 초과 전 | 답변 수: 0 | 0
0
답변질문
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); ...
5년 초과 전 | 답변 수: 1 | 1
1
답변질문
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...
5년 초과 전 | 답변 수: 0 | 0
0
답변질문
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...
5년 초과 전 | 답변 수: 1 | 0
1
답변질문
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...
5년 초과 전 | 답변 수: 0 | 0
0
답변질문
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...
5년 초과 전 | 답변 수: 1 | 0
1
답변질문
Using LAPLACE TRANSFORM Can it be possible to draw the attached problem
Using LAPLACE TRANSFORM Can it be possible to draw the attached problem.
5년 초과 전 | 답변 수: 0 | 0