Answered

Heat transfer differential equation

Call the function this way: [t,y] = ode45(@f,[0,5],1673)

Heat transfer differential equation

Call the function this way: [t,y] = ode45(@f,[0,5],1673)

environ 10 heures ago | 0

Answered

2D information to 1D colorplot

Simple example x = 0:10; y = x.^2; temp = rand(1,10); % temperature cm = jet(10); % colormap ind = 1+rou...

2D information to 1D colorplot

Simple example x = 0:10; y = x.^2; temp = rand(1,10); % temperature cm = jet(10); % colormap ind = 1+rou...

environ 10 heures ago | 0

Answered

Quiver with a system of ODEs.

For vector field u = x'; v = y'; To plot the curve [T,Y] = ode45 ... plot(Y(:,1),Y(:,2)) % (x,y) I suppose x and y depen...

Quiver with a system of ODEs.

For vector field u = x'; v = y'; To plot the curve [T,Y] = ode45 ... plot(Y(:,1),Y(:,2)) % (x,y) I suppose x and y depen...

environ 22 heures ago | 0

Answered

Solving nonLinear equation and Linearized Eqaution using ODE45

YOu have to express , and . Solve the system function du = nonlin(t,y) fi = y(1); psi = y(2); theta = y(3); ...

Solving nonLinear equation and Linearized Eqaution using ODE45

YOu have to express , and . Solve the system function du = nonlin(t,y) fi = y(1); psi = y(2); theta = y(3); ...

environ 22 heures ago | 0

Answered

How to plot magnitude and directions between 2 points?

