Issue in plotting of output after function

2 views (last 30 days)
Unable to plot output after function.
I gave my script below
syms theta
A = cos(theta).*[1 5 6; 2 9 3; 5 1 0];
B = tan(2*theta).*[8 5 1; 0 1 6; 2 4 7];
myfun = @(t,y)scriptname(t,y,A,B);
% dummy values for tspan and y0
tspan = [0 1];
y0 = [0; 0; 0];
% ode solver
sol = ode45(myfun,tspan,y0);
function dydt = scriptname(t,y,A,B)
V = [1; 5; 3];
% evaluation of A and B (numerical) with theta = y(3)
An = double(subs(A,y(3)));
Bn = double(subs(B,y(3)));
dydt = V-((An+Bn)*y);
end
for i=1.:3.
stator=y(:,i);
figure(1.);
subplot(2.,2.,i);
plot(t,stator);
xlabel('time');
ylabel('Current');
title(['stat current # ',int2str(i)]);
end
iam getting error like this
Function definitions in a script must appear at the end of the file.
Move all statements after the "scriptname" function definition to before the first local function definition.
How to resolve this??
Thanks in advance

Accepted Answer

KSSV
KSSV on 8 Sep 2021
Edited: KSSV on 8 Sep 2021
function myfun()
syms theta
A = cos(theta).*[1 5 6; 2 9 3; 5 1 0];
B = tan(2*theta).*[8 5 1; 0 1 6; 2 4 7];
myfun = @(t,y)scriptname(t,y,A,B);
% dummy values for tspan and y0
tspan = [0 1];
y0 = [0; 0; 0];
% ode solver
sol = ode45(myfun,tspan,y0);
%% Plot here
plot(sol.x,sol.y)
function dydt = scriptname(t,y,A,B)
V = [1; 5; 3];
% evaluation of A and B (numerical) with theta = y(3)
An = double(subs(A,y(3)));
Bn = double(subs(B,y(3)));
dydt = V-((An+Bn)*y);
Save the aboove code in myfun.m and run it.
  10 Comments
Bathala Teja
Bathala Teja on 8 Sep 2021
sorry to bothering you again, one last doubt.
i followed your 1st method(keeping whole script as function), i am getting all waveforms in one graph. How to plot each solution separatley.
KSSV
KSSV on 8 Sep 2021
figure
plot(sol.x,sol.y(1,:))
figure
plot(sol.x,sol.y(2,:))
figure
plot(sol.x,sol.y(3,:))

Sign in to comment.

More Answers (0)

Categories

Find more on Programming in Help Center and File Exchange

Products

Community Treasure Hunt

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

Start Hunting!