Output of function not variable I assign
古いコメントを表示
When I run this function I want that vector to be called u so I can use it in another script file. Instead it calls it "ans" which I can't use. What is going on? Here is the function, just run it for n=4.
function [u]=finitevolume_KG(n)
%n is the number of volumes
delta_x = 1/n;
A=zeros(n-1);
A(1,1)=1;
A(1,2)=-1;
for i=2:n-2
A(i,i-1)=-1;
A(i,i)=2;
A(i,i+1)=-1;
end
for i=n-1
A(i,i-1)=-1;
A(i,i)=2;
end
f=ones(n-1,1);
f_x=((delta_x)^2)*f;
u = (A)\(f_x);
回答 (2 件)
per isakson
2013 年 3 月 17 日
編集済み: per isakson
2013 年 3 月 18 日
Try the call
u = finitevolume_KG( 4 )
and make it a habit to end your functions with the keyword "end".
2 件のコメント
Jan
2013 年 3 月 18 日
I still omit the trailing "end" to support the small part of Matlab 6.5 users.
per isakson
2013 年 3 月 19 日
Yes, but "end" kind of indicates the end of the workspace of the function.
Jan
2013 年 3 月 18 日
Some simplifications:
function u = finitevolume_KG(n)
%n is the number of volumes
delta_x = 1/n;
A = zeros(n-1, n-1);
A(1,1) = 1;
A(1,2) = -1;
A(2:n-2, 1:n-3) = -1;
A(2:n-2, 2:n-2) = 2;
A(2:n-2, 3:n-1) = -1;
A(n-1, n-2) = -1;
A(n-1, n-1) = 2;
u = A \ (delta_x ^ 2 * ones(n-1, 1));
カテゴリ
ヘルプ センター および File Exchange で Whos についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!