try this (not tested) data = load ('Input.txt'); [row, col] = size (data); Y2 = data (:,1); X2 = data (:,2); Y1 = data (:,3...

How to plot magnitude and directions between 2 points?

try this (not tested) data = load ('Input.txt'); [row, col] = size (data); Y2 = data (:,1); X2 = data (:,2); Y1 = data (:,3...

1 jour ago | 0

| accepted

Answered

surface fixed point value Z

I want to know the value of Z from the surfc function Here is it

surface fixed point value Z

I want to know the value of Z from the surfc function Here is it

4 jours ago | 0

Answered

Where are the bugs for this ODE finite difference problem that solve using Newton Raphson method?

Here is my attempt clc,clear N = 10; % 10 Set parameters L = 16*0.3048; %meter b = 19.625 * 0.0254; d = 1.625 * 0.0254; ...

Where are the bugs for this ODE finite difference problem that solve using Newton Raphson method?

Here is my attempt clc,clear N = 10; % 10 Set parameters L = 16*0.3048; %meter b = 19.625 * 0.0254; d = 1.625 * 0.0254; ...

8 jours ago | 1

| accepted

Answered

I am receiving the error message Index exceeds the number of array elements (1) when trying to use ode45.

Try @(V,T)

I am receiving the error message Index exceeds the number of array elements (1) when trying to use ode45.

Try @(V,T)

10 jours ago | 0

| accepted

Answered

I have a 14 equation system of second order differential equations which I am trying to solve using dsolve. The program just keeps running without finishing ?

alpha and E should be numerical function main m = 15000; k = 150000; n = 14; A = zeros(n); A(n,n) = -k/m; A(n,n-1) = k/m;...

I have a 14 equation system of second order differential equations which I am trying to solve using dsolve. The program just keeps running without finishing ?

alpha and E should be numerical function main m = 15000; k = 150000; n = 14; A = zeros(n); A(n,n) = -k/m; A(n,n-1) = k/m;...

11 jours ago | 0

Answered

Euler 3D rotation between two vectors

Here is rotation using Rodriguez matrix Here is rotation about Y axis and Z. Here is the difference between them See ...

Euler 3D rotation between two vectors

Here is rotation using Rodriguez matrix Here is rotation about Y axis and Z. Here is the difference between them See ...

11 jours ago | 0

| accepted

Answered

Simultaneous differential equations - derivative of y^3 wrt t

Just get dydt(1)=(-3*D*Cs/rho/r0^2)*y(1)*(1-y(2)); dydt(1)=nthroots(dydt(1),3);

Simultaneous differential equations - derivative of y^3 wrt t

Just get dydt(1)=(-3*D*Cs/rho/r0^2)*y(1)*(1-y(2)); dydt(1)=nthroots(dydt(1),3);

11 jours ago | 0

Answered

Numerically derive a continous, non-symbolic function

Derivative is (if you have numerical data) dy = (y(i)-y(i-1)) / (t(i)-t(i-1)); Maybe diff (if you have a function) ? syms x ...

Numerically derive a continous, non-symbolic function

Derivative is (if you have numerical data) dy = (y(i)-y(i-1)) / (t(i)-t(i-1)); Maybe diff (if you have a function) ? syms x ...

13 jours ago | 0

Answered

3D Point Cloud - Gaussian Curvature

one way clc,clear % generate some data r = 3; t = linspace(0,10)'; x = r*cos(t); y = r*sin(t); z = sin(1*t); t1 = ...

3D Point Cloud - Gaussian Curvature

one way clc,clear % generate some data r = 3; t = linspace(0,10)'; x = r*cos(t); y = r*sin(t); z = sin(1*t); t1 = ...

14 jours ago | 0

Answered

Runge Kutta for system of eqs

I wrote function in the way similar to ode45 y_dot = @(x,dx)[dx; -4*dx-5*x]; %insert function to be solved for i = 1:n k1...

Runge Kutta for system of eqs

I wrote function in the way similar to ode45 y_dot = @(x,dx)[dx; -4*dx-5*x]; %insert function to be solved for i = 1:n k1...

15 jours ago | 0

Answered

The problem of finding intersection point of triangle mesh

To check every edge: for j = 1:size(face,2) i1 = 1; for i2 = [2 3 1] % checking edges [1 2], [2 3], [3 1...

The problem of finding intersection point of triangle mesh

To check every edge: for j = 1:size(face,2) i1 = 1; for i2 = [2 3 1] % checking edges [1 2], [2 3], [3 1...

16 jours ago | 0

| accepted

Answered

How can i solve a system which contains matrices

Here is the best way i know: A = [5 3 1 1 1 1 3 1 2]; B = [c1(:) c2(:) c3(:)]'; C = inv(A) * B;

How can i solve a system which contains matrices

Here is the best way i know: A = [5 3 1 1 1 1 3 1 2]; B = [c1(:) c2(:) c3(:)]'; C = inv(A) * B;

18 jours ago | 1

| accepted

Answered

Using a for-loop to score values in different intervals

Correct form if 0 <= x(i) && x(i) < 1

Using a for-loop to score values in different intervals

Correct form if 0 <= x(i) && x(i) < 1

19 jours ago | 0

| accepted

Answered

Intersection between line and circle when line ends inside circle (using geom2d)

Someone should copy their comment to the answers section so it can be accepted. It should be me % Position of line: x1 = [0...

Intersection between line and circle when line ends inside circle (using geom2d)

Someone should copy their comment to the answers section so it can be accepted. It should be me % Position of line: x1 = [0...

19 jours ago | 1

| accepted

Answered

Removing points from meshgrid (mask)

Try this simple script [X,Y] = meshgrid(1:20); n = [1 -1]; n = n/norm(n); % normal of a line dist = -7; ...

Removing points from meshgrid (mask)

Try this simple script [X,Y] = meshgrid(1:20); n = [1 -1]; n = n/norm(n); % normal of a line dist = -7; ...

20 jours ago | 0

Answered

How can I find the surface normal of a 8 by 3 matrix?

Use griddata to create a surface Use surfnorm to find normal vectors [nx,ny,nz] = surfnorm(x,y,z); surf(x,y,z) hold on quiv...

How can I find the surface normal of a 8 by 3 matrix?

Use griddata to create a surface Use surfnorm to find normal vectors [nx,ny,nz] = surfnorm(x,y,z); surf(x,y,z) hold on quiv...

20 jours ago | 0

Answered

How do I increase the circle size and animate a stick figure?

Use for loop and pause t = linspace(0,2*pi); x = cos(t); y = sin(t); axis([-10 10 -10 10]) hold on for i = 1:10 h = p...

How do I increase the circle size and animate a stick figure?

Use for loop and pause t = linspace(0,2*pi); x = cos(t); y = sin(t); axis([-10 10 -10 10]) hold on for i = 1:10 h = p...

20 jours ago | 0

Answered

Finding specific data point x for given y

I need to find the CB and CS values at t=5 Use interp1 to find i need to find the t and CB values when CS = 0.001 Use fsolve ...

Finding specific data point x for given y

I need to find the CB and CS values at t=5 Use interp1 to find i need to find the t and CB values when CS = 0.001 Use fsolve ...

20 jours ago | 0

Answered

Surface between 3D lines

You are concatenating in a wrong way XX = [D1 D1 D1]; YY = [F1 F2 F3]; ZZ = [z1 z2 z3]; surf(XX,ZZ,YY,'EdgeColor','none')

Surface between 3D lines

You are concatenating in a wrong way XX = [D1 D1 D1]; YY = [F1 F2 F3]; ZZ = [z1 z2 z3]; surf(XX,ZZ,YY,'EdgeColor','none')

20 jours ago | 0

| accepted

Answered

two body problem using ode45

- equation of equilibrium Is this condition fullfield?

two body problem using ode45

- equation of equilibrium Is this condition fullfield?

21 jours ago | 0

Answered

Numerically Integrating the differential equation below?

Any guidance on it is much appreciated! Helpfull page: ODE45 I do not know what I can do if I no tf. Try tf = 5. This should ...

Numerically Integrating the differential equation below?

Any guidance on it is much appreciated! Helpfull page: ODE45 I do not know what I can do if I no tf. Try tf = 5. This should ...

22 jours ago | 0

Answered

Streamline plot not complete

Works ok for me N = 15; xx = linspace(min(x(:)),max(x(:)),10); yy = linspace(min(y(:)),max(y(:)),10); [startx,starty] = mes...

Streamline plot not complete

Works ok for me N = 15; xx = linspace(min(x(:)),max(x(:)),10); yy = linspace(min(y(:)),max(y(:)),10); [startx,starty] = mes...

22 jours ago | 0

| accepted

Answered

2D Matrix revolution

Use cylinder for each curve separately: z = linspace(-1,1,20); % axis r = sqrt(1-z.^2); ...

2D Matrix revolution

Use cylinder for each curve separately: z = linspace(-1,1,20); % axis r = sqrt(1-z.^2); ...

22 jours ago | 1

Answered

How do I align the vertices and the adjacency graph of a Buckyball?

Here is a way ind = full(B); % convert to ordinary matrix [rr,cc] = find(ind); % find non-zero elements ix = ...

How do I align the vertices and the adjacency graph of a Buckyball?

Here is a way ind = full(B); % convert to ordinary matrix [rr,cc] = find(ind); % find non-zero elements ix = ...

24 jours ago | 0

Answered

Find the coordinates of a point chosen on a plot

Short example x = linspace(0,2,20)'; % generate some data y = sin(x); plot(x,y,'.-b') xy = ginput(1); ...

Find the coordinates of a point chosen on a plot

Short example x = linspace(0,2,20)'; % generate some data y = sin(x); plot(x,y,'.-b') xy = ginput(1); ...

26 jours ago | 1

| accepted

Answered

How do I create a sphere without sphere function that will inflate like a balloon utilizing the getframe function?

Use surf to create a 3D surface rr = linspace(0,0.5,20); tt = linspace(0,2*pi,40); [R,T] = meshgrid(rr,tt); [X,Y] = pol2cart...

How do I create a sphere without sphere function that will inflate like a balloon utilizing the getframe function?

Use surf to create a 3D surface rr = linspace(0,0.5,20); tt = linspace(0,2*pi,40); [R,T] = meshgrid(rr,tt); [X,Y] = pol2cart...

28 jours ago | 0