status = mkdir('D:\PK7\kpk'); cd D:\PK7\kpk
Fr = .1; M = .2; Kp = 0.50; lambda = 0.1; Nr = 0.1; Pr = 7; Rd = 0.5; Nb = 0.5; Nt = 0.5; H = 0.01;
Ec = 0.01; Le = 2; Sr = 1; D = 0.5; n = 1; E = 0.5; Bi = 0.5; Slip = 0.1;
allM = [0.2 0.3 0.5];
for k = 1 : numel(allM)
M = allM(k); % Get the M for this iteration.
% MINATI's computations:
ODE = @(x,y)[ y(2); y(3); - y(1)*y(3) + (1+Fr)*y(2)^2 + (M+Kp)*y(2) - lambda*( y(4) - Nr*y(6) );
y(5); -(Pr/(1+Rd))*( y(1)*y(5) + Nb*y(5)*y(7) + Nt*y(5)^2 + H*y(4) + Ec*y(3)^2 );
y(7); -Pr*Le*y(1)*y(7) - (Nt/Nb)*(-(Pr/(1+Rd))*(y(1)*y(5) + Nb*y(5)*y(7) + Nt*y(5)^2 + H*y(4) + Ec*y(3)^2)) + Le*Pr*Sr*(1 + D*y(4))^n *y(6)*exp(-E/(1 + D*y(4)))];
BC = @(ya,yb)[ya(1); ya(2)-1-Slip*ya(3); ya(5)+Bi*(1-ya(4)); ya(7)+(Nt/Nb)*ya(5); yb([2;4;6])]; xa = 0; xb = 6; x = linspace(xa,xb,101);
solinit = bvpinit(x,[0 1 0 1 0 1 0]); sol = bvp5c(ODE,BC,solinit); S = deval(sol,x); f0 = deval(sol,0);
% Save the values for this loop
Cf(k) = f0(3); Nu(k) = - (1+Rd)*f0(5); Sh(k) = - f0(7);
fprintf( 'for k = %d, Cf = %f, Nu = %f, Sh = %f.\n', Cf(k), Nu(k), Sh(k) )
end
output = [Cf(:), Nu(:), Sh(:)]
writematrix(output, 'PriyaSRM11.xlsx','Sheet', 'Results');
%% I want to save the output in specific columns of the Excel sheet, say, First column output in A1:A3, Second in B4:B6, and third in C7:C9

2 Kommentare

Dyuman Joshi
Dyuman Joshi am 26 Dez. 2023
It can be done, but it won't be a good idea as you will problems while reading the data.
Any particular reason why you want to do this?
MINATI PATRA
MINATI PATRA am 26 Dez. 2023
Actually this will be helpful for long data

Melden Sie sich an, um zu kommentieren.

 Akzeptierte Antwort

Dyuman Joshi
Dyuman Joshi am 26 Dez. 2023

0 Stimmen

Well, here's a way -
Fr = .1; M = .2; Kp = 0.50; lambda = 0.1; Nr = 0.1; Pr = 7; Rd = 0.5; Nb = 0.5; Nt = 0.5; H = 0.01;
Ec = 0.01; Le = 2; Sr = 1; D = 0.5; n = 1; E = 0.5; Bi = 0.5; Slip = 0.1;
allM = [0.2 0.3 0.5];
for k = 1 : numel(allM)
M = allM(k); % Get the M for this iteration.
% MINATI's computations:
ODE = @(x,y)[ y(2); y(3); - y(1)*y(3) + (1+Fr)*y(2)^2 + (M+Kp)*y(2) - lambda*( y(4) - Nr*y(6) );
y(5); -(Pr/(1+Rd))*( y(1)*y(5) + Nb*y(5)*y(7) + Nt*y(5)^2 + H*y(4) + Ec*y(3)^2 );
y(7); -Pr*Le*y(1)*y(7) - (Nt/Nb)*(-(Pr/(1+Rd))*(y(1)*y(5) + Nb*y(5)*y(7) + Nt*y(5)^2 + H*y(4) + Ec*y(3)^2)) + Le*Pr*Sr*(1 + D*y(4))^n *y(6)*exp(-E/(1 + D*y(4)))];
BC = @(ya,yb)[ya(1); ya(2)-1-Slip*ya(3); ya(5)+Bi*(1-ya(4)); ya(7)+(Nt/Nb)*ya(5); yb([2;4;6])]; xa = 0; xb = 6; x = linspace(xa,xb,101);
solinit = bvpinit(x,[0 1 0 1 0 1 0]); sol = bvp5c(ODE,BC,solinit); S = deval(sol,x); f0 = deval(sol,0);
% Save the values for this loop
Cf(k) = f0(3); Nu(k) = - (1+Rd)*f0(5); Sh(k) = - f0(7);
% fprintf( 'for k = %d, Cf = %f, Nu = %f, Sh = %f.\n', Cf(k), Nu(k), Sh(k) )
end
output = [Cf(:), Nu(:), Sh(:)]
output = 3×3
-1.1266 0.5060 -0.3373 -1.1557 0.5039 -0.3359 -1.2112 0.4997 -0.3331
%Run a for loop through the columns of output
for k=1:size(output,2)
%Generate the range to save the data
str = sprintf('%s%d:%s%d', 64+k, 3*(k-1)+1, 64+k, 3*k)
%Save data via writematrix() specifying range
writematrix(output(:,k), 'PriyaSRM11.xlsx', 'Sheet', 'Results', 'Range', str)
end
str = 'A1:A3'
str = 'B4:B6'
str = 'C7:C9'
%Check the contents of the file
readcell('PriyaSRM11.xlsx', 'Sheet', 'Results')
ans = 9×3 cell array
{[ -1.1266]} {[<missing>]} {[<missing>]} {[ -1.1557]} {[<missing>]} {[<missing>]} {[ -1.2112]} {[<missing>]} {[<missing>]} {[<missing>]} {[ 0.5060]} {[<missing>]} {[<missing>]} {[ 0.5039]} {[<missing>]} {[<missing>]} {[ 0.4997]} {[<missing>]} {[<missing>]} {[<missing>]} {[ -0.3373]} {[<missing>]} {[<missing>]} {[ -0.3359]} {[<missing>]} {[<missing>]} {[ -0.3331]}

5 Kommentare

MINATI PATRA
MINATI PATRA am 27 Dez. 2023
Bearbeitet: MINATI PATRA am 27 Dez. 2023
Wow
Thanks Dyuman for your input.
I tried to change the following line to find the output in A1:A25, B1:B25, C1:C25 ( Serially )
%Generate the range to save the data
allM = linspace(0,0.5,25);
str = sprintf('%s%d:%s%d', 64+k, 25*(k-1)+1, 64+k, 25*k)
but unable. Please have a last try.
Thank you again.
If you want to store the data serially, just call readmatrix() directly, no need to specify range -
writematrix(matrix, filename)
Did you mean to write A1:A25, B26:B50, C51:C75 ?
MINATI PATRA
MINATI PATRA am 27 Dez. 2023
Bearbeitet: MINATI PATRA am 27 Dez. 2023
Yes since there are 03 outputs, I want to put in three different columns A1:A25, B1:B25, C1:C25
Dyuman Joshi
Dyuman Joshi am 27 Dez. 2023
Concatenate the outputs in a single matrix and call the readmatrix() directly.
As I mentioned, there's no need to specify range if you want to store data serially.
MINATI PATRA
MINATI PATRA am 27 Dez. 2023
ok
Thank you Dyuman

Melden Sie sich an, um zu kommentieren.

Weitere Antworten (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by