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

Dyuman Joshi
Dyuman Joshi 2023년 12월 26일
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 2023년 12월 26일
Actually this will be helpful for long data

댓글을 달려면 로그인하십시오.

 채택된 답변

Dyuman Joshi
Dyuman Joshi 2023년 12월 26일

0 개 추천

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

MINATI PATRA
MINATI PATRA 2023년 12월 27일
편집: MINATI PATRA 2023년 12월 27일
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 2023년 12월 27일
편집: MINATI PATRA 2023년 12월 27일
Yes since there are 03 outputs, I want to put in three different columns A1:A25, B1:B25, C1:C25
Dyuman Joshi
Dyuman Joshi 2023년 12월 27일
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 2023년 12월 27일
ok
Thank you Dyuman

댓글을 달려면 로그인하십시오.

추가 답변 (0개)

카테고리

질문:

2023년 12월 26일

댓글:

2023년 12월 27일

Community Treasure Hunt

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

Start Hunting!

Translated by