آخرین مطالب
امکانات وب
r = [ 1, cos(t), cos(t)^2, cos(t)^3, cos(t)^4, cos(t)^5, cos(t)^6]
A = [ r; diff(r,t,1); diff(r,t,2); diff(r,t,3); diff(r,t,4); diff(r,t,5); diff(r,t,6)]
W = det(A) سوال چهارم...
برچسب : نویسنده : matlabaachen بازدید : 61
l=input('please input your desired length')
w=input('please input your desired width')
fprintf('your desired length is %.2fn',l)
fprintf('your desired width is %.2fn',w)
area=l*w
fprintf('therefore your desired area will be %.2fn',area)
end سوال چهارم...
برچسب : نویسنده : matlabaachen بازدید : 57
y=exp(x)
plot(x,y,'r')
title('exponential e')
xlabel('x')
ylabel('y')
hold on
S=sin(x)
plot(x,S)
Q=cos(x)
plot(x,Q)
legend('expo','sin(x)','cos(x)') سوال چهارم...
برچسب : نویسنده : matlabaachen بازدید : 50
p=.001
k=.1
equ=(sqrt((p^2)-(4*m*k))/(2*m))==(4*pi*sqrt(-1))
a=solve(equ,m)
a=eval(a)
سوال چهارم...برچسب : نویسنده : matlabaachen بازدید : 46
close all
clc
A = [0 1;-10000 2000]
u=eig(A)
stiffness=max(u)/min(u)
سوال چهارم...برچسب : نویسنده : matlabaachen بازدید : 55
clear all
close all
h = 0.1;
y = 1;
step = 0;
t = 5;
%y'=2*y/t -> t and y increases -> y as a result of t so increase t.
while step <= t
y = y + h*(2*y/(y+0.1))
step = step+h;
end
برچسب : نویسنده : matlabaachen بازدید : 57
clear y % y values from previous runs
a=0 % Initial time
b=1 % Final time
N=10 % Number of time steps
y0=0 % Initial value y(a)
h=(b-a)/N; % Time step
t(1)=a;
y(1)=y0;
for n=1:N % For loop, sets next t,y values
t(n+1)=t(n)+h;
y(n+1)=y(n)+h*f(t(n),y(n)); % Calls the function f(t,y)=dy/dt
end
plot(t,y)
title('Euler Method using N=’,num2str(N),’ steps, by MYNAME’')
% Include your own name سوال چهارم...
برچسب : نویسنده : matlabaachen بازدید : 61
clear all
close all
clc
x_0=0;
y_0=1;
h=0.1;
X=0:0.1:5;
L_X=length(X);
Y=-zeros(1,L_X)
Y(1)=y_0;
y_n=y_0;
for mm=2:L_X
Y(mm)=y_n+h*-2*(y_n);
y_n=Y(mm);
end
plot(X,Y)
xlabel('x')
ylabel('y') سوال چهارم...
برچسب : نویسنده : matlabaachen بازدید : 62
%as an example copy -> [1 2 3; 2 4 6;3 6 9]
%-> paste into command window
matrixvalue(A)
---------------------------------------------------------
function [y] = matrixvalue(A)
[m,n] = size(A);
m
n
if m==n
fprintf('The Matrix is a symmetric so it is a square-matrix.n')
fprintf('The rank of the matrix is:t')
rank(A)
fprintf('The Eigenvalues are:t')
[vektor,value] = eig(A);
diagonal = eig(A);
vektor
value
maxEigenvalue = max(diagonal)
else
fprintf('The Matrix is not symmetric, so the Eigenvalue can not be given.')
end
برچسب : نویسنده : matlabaachen بازدید : 64
برچسب : نویسنده : matlabaachen بازدید : 61